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11.  AMTAACT  (A)AaMmM«200«MnMI 

^"Cesium,  hydrogen  and  oxygen  adsorption  on  beryllium  clusters  are-studied  using 
restricted  Hartree-FocK  (RHF)  calculations  and  ab  initio  relativistic  effective 
core  potentials.  The  clusters  are  taken  as  cylindrical  plugs  from  Be  wafers. 

Cs-,  H-  and  0- to-suos t ra te  internuclear  distances  are  optimized.  For  each  system 
numerous  low-lying  electronic  states  are  investigated  and  Mul liken  electron 
populations  analyzed.  RHF  calculations  show  that  Bcj',  with  three  layers  of 
atoms,  is  too  small  to  adequately  model  the  Be  surface,  while  Be,,,  a  five- 
iavered  system,  and  Be.!-,  a  seven- 1  ayered  system,  are  more  accurate  represen- 
rations  of  the  bulk  metal.  The  emitted  electron  is  clearly  seen  as  vacating  a 
-olecular  orbital  v^hich  is  localized  in  the  surface  layer  of  the  cluster,  thereby 
giving  further  credence  to  the  model.  RHF  calculations  are  completed  for  Pbl  and 
Dil  semi conqijctor  clusters.  Blue  shifts  in  optical  spectra  and  geometry  changes 
are  shown  to  be  due  to  quantum  size  effects.  Scanning  tunneling  microscopy  is 
used  to  investigate  the  nature  of  colloidal  particles  in  the  15  nm  diameter  size 
f'anqe.  Images  snow  a  near  monodispersion  of  small  Au  clusters.  A  model  of  the 
STM  tip  as  a  oolvatomic  crystalline  sijrfarLe_  Is  shown  '  *  ’  — ’  ■ 


1A  SUAilCT  TMAAi  observat  ions  ot  anomalous  long  range  order. 
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Final  Report 


One  area  of  concentration  during  the  first  year  of  the  grant  involved  a 
series  of  ab  initio  self-consistent  field  (SCF)  calculations  on  clusters  of 
beryllium  atoms  together  with  associated  cesium  and  hydrogen  atoms.  They  were 
completed  by  Maria  M.  Marino,  the  doctoral  student  whose  stipend  and  tuition 
were  covered  under  the  conditions  of  this  grant.  An  analysis  of  both  the 
enpefSental  and  theoretical  aspects  of  exposing  Be  metal  surfaces  to  Cs  and  H 

'  ■,  -It' 

^  was 

y 

* 

SCF  calculations  were  carried  out  on  clusters  of  19,  33  and  45  Be  atoms. 
The  clusters  were  defined  as  cylindrical  plugs  from  a  thin  wafer  of  Be  metal; 
an  arrangement  that  is  consistent  with  experimental  measurements  taken  by  Tompa 
and  Seidl.  Adatoms  were  made  to  approach  the  clusters  along  the  three-fold 
symmetry  axis  from  above  and  below  the  cylinder  faces.  This  is  also  consistent 
with  the  experimental  procedure  of  exposing  a  foil  of  Be  to  Cs  an-j  H.  The 

1  car  Hict-anr>o5  ♦'he  Be  atOmS  W^re  held  fi"ar  a*- 

corresponding  to  the  hep  metal  lattice  constants. 

The  calculations  were  carried  out  to  examine  the  nature  of  the  effects  of 
cesium  and  of  hydrogen  adsorption  on  the  work  function  of  Be  metal.  While  it 
IS  well  known  that  cesium  adsorption  appreciably  lowers  the  work  function  of 
many  metals,  details  of  the  mechanism  by  which  it  occurs  is  still  not  well 
understood.  In  the  case  of  hydrogen,  the  interaction  with  a  metal  such  as  Be 
provides  a  good,  simple  model  for  the  process  in  which  we  are  most  interested, 
that  cf  neaative  hydrogen  production  from  surfaces.  Theoretical  studies  of  the 


adsorption  process  can  provide  detailed  information  about  the  interaction  that 
is  not  obtainable  from  experiment.  For  example/  results  related  to  charge 
distributions  and  orbital  structure  at  the  interface  are  calculated. 

A  good  model  of  the  interaction  between  adsorbate  and  substrate  requires 
the  use  of  clusters  large  enough  to  adequately  desccibs  surface, 

while  remaining  small  enough  for  the  application  of  aco«ai^!!^B|Moretical 
methods.  In  the  study  recently  completed  we  used  relativ«l^U|H  clusters  and 
made  use  of  relativistic  pseudopotentials  and  the  high  ^pMtry  of  the 

systems  to  render  the  calculations  tractable.  The  calculations  were  carried 
out  on  a  Cray  X-MP  supercomputer  at  the  NSF's  Pittsburgh  Supercomputing  Center 
and  at  the  Ohio  Supercomputer  Center.  Closed-  and  open-shell  SCF  calculations 
were  completed  for  each  cluster  system  for  a  range  of  low-lying  electronic 
states  and  positive  ion  states. 

Activities  during  the  second  year  of  the  grant  proceeded  essentially  along 
the  directions  described  in  the  original  proposal.  A  series  of  ab  initio 
self-consistent  field  (SCF)  calculations  on  clusters  of  beryllium  atoms 
together  with  associated  cesium,  hydrogen,  and  oxygen  atoms  was  completed.  An 
analysis  of  both  the  experimental  and  theoretical  aspects  of  exposing  Be  metal 
surfaces  to  Cs  and  H  was  published  in  Surface  Science.  Prof.  Milos  Seidl  of 
the  Stevens  Physics  Department  and  his  graduate  student  Gary  Ton^a  were 
co-authors  of  the  article.  In  addition,  studies  detailing  the  results  of  SCF 
calculations  on  the  simultaneous  cesiation  and  hydrogenation,  as  well  as  the 
simultaneous  cesiation  and  oxygenation,  of  the  Be^^  surface  were  initiated. 

The  SCF  calculations  were  carried  out  on  clusters  of  19,  33  and  45  Be 
atoms.  The  clusters  were  defined  as  cylindrical  plugs  from  a  thin  wafer  of  Be 
mof -3 1  ;  an  arrrinoement  that  is  consis<"9nt  wit’h  P"p®rim®n^a1  TT'“as'i’'on(on+’«s  f^v^n 
by  Prof.  Seidl' s  research  group.  Adatoms  were  made  to  approach  the  clusters 
along  the  three-fold  symmetry  axis  from  above  and  below  the  cylinder  faces. 
This  is  also  consistent  with  the  experimental  procedure  of  exposing  a  foil  cf 
Be  to  Cs,  H  and  0.  The  internuclear  distances  between  the  Be  atoms  were  held 
fixed  at  values  corresponding  to  the  hep  metal  lattice  constants. 

The  calculations  were  carried  out  to  examine  the  nature  of  the  effects  of 
oesium,  hydrogen  and  oxygen  adsorption  on  the  work  function  of  Be  metal.  While 
it  IS  well  known  that  cesium  adsorption  appreciably  lowers  the  work  function  of 


«.»*■  •( 


many  metals,  details  of  the  mechanism  by  which  it  occurs  is  still  not  well 
understood.  In  the  case  of  hydrogen,  the  interaction  with  a  metal  such  as  Be 
provides  a  good,  sin?>le  model  for  the  process  in  which  we  are  most  interested, 
that  of  negative  hydrogen  production  from  surfaces.  Theoretical  studies  of  the 


adsorption  process  can  provide  detailed  information  about  tba  interaction  that 
is  not  obtainable  from  experiment.  For  exang>le,  rasttlta to  charge 
distributions  and  orbital  structure  at  the  interface  are  oilHKed. 


The  adsorption  energies  of  Cs  on  Be^^  and  on  Be^^  were  calcttlated  as  21.8 
and  15.1  kcal/mol'atom,  respectively,  whereas  the  adsorption  of  H  on  Be^^  and 
on  Be^^  was  not  predicted  to  occur  at  all.  Optimum  Cs-to-Be-surface  and 
H-to-Be-surface  distances  were  3.70  and  0.86  A,  with  the  latter  value 


corresponding  to  a  local  energy  minimum.  The  decreases  in  the  ionization 
potentials  of  Be23Cs2  relative  to  Be^^  and  of  Be^3CS2  relative  to  Be^^  were  1.5 
and  0.54  ev,  respectively.  This  con^ares  favorably  with  the  experimental  bulk 
work  function  lowering  of  2.3  eV.  Theoretically,  H  is  predicted  to  adsorb  on 
Bej^g  but  not  on  Be^^  nor  on  8645/  consistent  with  the  fact  that  there  is  no 
experimental  observation  of  adsorbed  hydrogen  on  beryllium  metal.  Analysis  of 
the  Mulliken  populations  for  the  19  Be  atom  system  indicates  that  this  cluster 


,  is  too  small  to  model  the  surface  of  bulk  beryllium.  However,  the  population 
^analysis  for  Be^^  gives  strong  indications  that  this  cluster  is  an  adequate 
model  for  surface  interactions.  Simultaneous  oxygenation  and  cesiation  of  Be^^ 
results  in  a  greater  work  function  lowering  than  does  cesiation  alone,  while 


hydrogenation  of  Be45Cs2  results  in  a  work  function  increase.  The  electron 
population  analysis  for  Be^j02Cs2  indicates  that  the  emitted  electron  is 
predominantly  vacating  a  cesium  orbital.  In  Be^jH2Cs2,  however,  the  Be 
cylinder  is  more  involved  in  ionization,  donating  approximately  40%  :f  the 


emitted  charge. 


Two  additional  graduate  students  worked  on  certain  aspects  of  the  project. 
Doctoral  student  Makoto  Sawamura  investigated  quantum  confinement  phenomena  in 
the  context  of  lead  iodide  clusters  at  the  Hartree-Fock  level  of  appro.ximaticn. 

The  work  is  described  in  detail  in  his  doctoral  dissertation  and  in  an 
articles  published  in  The  Journal  of  Physical  Chemistry,  Physical  Review  B  and 
Chemical  Physics  Letters.  Masters  thesis  student  Ping  Wang  completed  a  project 
involving  a  theoretical  study  of  the  cesium  suboxides  and  the  nature  of  the 
large  work  function  lowering  (as  low  as  0.9  eV)  from  the  standpoint  of  ab 
initic  ralculations  using  the  finite  cluster  model.  She  completed  her  degree 


requirements  in  May,  1990  and  the  results  of  her  study  are  being  published  in 
The  Journal  of  Chemical  Physics. 

During  the  third  year  of  the  grant  permission  was  obtained  from  the  AFOSR 
to  use  certain  of  the  grant  funds  to  assist  in  the  purchase  of  a  scanning 
tunneling  microscope  (STM)  .  The  instrument  was  acquired  ud  '<i^4eboratory  was 
set  up.  The  initial  research  was  carried  out  by  a  third  docfeB^^^tudent,  John 
F.  Womelsdorf.  He  completed  benchmark  studies  on  small  ("’ISIlimilameter)  gold 
colloidal  particles.  The  images  were  striking  and  the  work  wan  published  in 
The  Journal  of  Physical  Chemistry.  Following  this  study,  it  was  planned  to 
study  a  series  of  cesium  suboxides  and  attempt  to  measure  their  work  functions. 

This  work,  together  with  a  study  of  the  fundamental  interaction  between  the 
STM  tip  and  substrate,  are  now  underway  and  are  being  supported  by  a  one  year 
e.xtension  from  the  AFOSR.  A  paper  describing  our  findings  on  tip-surface 
interactions  was  published  in  Surface  Science. 

Publications  1-5  and  9-13  (Section  4)  comprise  the  remainder  of  this  Final 
Report.  Copies  of  the  three  student  theses  (Publications  6,  7  and  8,  Section 
4)  were  sent  separate  from  this  Report. 
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IVFICIS  OF  HYMOGEN  AND  OF  CESIUM  ADSORPTION 

ON  A  innniJUM  surface.- 

A  IIBOmiCAL  AND  EXPERIMENTAL  STUDY 


MM.  IIAROIO.  W.C  ERMLER 

Dtfrnmmt  mti  Otuuce/  Engmtehmg,  Sitveta  Insnnae  of  Ttt 

K/tfTWa  USM 


G.S.  TOMPA  *  and  M.  SEIDL 


fkysiet,  Sttoms  Imstauu  of  Tt 


MJtmtL  USA 


Itemed  29  Fehnnry  tSSI:  aoeepled  for  pubhcetioo  19  September  19*8 


Cetem  eud  hydrofa  edeetpiiaa  on  beryllinm  dusters  ooBtaiaiiig  19,  33  and  43  atoms  aie 
smdM  aains  smnieiad  Haniue  Fock  (RHF)  calculations  and  ab  initio  rdaiivistic  effective  cote 
potaniialL  The  le  ctnaian  am  taken  as  cyimdiricai  pluss  tern  a  Be  metal  wafer.  Be-Cs  and  Be-H 
intemutter  dmumem  ate  opiimted,  while  Be-Be  inteniuckar  distanrte  ate  ftonen  at  the  bulk 
BMtal  vaIntB.  For  each  ijraMn,  aaaamN  low-lyina  eiectranic  states  ate  invesiiiatod  The  calcuU' 
none  ate  catiiad  oat  in  Ihs  eonstti  of  an  axpstimental  study  to  determine  the  effectt  of  Cs  and  of 
H  adsosptian  on  the  twaifc  fMMtian  of  Be  metaL  Auper  electron  spectioecopy  and  experimental 
wash  teniiiw  maasuiamaMs  indieate  that  Hj  does  not  adaotb  on  poiyaystalline  Be.  wlule 
phniaemimien  mi  thick  Cs  osetkyer  ■eeeuwsnts  show  a  2J  eV  Inwetiin  in  the  ewrit  functioo 
of  Be  nwai  npon  Cs  atterptitm.  fUlF  calruletinns  indicate  that  Be,*,  with  three  layen  of  atoms, 
is  tan  aael  la  aBthaasaly  saadal  tha  Ba  settee.  Ctealatiens  on  Be3>,  a  fivnlayered  system,  and 


i  H  dost  not  chasnisotb  on  the  surfaces  of  these  dusters, 
on  potential  of  Btf,  by  13  eV.  The  easittad  eiectroa  from 
otUtal  which  is  localitnd  in  the  suttee  layer. 


It  if  well  known  that  adsoqstion  of  cesiuin  on  a  metal  surface  lowers  the 
work  function  of  that  metal  [1,2].  In  particular,  the  reduction  in  the  work 
functioo  of  beryllium  arising  from  cesium  adsorption  has  been  the  focus  of  a 
very  recent  study  [3],  In  the  present  work,  cesium  and  hydrogen  adsorption  on 
a  beryllium  surface  are  studied  using  ab  initio  quantum  mechanical  calcula¬ 
tions.  adsorption  on  polycrystalline  Be  is  also  studied  by  monitoring  the 
work  function  and  the  Be  (104  eV)  differentiated  Auger  electron  signal 
intensity. 


*  Current  address:  EMCORE  Corporation.  35  Elizabeth  Avenue.  Sotnerset.  NJ  08873.  USA 
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M.M.  Marmo  ttaL  /  Hydn/gtn  and etsuan  on  Be  chaten 


Thecmtical  studies  of  the  adsoipuon  process  yield  detailed  information 
about  the  metal  surface  which  is  not  easily  obtained  experimentally.  For 
example,  results  related  to  charge  distributions  and  orbital  structure  are 
Many  such  studies  have  used  clusters  as  models  of  the  surface 
{4-7].  A  good  model  of  the  adsorption  process,  however,  requires  the  use  of 
dusters  large  enough  to  aocuratdy  describe  the  metal  surface  and  its  interac¬ 
tion  with  the  substrate,  while  remaining  small  enough  for  application  of 
aocuraie  theoretical  methods.  In  this  study,  the  dusters  are  taken  as  cyfadrical 
phigs  from  a  Be  wafer  having  Be-Be  intemuclear  distaiMies  equal  Ikl' 
the  bulk  hep  metal,  Larger  spherical  dusters  were  the  focus  ci 
[8].  The  use  of  effective  core  potentials  to  describe  the  core 
and  beryllium  aBow  the  calculations  to  remain  tractable  wliBe 
well-defined  levd  of  theory.  The  complexity  of  the  calodatioos  is 
further  by  preserving  full  point-group  symmetry  b  each  cyliadrical 
duster. 


im  ivicDn0QiQ|y 

2.1.  ThearttictU  model  and  calculations 

The  modd  systems  are  taken  as  cyhndrical  “plugs”  from  Be  wafers  havmg 
surfaces  corresiModtiig  to  the  (0001)  hep  metal  faces.  A  view  of  this  wafer 
along  the  c  directian  is  illustrated  b  fig.  1.  Tabe  1  defines  all  the  cyhnders  of 
a  given  radius  wfaiefa  may  be  formed  fiom  a  wafer  of  a  given  thickness.  In  this 
study,  three  such  cybiders  are  treated.  The  first  has  a  thickness  of  three  layers 
aud  is  comprised  of  19  Be  atoms;  the  second  is  five  layers  thick  and  contains 
33  Be  atoms,  while  the  third  contains  45  atoms  and  seven  layers.  All  three 
cyhpden  have  a  radius  Rj  which  bdudes  a  one-unit  cell  step  along  the  a 
directian.  Be-Be  intemuciw  distances  are  equal  to  those  in  the  bulk  hep 
metal  (a  -  2.29  A,  c  •  3.58  A)  [9].  The  resultmg  cyhnders  have  pomt-group 
symmetry,  and  the  adsorption  of  atomic  hydrogen,  one  on  the  top  sun'ere  and 
one  on  the  bottom  surface  of  all  three  cyhnders.  is  modeled  such  that  the 
three-fold  symmetry  is  preserved.  The  same  apphes  to  the  adsoipuon  of 
atomic  cesium.  In  all  cases,  the  Be-adsorbate  mtemuclear  distances  have  been 
optimized  m  self-consistent  field  (SCF)  calculauons.  The  19-  and  33-atom 
cetiated  and  hydrogenated  cyhnders  are  shown  m  figs.  2  and  3.  The  45-atom 
cyhoder  is  not  shown,  but  it  is  similar  to  the  one  picture  in  fig.  3.  except  that  it 
contains  seven  layers  instead  of  five,  and  the  two  surface  layers  are  identical  to 
the  surface  layers  of  Be, 9  (fig.  2).  Results  obtained  using  these  systems  are 
compared  to  those  calculated  by  treaung  idenucal  bare  Be  cyhnders. 

Calculations  were  accomphshed  on  a  Cray  X-MP  supercomputer  using 
programs  based  on  the  “equal  contribution  theorem"  for  iwo-eiectron  in- 
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Fig.  1.  A  view  of  bulk  Be  meul  along  the  c  direction.  Nuclei  are  situated  at  the  apices  of  the 
heaavons.  as  well  as  at  their  centers.  Layers  tcpiesented  by  dashrs  be  at  c/2  on  the  c  axis  (tee 

table  1.) 


tqtnls  [10].  Ab  initio  restricted  closed-shell  and  restricted  open-shell 
Htrtree-Fock  calculations,  each  corresponding  to  an  average  energy  of  con- 
nguration.  were  carried  out  on  numerous  low-lying  sutes  of  each  cluster.  (The 

'  tbMal 

'#a  ehiftcn  by  ooordiiiatioa  cylinder  and  hep  layer 
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1 /•coordinate  1 

Rq 

«i 

*  0 

aoc/2 

ic/2 

0 

3 

3 

6 

c 

1 1 

1  1 

7  7 

7  7 

c/2 

000 
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1111 
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7  7  7  7 

-  c/2 

000 
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I  I 

I  I 

7  7 

7  7 

-3c/2 

0 

3 

3 

6 

1 

8 

K» 

No.  atoms 

113  3 

1  79  15 

7  13  27  31 

7  19  33  45 
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•vcrate  «aer|y  of  configuntioo  is  defined  as  the  weighted  mean  of  the 
of  all  the  nuihiplets  for  the  configuration.)  The  largest  clusters 
Hguirwi  about  one  hour  of  Cray  X-MP  time  for  each  geometric  orientation. 

foOowing  baau  sets  of  contracted  Gaussian-type  functions  were  used  for 
;  liiyUiani.  cesium  and  hydrogen,  respectively:  (3s2p)/[2slp]  (11],  (5sSp)/(3s2pl 
^'^2).  and  (4slp)/[2slp]  [13].  Ab  initio  effective  potentials  (EP)  (14)  werr  used 
10  represent  tte  Is  core  electrons  on  Be  [11]  and  the  ls-4s.  2p-4p.  and  3d-4d 
core  electrons  in  Cs  [12].  In  the  case  of  Cs.  relauvistic  effects  were  incorpo¬ 
rated  into  the  EP  [14],  Binding  energies  calculated  relative  to  the  completely 
dissociated  clusters  are  given  in  uble  2.  (Valence  SCF  energies  for  atomic  Be. 
Be".  Cs.  Cs",  H  and  H'  are  -0.95083.  -0.65457.  -19.84225.  -19.72876. 
-0.49928  and  -0.44815  hanrees.  respectively,  for  the  basis  sets  used  in  this 
study.)  The  total  adsorption  energies  given  in  table  2  for  the  cesiated  and 
hydrogenated  clusters  were  calculated  relative  lo  the  SCF  energies  of  the 
naked  clusters  plus  that  of  two  adsorbate  atoms. 

Theoretical  work  funcuons  are  given  in  table  3.  The  first  values  in  each  row 
were  calculated  as  the  difference  between  the  total  valence  SCF  energy  of  a 
neutral  cluster  and  that  of  the  ion  generated  by  removing  one  electron  from 
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Table  2 

Restricted  Hartree-Fock  energies  of  suies  of  Be„  and  Be;„X;  clusters 
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Cluster 

Occupied  MO's  ■" 

Cluster 

BE 

(kcsl/mol) 

Total 

adsorption  energy 
(kcal/mol-atomi 

al 

a  2 

a  ^ 

c' 

al' 

Be,,"'’ 

4 

1 

3 

4 

0 

325 

4 

0 

3 

4 

•> 

0 

294.70 

6 

1 

5 

5 

X 

0 

402.99 

38.87 

6 

I 

5 

5 

3 

0 

405.37 

Be„H, 

4 

I 

3 

4 

2 

1 

393.76 

9425  . 

4 

! 

4 

2 

0 

579.66 

Be, 3 

6 

1 

5 

6 

4 

1 

724.94 

lBe„r 

6 

1 

6 

6 

4 

1 

815.76 

Be„Cs2'’*' 

8 

1 

7 

7 

6 

1 

768.46 

21.76  .'V^- 

IBejjCs,]* 

8 

1 

7 

7 

6 

1 

779.74 

7 

1 

5 

6 

4 

1 

HIM 

1.19 

(BenHjl*  •>' 

7 

1 

5 

6 

4 

1 

818.70 

Be.,'*^' 

8 

7 

8 

6 

2 

1078.74 

IBe,,]* 

8 

7 

7 

6 

■> 

1196.77 

Be.,H2 

8 

2 

7 

8 

6 

A 

1115.07 

18.16 

IBe.,H:l* 

8 

*» 

7 

8 

6 

1228.41 

**  Entfies  corre^KNid  to  total  number  of  occupied  MO's  of  each  entry.  Bold  numbers  denotes 


orbital  from  which  the  electron  was  removed. 

Sute  correipoods  to  the  weighted  average  of  configuration. 
“  Open'thell  state:  (a'^l'fa; )'. 

Open-shell  sute:  (a'jl'fa,")' 

*'  Open-shell  sute:  fe")’. 

'*  Opan-sheU  stau:(  (aD'Cai' )'. 

**  Open-shell  stau:  fa'jlVe')'. 

Open-shell  sute;  (ajlVe')'. 


the  highest  occupied  orbital  of  that  cluster.  The  values  immediately  following 
are  due  to  Koopmans’  theorem  and  correspond  to  the  negative  of  the  energy 
of  the  molecular  orbital  (MO)  from  which  the  electron  was  lomzed.  Expen- 
mental  values  obtained  as  discussed  below  are  also  shown. 

Atomic  net  electron  charge  values  arc  given  in  tables  4-6.  as  calculated 
usmg  a  MuUiken  population  analysis  |15].  Values  of  R  appearing  in  these 
tables  correspond  to  optimum  atom-to-surface  distances.  These  distances  were 
calculated  from  the  minimum  of  a  curve  fit  to  SCF  energies  for  three  or  more 
distances. 

2.2.  Experimental  procedure 

Expenmental  measurements  were  performed  in  an  ultrahigh  vacuum  system 
with  a  base  pressure  of  5  x  10"”  Torr.  Pressures  were  monitored  with  a  nude 
lomzation  gauge  whose  controller  allowed  calibration  for  H;.  The  general 
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T*bic3 

fuBcoom of  and  B.,X;  clusien 


♦  «eV)  *• 


BbmH] 

••jj 

BojjCa: 

■•mH, 

*4.5 

Befisipit 

BcCtfeapil 

Mifoipt) 

A4(i.-»iCaI^‘' 


5.35.  5.04 
3.41. 2.99 
4.30.  3.96 
4.41.4.11 
2.92. 2.60 
4.37. 4.09 
3.16.  2.94 
3.40.  3.15 
3.92 
1.6 
3.92 

1.94.2.05 
1.05. 1.08 
1.49. 1.51 
0.05. 0.01 
-0.24.  -0.21 
2.3 
0 


*'  WImk  two  vafaMi  appear,  theie  are  Koopmans'  theorem  and  vaiuea.  retpeciively. 

At  optiaana  Cs  rewarapr 

“  Vahiet  roriaapoad  to  a  work  fuaetioB  lowering. 

Valaaa  conaipoBd  lo  a  work  fttaction  increase. 


expennHOttl  appintus  was  described  previously  (3.16].  The  sample  was 
mouated  on  a  carousel  that  allowed  rotation  to  either  a  retarding  fleld  electron 
diode  station  (17.18],  where  work  function  shifu  were  measured,  or  lo  a 
.,;^2cyliadrKal  ■niror  uudyaer  tutkm  having  a  center>mounted  electron  gun  with 
^whidi  Aiifar  deciron  spectroacopy  (AES)  measurements  were  made.  A  sputter 
i  gun  was  oo-focused  onto  the  sample  with  the  AES  system  and  was  used  to 
stter  clean  the  surface. 

AES  showed  diat  no  O  or  C  contaminants  were  present  after  sputter 
cleanif^  The  only  contaminants  present  were  Ar  ( <  1.1%)  and  N  (  <  2.4%). 
_  ^  and  these  were  due  to  the  sputter  cleaning.  After  cleaning.  AES  and  work 
function  shifts  indicated  the  ^  surface  to  be  free  of  contammants  for  a  penod 
in  excess  of  one  hour. 

The  expenmental  procedure  was  as  follows.  The  sample  was  first  sputter 
cleaned.  H .  could  be  introduced  into  the  system  at  any  time  before,  dunng  or 
after  sputter  cleaning.  Cs  was  deposited  onto  the  surface  from  a  low-energy  Cs 
(S  eV)  source  (3)  before  or  after  sputter  cleaning  at  a  separate  deposition 
station  adjacent  to  the  work  function  station.  The  Auger  signal  was  generally 
monitored  during  cleaning.  After  cleaning,  the  Auger  signal  either  continued 
to  be  monitored  or  the  sample  was  rotated  from  its  position  facing  the 
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AES/^tter  deaning  sution  to  a  new  position  facing  the  work  function 
diode,  thus  altowing  the  work  function  shift  to  be  monitored.  The  reposition¬ 
ing  was  accomplished  in  approximately  twenty  seconds. 


Tottl  adaocpttoo  f  rgiri  for  the  cesiated  and  the  hydrog— ted 
appear  ia  tabk  2.  As  was  previously  mentioned,  the  energies  awe 
rdiutve  lo  the  total  SCF  energies  of  the  naked  clusters  phu  that  M 
adsorbate  atoms  ia  their  ground  sutes.  These  energies  are  38.9. 
and  18J  kcal/mol  •  atom  for  Be„Cs2<  Be}3Cs2,  Be,,H2.  BejiM 
respacdaply.  Pot  dissodatiem  to  Bef,  +  Cs^  +  Cs  reqtiilri5S.fflDeal/ 

mol  •  atom.  PiMOCiatioB  ot  Be,9H2  to  yield  661^9  -t-  H  -f  H~  and  Ber«  -i-  H  +  H" 
requiem  101.6  and  ITS  J  kcal/mol  •  atom,  respectively.  Dissociation  of  Be^iCs. 
to  Be3'^4-Cs'f  Cs''  needs  49.4  kcal/mol  ■  atom.  Energies  of  149.8.  and  64.6 
keak/mai •  atom  are  naadad  to  dissociate  Be,,H2  to  Befj  +  H  +  H"  and  Bei% 
-»■  H  -f  H~,  respectively.  Finally,  the  dissociation  of  Be^jH;  to  yield  BeTs  +  H 
-f-  H~  Md  lais  4-  H  H''  requires  68.1  and  150.12  kcal/mol  •  atom,  respec¬ 
tively.  Therefore,  the  lowest*cnergy  dissociation  limit  is  of  the  form  Be,  +  2X. 
where  X  is  an  adatom  in  the  ground  sute.  This  procedure  will  be  named 
method  1.  There  are  two  ahemative  ways  of  calculating  the  adsorption  energy, 
methods  2  and  3.  Method  2  involves  calculaticm  of  the  adsorption  energy 
reiaitve  to  the  total  enec|y  of  the  lowest  electronic  sute  and  MO  configuration 
for  the  equilibrium  geoasetry  of  the  cluster-adsorbate  system  and  the  same 
rlertrimsT  sute  for  a  duster-adsorbate  separation  of  approximately  10  A. 
^taken  as  mtm»e  sqiaraticm.  This  method  assumes  no  crossmgs  among 


Be,9a,.  BertO,,  Be|9H2  and  Be45H2,  this  assumption  leads  to 
4|lM|Me  erron.  In  these  cases,  t^  lowest  eneigy  electronic  state  corresponding 
'■WWt  calculated  equilibrium  geometry  is  not  the  lowest  energy  sute  calculated 
ii  infinite  cluster-adsorbate  separation.  'This  method  yields  adsorption  en- 
meffm  of  53.7.  56.3.  136.4  and  177.0  kcal/mol  -  atom  for  Be,9Cs->.  Be^^Cs^. 
Be,9H2  and  Be«3H2.  respectively.  Clearly,  method  2  overestimates  the  adsorp¬ 
tion  energies  of  Be, 902-  BctXs:.  Be,9H2  and  Be4,H2  (by  14.8.  23.8.  102.1 
and  158.8  kcal/mol  -  atom,  respectively,  relative  to  method  1).  No  electronic 
sute  crossings  were  found  for  Be^H,.  However,  the  ground  state  of  Be,iH.  is 
open-shell  and  nearly  degenerate  (0.01  eV  lower  in  energy)  with  a  closed-shell 
sute  that  is  3.65  eV  higher  in  energy  at  dissociation  than  the  ground  state. 

Method  3.  the  second  alternative  method  for  calculation  of  the  adsorption 
energy,  involves  taking  the  energy  difference  between  the  lowest  electronic 
state  of  the  cluster-adsorbate  system  at  equilibrium  and  the  lowest  electronic 
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sute  of  the  same  system  at  a  cluster-adsorbate  separation  of  approximately  10 
A.  This  method  yielded  adsorption  energies  of  43.8.  36.8.  99.0.  47.1  and  64.2 
kcal/mol  ■  atom  for  Be^oCs^.  Be^tCs^.  Bei^H^.  Be^^H^  and  Be45H;;.  respec¬ 
tively.  Relative  to  method  1.  this  method  overestimates  the  adsorption  energy 
of  Be,9Cs;.  BC33CS2.  Be^qH;.  Be3jH2  and  Be45H2  by  4.9.  15.0.  64.7.  45.9  and 
46.0  kcal/mol  ■  atom,  respectively.  No  electronic  state  crossings  were  found 
for  Be33H2.  Consequently,  methods  2  and  3  yield  the  same  value  of  the 
adsorption  energy  for  this  cluster.  Ideally,  methods  1  and  3  should  agree, 
provided  that  the  lowest  energy  electronic  state  at 
geometry  and  at  infinite  cluster-adsorbate  separation  has 

There  are  several  problems  inherent  in  using  method  i. 
calculations  presented  here  all  involve  restricted  Hartree-Fodc  lhea^fl|tt||H 
is  not  generally  accurate  for  the  description  of  systems  tuKvipf  JOOinetfMli^V 
from  equilibrium.  As  was  previously  suted.  method  3  requirta  calcufaaion  of 
the  total  energy  for  the  system  at  a  cluster-adsorbate  separation  of  approxi¬ 
mately  10  A.  Inspection  of  the  gross  atomic  populations  (GAPl  for  Be,„Cs.  at 
this  distance  reveals  improper  dissoaauon  to  neutral  atoms.  The  GAP  of 
cesium  for  this  system  should  be  18.  9  electrons  per  Cs  atom.  In  this  case,  the 
GAP  is  17.  with  approximately  1/2  electron  from  each  cesium  6s  orbital  being 
donated  to  the  cluster.  There  should  be  no  such  ionization  at  this  cluster-ad¬ 
sorbate  separation,  unless  an  electron  is  transferred  from  Cs  to  Be  upon 
dissociation,  yielding  Be,',  and  Cs-t-Cs'.  However,  dissociation  to  these 
species  requires  40.1  kcal/mole  more  energy  than  dissociation  to  Be,,  -f  2  Cs. 

RHF  calculations  predict  a  net  eneigy  gam  of  18.2  and  1.2  kcal/mol  upon 
hydrogen  approach  to  the  8043  and  8633  surfaces,  respectively.  The  approach 
of  hydrogen  on  the  Be,,  surface,  on  the  other  hand,  results  in  an  energy  gam 
of  34.3  kcal/mol.  Cesium  adsorption  on  both  the  Be,,  and  Be33  surfaces  yields 
energy  gains  of  38.9  kcal/mol  atom  and  21.8  kcal/mol  ■  atom,  respectively. 

,  Therefore,  while  cesium  is  predicted  to  adsorb  on  both  Be,,  and  Be,,  (i.e.  at  a 
hole  and  at  a  head-on  site)  at  the  SCF  level  of  theory,  hydrogen  adsorbs 
strongly  on  Be,,  and  appears  to  adsorb  moderately  on  864,  (i.e..  at  a  hole  sue) 
«  but  only  negligibly  on  Be,,  (i.e..  at  a  head-on  site). 

Table  3  contains  calculated  values  for  the  ionization  potentials  of  all  the 
clusters  studied,  as  well  as  values  for  the  shift  in  ionization  potential  (A<{>) 
resulting  from  Cs  or  H  adsorption.  Experimentally  determined  values  of  the 
work  functions  due  to  Cs  or  H-  adsorption  on  Be  metal  are  also  included. 
Results  for  the  864,  system  indicate  a  slight  nse  by  0.2  eV  in  the  ionization 
potential  of  the  cluster  due  to  H  adsorption,  while  results  for  the  Be,,  system 
indicate  a  lowenng  by  1.5  eV  in  the  ionization  potential  of  the  cluster  due  10 
Cs  adsorption,  and  no  lowenng  due  to  H  adsorption.  The  expcnmenial 
findings  are  2.3  and  0.0  eV  lowenng  in  the  work  function  of  Be  metal  due  to 
adsorption  of  Cs  and  of  H3.  respectively.  Adsorption  of  Cs  and  of  H  on  Be,,, 
however,  results  in  ionization  potential  lowenngs  of  approximately  2.0  and  1.0 
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Table  4 

Election  populinort  of  Bei,.  Be„C&,  and  Be,,!!; 

Qnsicr  i<«aoidi>  Atom  No.  Net  char^  per  atom 

Mte*'  label  aw«»'»  clutter  [auaterj  ' 


Total  net  charge 
ditfcfCBce 

((Qutter) '  -  Cluster) 


Bcia 


BoiaCt; 


0 

B0O 

1 

1.43 

1.41 

e/2 

6 

-0.08 

-0.07 

e/2 

BcC 

6 

-0.09 

0.00 

0 

BcB 

6 

-0.07 

0.00 

0 

BaO 

1 

1.36 

1.35 

e/2 

BcA 

6 

-0.05 

-0.07 

c/2 

BcC 

6 

-0.17 

-0.11 

0 

Bd 

6 

-0.14 

-0.07 

B*' 

a 

2 

0.39 

0.55 

0 

BcO 

1 

1.38 

1.33 

e/2 

BcA 

6 

-0.16 

-0.15 

e/2 

BcC 

6 

-0.06 

0.07 

0 

BcB 

6 

-0.02 

0.01 

K*' 

H 

2 

0.04 

0.05 

*'  r-3^A. 

Number  of  symmetry  eqnivalant  atoms. 

"  B  —  3.70  A  for  Be  plane  to  Cs  thstaitcc.  (Be  to  Cs  distance  is  3.93  A.) 
B  •  0.8S  A  for  Be  plane  to  H  distance.  (Be  to  H  /ti«»«iw»  u  1.57  A.) 


eV,  ropoctivciy.  AiUioufh  the  first  value  is  in  close  agreement  with  expeh- 
niMt.  the  aetsaad  is  not. 

ThMC  discrepancies  are  explained  using  the  electron  populations  listed  in 
4  and  S.  In  the  case  of  Be,,,  total  net  charge  differences  indicate 
in  charge  redistribution  by  both  the  surfwe  (BcC.  r-coordinate 
SMl  middle  (BeB.  r-ooordinate  —  0)  layers  upon  ionization  of  the 
la  Be,,,  ioaizatiaa  of  the  cluster  results  in  a  contribution  of  0.4 
I  from  k  aumis  (BeB)  in  the  middle  layer  compared  to  O.S  electron  for 
,pi  a|kmis  on  the  surface.  Thus,  the  middle  and  surface  layers  of  the  cluster  arc 
Imitfved  almost  equally  in  the  electron  ionization  process.  For  Be^pCs;.  both 
the  surface  (BeC)  and  middle  fBeB)  layers  contribute  0.4  electron,  while  Cs 
donates  0.3  elearon.  Since  electron  emission  is  from  the  surface  of  the  bulk 
metal,  it  is  concluded  that  Be,p  is  too  small  a  cluster  to  model  the  Be  metal 
surface  involved  in  Cs  adsorption. 

On  the  other  hand,  the  middle  layer  of  BcipH;  contributes  0.2  electron, 
with  the  surface  layers  (BeC)  being  predominantly  involved  in  ionization.  H 
donates  only  a  minimal  amount  (0.02)  of  charge.  Part  of  the  charge  contrib¬ 
uted  by  the  middle  (BeB)  layer  is  shifted  toward  the  BeC  group  of  atoms  in 
Be, pH..  Since  the  greatest  charge  perturbation  occurs  on  the  surface,  this 
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Tables 

Electron  populations  of  Bcsi. 


BeiiCs-)  and  BciiH-i 


Cluster 


Be, 


Be„H, 


Be3,Cs; 


z-coordi- 

Mte*' 

Atom 

label 

No. 

atoms 

Net  charge  per  atom 

Ciusi«  |Clutter| ' 

Total  net 
difference 
((Ouster) 

0 

BeO 

1 

o.% 

0.97 

0.01 

c/2 

BeA 

6 

0.55 

0.53 

-0.14 

c/2 

BeC 

6 

-0.10 

-0.07 

0.14 

0 

BeB 

6 

-0.32 

-0.29 

0.21 

C 

BeH 

12 

-0.18 

-0.10 

0J9 

c 

BeD 

2 

0.20 

0.14 

-0.11 

0 

BeO 

1 

0.94 

0.93 

-OjOI*'" 

c/2 

BeA 

6 

0.57 

0.52 

-025 

c/2 

BeC 

6 

-0.12 

-0.07 

027 

0 

BeB 

6 

-0.34 

-0.33 

Oi0» 

c 

BeH 

12 

-0.28 

-0.22 

0.67 

c 

BeD 

0.44 

0.40 

-0.08 

Cs 

‘y 

0.46 

0.60 

0.29 

0 

BeO 

\ 

1.00 

1.00 

0.00 

c/2 

BeA 

6 

0.45 

0.44 

-0.12 

c/2 

BeC 

6 

-0.06 

-004 

0.12 

0 

BeB 

6 

-0.30 

-0.26 

0.24 

c 

BeH 

12 

-0.23 

-0.16 

0.84 

c 

BeO 

0.87 

0.82 

-0.10 

H 

2 

-0.24 

-0.23 

-0.02 

c-3.58  A. 

Number  of  symmetry  equivalent  atoms. 

“  X  *  3.77  A  for  Be  to  Cs  distance- 

/I  •  1.38  A  for  Be  to  H  distance  at  local  minimum  (see  text). 


model  appears  to  be  adequate  for  the  description  of  H  adsorption.  However, 
further  analysis  reveals  that  this  is  not  the  case.  As  was  stated  previously,  the 
middle  layer  (BeB)  of  Be,,  participates  appreciably  (donating  0.4  electron)  in 
the  ionization  process,  indicating  that  this  cluster  does  not  reasonably  model 
the  bulk  surface.  Adsorption  of  Cs  or  H  onto  this  surface,  therefore,  yields 
information  about  the  interaction  of  these  two  species  with  a  Be  cluster,  but 
this  information  cannot  be  interpreted  as  also  applying  to  the  process  wherebv 
a  Cs  or  H  atom  adsorbs  onto  the  bulk  surface.  Such  an  extension  may  be 
made  only  in  cases  where  the  bare  cluster  is  a  good  model  of  the  bulk  metal. 

Be,,,  a  five-layer  cluster,  and  Be45,  a  seven-layer  cluster,  appear  to  model 
the  Be  metal  surface  more  appropriately.  Ionization  of  these  clusters  indicates 
that  the  surface  layers  (BeH  for  Bcjj  and  BeG  for  Be«,)  donate  most  of  the 
emitted  electron.  Although  the  middle  layers  also  contribute  charge,  the 
involvement  of  the  surface  layers  is  three  to  four  times  as  great.  The  adsorp- 
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the  Be  surface  awL  in  fact,  does  not  in  the  case  of  Be,]  (the  adsorption  energy 
is  1  kcai/nol'atoin).  The  adsorption  energy  of  H  on  Be*;,  although  greater 
(18  kcal/mol  ■  atom),  is  attributed  to  a  local  minimum  since  both  the  popula- 
non  analysis  and  A4  values  indicate  no  effect  on  the  Be*,  surface  due  to  H 

adsocptiOB. 

H  epfMaadMS  the  Be,,  surface  directly  above  a  Be  atom  (directly  overhead 
site)  in  the  priaeat  calculations.  There  are  three  other  poshhie  sites  for 
■demiiiinn  TiHee  are  the  Be-Be  midpoint,  eclipsed  asd;^ 
has  llte  adeochaie  approaching  the  Be  surface  directly 
focaiad  on  the  surface.  The  second  involves  adsorptiodl|[A 
trianfle  of  Be  atoms  located  on  the  surface,  the  center  ^ 
shove  a  Be  atom  sitoated  in  a  layer  next  to  the  surface  (sWiicii  it' 

Bgp,}.  The  open  site  has  the  adsorbate  approaching  the 'Be  aorfaCe  (flioctly 
ahovc  the  oeaam  of  a  triangle  of  Be  atoms  located  on  the  surface,  directly 
bdow  whid)  thMO  is  no  Be  atom.  According  to  a  study  by  Bagus  et  al..  the 
dinctiy  overhand  site  »  the  least  stable  for  the  adsorption  of  H  onto  a  Be 
surface  (l^l-  Hmm  adsorption  energies  were  obtained  in  SCF  calculations 
uainf  dmten  nmtainini  two  to  three  layers.  Be,,  is  a  five^layer  cluster  and 
thus  should  yield  a  reasonable  value  for  the  adsorption  energy  associated  with 
the  directly  overhead  approach  of  an  atom  to  the  surface,  while  adsorption  on 
an  edipsed  site  is  a^uatdy  modded  by  the  seven>layer  Be*,  cylinder. 
Method  3  yields 47.1  kcal/mol  •  atom,  while  method  1  predicts  only  negligible 
adsorption  for  H  on  Be,,.  For  reasons  already  suted.  method  1  is  expected  to 
he  more  accuraie.  Therefore.  H  is  predicted  to  ads^  very  weakly,  if  at  all  on 
Be,,.  The  lowen  energy  distance  is  1.58  A.  which  represents  a  local 

mininHim  since,  as  discussed  above.  H  does  not  bind  to  Be,,.  Bagus  et  al 
reponad  an  SCF  adsorption  energy  of  32  kcal/mol  for  H  on  Be^  (a  three-layer 
dMter  with  14  atoms  on  the  top  and  bottom  layers  and  8  atoms  in  the  middle) 
and  a  Be^  -H  of  1.38  A  (19].  These  values  also  correspond  to 

adsorption  on  a  directly  overhead  site. 

AdKxption  of  H  on  Be4,  ammnts  to  18.2  kcal/mol  ■  atom.  The  optimum 
H-to-Be-plane  is  0.86  A.  which  is  also  a  local  minimum  since  the 

electron  populations  of  Be*,  and  are  nearly  idenucal.  indicating  that  H 

has  no  effect  on  -  i.e..  docs  not  chemisorb  to  -  the  Be,,  surface  The 
corresponding  values  for  adsorpuon  onto  a  Be,,,  cluster  reported  by  Bagus  et 
al.  are  42.3  kcal/mol  for  the  adsorpuon  energy  and  0.95  A  for  the  distance 
from  H  to  the  surface. 

It  is  also  the  surface  layers  (BeH)  of  Be,,Cs,  that  are  predomjnanii\ 
involved  in  electron  removal.  These  layers  contribute  0.7  electron,  whereas  the 
greatest  charge  contributed  by  the  inner  layers  is  0.3.  The  total  net  charge 
differences  of  Be„Cs,  differ  from  those  of  Be,,,  indicating  that  Cs.  unbke  H 
IS  affecting  electron  removal.  In  fact.  Cs  donates  0  3  electron  to  this  process 
The  effects  of  charge  redistribution  resulting  from  lomzation  are  greater  fvsr 
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the  inner  layers  (relative  to  the  surface  layers)  of  Be^iCs^  than  for  those  of 
Be3i3H2.  The  BeC  group  of  atoms  contribute  0.3  electron,  while  the  BeA  atoms 
accept  0.3  electron  in  Be^jCs^.  These  amounts  correspond  to  less  than  half  of 
the  charge  contributed  by  the  surface  layers.  In  BcjjH^.  the  greatest  contribu¬ 
tion  by  the  inner  layers  is  approximately  one  fourth  of  the  charge  contributed 
by  the  surface.  The  small  but  non-negligible  involvement  of  the  inner  layers  of 
Be33Cs2  in  the  ionization  process  indicates  that  the  discrepancy  between  the 
calculated  and  measured  work  function  lowering  (uble  3)  may  be 
treating  a  cluster  corresponding  to  a  larger  cylinder  height:  for  exaggp^^^ 
atom  cylinder  of  uble  1.  Another  possibility  is  the  treatment  dT  IpIpifl 
having  the  same  height  as  Be33  or  Bets  but  a  greater  radius  (see  uble 
1 ).  in  addition  to  improvements  in  Ao.  adsorpuon  of  Cs  on  an  rrlipend  laeji 
(for  example,  Be4S)  is  predicted  to  be  more  suble.  Unlike  the  hendwi  ike; 
this  position  allows  the  Cs  to  interact  more  closely  with  the  surface  triangles  of 
Be  since  there  is  no  center  atom  present. 

It  is  noted  that  the  present  study  contains  certain  constraints  and  assump¬ 
tions.  First,  no  geometry  relaxation  was  attempted  for  the  Be  clusters.  Al¬ 
though  the  Be-adatom  distances  were  optimized,  the  full  D31,  symmetry  of  the 
cluster  was  maintained,  and  the  Be- Be  distances  were  held  fixed  at  the 
experimental  lattice  constants  of  Be  metal.  For  clusters  of  this  size,  it  is 
probable  that  the  optimum  geometry  may  not  match  precisely  that  of  the  bulk 
metal.  Surface  geometries  do  not  generally  match  the  corresponding  bulk 
lattice  constants.  Most  materials  are  known  to  expenence  some  degree  of 
reconstruction  at  the  surface  due  to  the  bending  of  outward-directed  (from  the 
surface)  unpaired  orbitals.  Second,  adsorption  of  Cs  and  of  H  was  modeled 
using  only  one  atom  above  the  top  and  bottom  layers  of  the  cluster.  Since  the 
shift  in  the  work  function  of  Be  metal  depends  on  the  degree  of  surface 
coverage  (3].  the  results  may  be  altered  by  adsorption  of  more  than  one  atom 
of  Cs  or  H  on  each  cluster. 

;  Third,  the  calculations  do  not  include  electron  correlation  effects.  As 

i  previously  mentioned,  all  low-lying  electronic  states  of  each  cluster  were 
studied.  Several  of  these  states  were  nearly  degenerate  (within  01  eVi.  making 
a  reordering  of  states  possible  upon  inclusion  of  electron  correlation.  Such 
effects  may  be  accounted  for.  at  least  in  part,  through  the  use  of  local-dcnsitv- 
functional  approaches.  These  methods  replace  the  exchange  terms  in  the 
Hartrec-FcKk  equations  by  sphencal  potentials,  otherwise  known  as  muffin-tin 
potentials,  through  which  some  electron  correlation  corrections  arc  introduced 
into  the  svstcm  Metal  surfaces  such  as  that  of  tungsten  have  been  studied 
using  these  approaches  [20).  However,  the  results  presented  here  are  ab  initio. 
unlike  those  obtained  using  densitv-functional  theors  Furthermore,  correla¬ 
tion  energy  cancellation  is  expected  between  the  neutral  and  cationic  cluster 
svstems  .Although  this  cancellation  is  not  complete  because  the  neutral  svstem 
contains  a  greater  number  of  electrons,  the  clusters  treated  here  are  large 
enough  (e  g..  Be,,Cs.  contains  84  electronsi  so  as  to  minimize  this  error 


Fi|  4.  (1)  The  work  function  shift  of  sputter  cleaned  Be  in  time  exposed  lo  the  residual  gas  in  ihe 
chamber  (7  6x  10""  Torr).  (b)  The  work  function  shift  of  spuiier  cleaned  Be  exposed  to  H 
(;0xl0'"  Tom  at  /  -SO  s.  to  The  change  in  spuiier  cleaned  Be  AES  signal  in  nme  upon 
exposure  lo  ihe  residual  gas  (7  ^  »  m"  "  Tom  in  the  chamber  idark  curvet  and  upon  exposure  n 

H  .  1 3  0  X  1(1  ■"  T  orr  i 
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The  states  were  chosen  by  analyzing  the  highest  occupied  and  lowest  unoc¬ 
cupied  molecular  orbitals.  Finally,  independent  SCF  calculations  were  used  to 
determine  the  ionization  potenuals  of  the  clusters  in  addition  to  the  use  of 
Koopmans'  themem. 

The  experimental  findings  of  Cs  adsorption  on  Be  have  been  presented  and 
discussed  elsewhere  (3].  Exposure  of  the  surface  to  yielded  no  observable 
change.  That  is.  the  Be  substrate  Auger  signal  was  not  attenuated  and  there 
was  mo  disceniable  wotk  function  shift.  Typical  results  are  shown  in  fig.  4.  ip. 
curve  a.  the  change  in  the  work  funcuon  with  time  is  shown  for  tfMtdtr'- 
cleaned  polycrystalline  Be  exposed  to  the  residual  in  the 
7.6  X  10~  "  Torr.  Curve  b  corresponds  to  the  same  sample  with  the  ant^lpS 
that  at  time  /  -  SO  s.  was  let  into  the  chamber  at  2.0  x  10'“  Torr. 
corresponds  to  a  flux  of  3  x  10'^  molecules/cm'  -  s  at  room  tantperantre.  No 
change  in  the  work  function  was  seen  for  a  period  in  excess  of  10^  s. 
corresponding  to  a  fhience  greater  than  3  x  10'“  molecules/cm*.  Curve  c 
shows  the  Be  (104  eV)  Auger  electron  signal  for  a  sample  first  sputter  cleaned 
and  exposed  to  the  residual  gases  in  the  vacuum  (solid  symbol  curve)  and  for 
the  same  sample  sputter  cleaned  in  and  maintained  in  H.  at  3  x  10'“  Torr. 


Analysis  of  the  MuUiken  populations  for  Be,,  indicates  that  this  cluster  is 
too  small  to  model  the  bulk  surface.  However,  the  population  analyses  for  Be,, 
and  Be;t5  indicate  that  these  clusters  are  adequate  models  of  the  Be  metal 
surface.  The  adsorption  energy  of  Cs  on  Be, 3  is  calculated  as  21.8  kcal/  mo)  - 
aipm.  The  decrease  in  the  ionization  potential  of  BejtCs,  relative  to  Be,,  is 
'  dl||culated  to  be  1.S  eV.  The  expenmental  work  function  lowering  resulting 
jjife  Cs  adsorption  is  measured  as  2.3  eV  [3].  The  adsorption  energies  of  H  on 
and  Be«5  are  1  and  18  kcal /  mol  ■  atom,  respectively.  The  first  value  is 
nagligible.  while  the  latter  is  attnbuted  to  a  local  minimum  since  Mulliken 
population  analysis  and  dp  values  for  both  Be,,  and  Be,,  indicate  no  effect  on 
these  clusters  due  to  H  adsorption.  Expenment  reveals  no  adsorption  of  H  -  on 
Be  metal,  thus  ruling  out  dissociative  chemisorption 
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Bi«lii  IS  studied  using  ab  iniuo  Haruee-Fock  calculations  including  relativistic  eflects.  EJearonic  spears,  geometnc  parame- 
ten  and  na  aimme  cbargea  am  reported.  Biiie^ifu  in  the  opucal  absorption  spectra  relative  to  bulk  Bil.  are  observed.  Increases 
in  mterlaycr  rtiinnrri  result  in  atianenad  Bi-Bi  separations.  Resulu  are  consistent  with  quantum  size  efTeas  in  small  semicon- 
duaor  dusien  and  with  acaaaing  umaeiiiiB  mtcroacope  images. 


The  propenies  of  dusters  are  attracting  increased 
attention  because  they  represent  an  intermediate 
physical  regime  where  neither  atomic  nor  solid-sute 
descriptions  are  independeotly  appropriate  [1-3]- 
Although  numerous  experimental  techniques  to  syn¬ 
thesize  clusters  both  in  gaseous  and  condensed  phases 
are  in  use  (4-9].  there  is  no  generally  available  ex- 
pcriiaeatti  RMpdure  to  determine  the  struaure  and 
speqn  «f  llwijfcuura  produced.  Hence,  theoretical 
insiaa  liii|l4ilner  structure  as  a  function  of  cluster 
sizelafcei  odiraMjor  role. 

The  conipiiiiioiial  difficulty  inherent  in  treating 
large  syateim  comprised  of  heavy  elements,  how¬ 
ever.  has  precluded  rigorous  ab  initio  studies  on  such 
clusters.  Thaoretical  treatments  on  systems  of  this 
type  have  largely  negleaed  relativistic  effects,  as  well 
as  certain  exchange  interaaions  among  electrons. 
Consequently,  the  majonty  of  ab  initio  work  in 
semiconductor  systems  has  been  limited  to  small 
clusters  of  light  elements  such  as  silicon  [  lO.I  I  ].  In 
such  strongly  covalent  systems,  it  is  found  that  small 
isolated  clusters  tend  to  reconstruct  to  form  tightly 
bound  symmetrical  structures  with  a  high  degree  of 
coordination.  Moreover,  the  tendency  to  minimize 
the  number  of  exposed  dangling  bonds  means  that 
the  bulk  silicon  structure  is  recovered  only  for  vert 
large  clusters. 

The  structure  of  the  layered  semiconductor  inves¬ 
tigated  here  exhibits  a  radically  different  depen¬ 


dence  on  cluster  size  than  do  those  of  isotropic  sil¬ 
icon.  Unlike  clusters  of  Si.  Bitl,i  essentially  reuins 
the  structure  of  the  bulk  crystalline  lattice  in  the  lat¬ 
eral  direaion  for  the  range  of  geometries  studied.  We 
find,  however,  that  there  can  be  significant  elec¬ 
tronic  coupling  between  interlayer  perturbations  and 
intralayen  lattice  dimensions  in  this  cluster,  indeed, 
for  large  increases  in  the  Bi-1  interlayer  spacinp.  we 
see  in-plane  lattice  consunts  decrease  by  as  much  as 
10%.  With  these  results,  it  might  be  possible  to  ex¬ 
plain  recent  compelling  scanning  tunneling  images 
of  Bil,  clusters  [12]. 

The  Bi»l,«  cluster  model  studied  is  shown  in  fig. 
I .  It  has  the  layer  symmetry  of  the  bulk  crystal  and 
IS  compnsed  of  a  plane  of  bismuth  atoms,  above  and 
below  which  is  situated  a  triangular  layer  of  iodine 
atoms.  This  model  was  chosen  based  on  expenmen- 
tal  data  on  small  clusters  of  the  layered  semicon¬ 
ductors  PbU  and  Bil,  grown  as  colloidal  suspen¬ 
sions.  Clusters  of  the  type  Bil,  were  found  to  be 
compnsed  of  two  hexagonally  closed-packed  layers 
of  iodines  between  which  was  sandwiched  a  layer  of 
metal  atoms  in  a  honeycombed  arrangement  [8] 
Bi«I,g  has  D,e  point  group  symmetry  (cubic  close¬ 
packing  >  and  represents  a  fragment  of  the  bulk  cry  s¬ 
talline  solid  lone  unit  cell  thick)  having  iodine  lay¬ 
ers  slightly  penurbed  from  the  bulk  arrangement  to 
preserve  the  symmetry  of  the  cluster.  The  coordi¬ 
nation  numbers  for  our  cluster  and  for  the  bulk  crys- 
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C, 

Fit.  I  BiJii  clMMr  baviat  Dm  point  graap  symmetry.  (Tbe 
threefold  rouuin  uis  u  indicatad  by  C|. )  Tbrce  tymmetrydis- 
iinct  I  atom  ai«  wdicnted.  All  intcntonuc  distances  correspond 
10  bvik  laoieeeaMttna  The  positKins  of  the  ID  atemts.  however, 
are  lotated  by  49.2*  relative  to  the  bulk. 


Table  I 

Coordinabon  numbers 


Atom 

I 

Bi 

Total 

bulk  •' 

Bi 

6 

0 

6 

I 

6 

2 

8 

cluster 

Bi 

6 

0 

6 

lA 

5 

2 

7 

IC 

2 

4 

ID 

3 

1 

5 

Ref,  (131  •’Seefit.lll 


tal  arc  given  in  table  I ,  It  is  interesting  to  note  that 
these  values  arc  equal  for  Bi.  This  results  from  the 
fact  that  our  model  cluster  was  constructed  such  that 
each  bismuth  lies  at  the  center  of  a  distorted  octa¬ 
hedron.  ( Bi-ID  distances  are  38.5%  longer  than  Bi- 
I  distances  in  crystalline  Bil,. )  This  arrangement  is 


similar  to  that  of  bulk  Bih.  where  each  bismuth  lies 
at  the  center  of  a  nearly  perfect  octahedron  of  io¬ 
dines.  Bi-Bi.  l-I  and  Bi-1  intemuclear  distances  were 
initially  determined  fiom  the  lattice  constants  of 
crystalline  Bil  j  1131.  However,  the  ID  group  of  at¬ 
oms  (see  fig.  1 )  were  routed  by  49.2°  counterclock¬ 
wise  with  respea  to  theonpaal  Mjypgninaeinent  to 
obtain  a  cluster  of  DMiQjpuaetfji^ijpBilWltt  in  IC- 
ID  separations  which  are  the  cor¬ 

responding  bulk  cr>ml  valiMnL||^^BiMidear  dis¬ 
tances  were  then  opitimiMl  in  W MN|Ptamee-Fock 
calculations,  consaaicai  with  the  va^Mon  principle 
and  preserving  the  cylindrical  structure  of  the  clus¬ 
ter.  All  initial  and  optimum  nearest-neighbor  inter¬ 
atomic  distances  in  BiJ,,  arc  given  in  ublc  2.  This 
table  also  contains  the  corresponding  bulk  crystal 
values,  where  apphcable. 

Ah  initio  effective  core  potentials  which  include 
relativistic  effects  (14)  were  used  to  represent  the 
1  s-4s.  2p-4p  and  3d-4d  core  electrons  in  1  (15)  and 
the  Is- 5s.  2p-5p.  3d-5d  and  4f  core  electrons  in  Bi 
(16).  Basis  sets  of  contracted  CaLssian-type  func¬ 
tions  were  used  to  represent  the  valence  orbitals  of 
Bi  and  I  (15.16).  Self-consistent  field  (SCF)  ener¬ 
gies  for  the  lowest-lying  molecular  orbital  electron 
configurations  were  calculated  for  each  geometnc 
orienution  of  the  cluster.  These  ab  initio  I'estiicted 
closed-shell  and  restneted  open-shell  Hanree-Fock 
linear  combination  of  atomic  orbiul  -  molecular  or¬ 
bital  ( LCAO-MO )  calculations  were  accomplished 
using  programs  based  on  the  "equal  contnbution 
theorem”  for  two-elearon  symmetry  orbital  inte¬ 
grals  (17).  The  lowest  unoccupied  MOs  ( LUMOs  I 
of  all  the  electron  configurations  studied  for  neutral 
Bi«l|«  were  quite  low-iying.  For  example,  the  six 


Table  2 

nearesi-neighbor  interamon  disunces  ( A ) 


Atom  pair 

Bulk*' 

Initial 

Optimum 

Bl-IA 

3,09 

3.09 

3.20 

Bi-IC 

3  09 

3  09 

3.20 

Bi-ID 

- 

4.20 

4.28 

lA-IA 

4.37 

4  37 

4.35 

lA-IC 

4.37 

4.37 

4.35 

I.A-ID 

- 

4.37 

4  35 

IC-ID 

- 

6  18 

6.15 

Bi-Bi 

4.37 

4  37 

4.35 

•’  Ref-  (13) 
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LUMOs  for  the  lowest-energy  eiearon  configuration 
studied  had  energies  below  —  2.7S  eV.  indicating  that 
the  neutral  system  is  eiearon  deficient. 

This  conclusion  is  supptmed  by  Koopmans'  theo¬ 
rem  [  18].  If  the  cluster  wavefunaion  1P(2N).  a  2A'- 
elearon  antisymmetnzed  product  of  MOs  {o„|.  is 
sutionary  and  •«  is  an  eifenfunction  of  the  Fock  op¬ 
erator.  then  V(2N-^  1 ).  where  an  eiearon  is  added 
to  Pm-  is  also  sutionary  with  respea  to  any  further 
variations  in  the  orbitals  to  first  order.  This  indi¬ 
cates  that  systems  possessing  very  low-lying  LUMOs 
are  eiearon  deficienu  in  that  addition  of  one  or  more 
elearons  to  these  orbiuls  will  increase  the  overall 
subility. 

To  satisfy  this  eiearon  deficiency,  we  studied  a 
doubly  charged  anionic  BuIk  duster.  We  also  be¬ 
lieve  that  the  anion  is  the  most  natural  cluster  to 
consider  from  the  experimental  point  of  view  since 
clusters  of  this  type  are  prepared  in  solution  as  col¬ 
loidal  particles  in  the  presence  of  excess  iodide  ion 
[19].  The  total  energy  of  Bi«Ifr  was  found  to  be 
more  than  10  eV  lower  than  the  corresponding  neu¬ 
tral  spedes.  In  this  case,  nearly  all  the  LUMOs  have 
energies  greater  than  zero,  indicating  that  the  sys¬ 
tem's  electrophilicity  has  been  quenched.  A  geome¬ 
try  optimization  of  Bi«lii  indicated  that  the  equi¬ 
librium  interlayer  distances  were  expanded  vertically 
m  by  9.Sli«dative  to  bulk  Bil3.  This  vertical  expansion 
•  3.6%  increase  in  Bi-IA  and  Bi*IC  dis- 
ytawei  .jnd  to  a  1.9%  increase  in  Bi-ID  interatomic 
'I  sepnrglibnt  ( uMe  2 ).  Intralayer  distances,  however. 
Wife  contracted  by  only  0.5%  relative  to  those  in  the 
bulk.  This  system  tends  to  contract  laterally  as  the 
interlayer  separations  are  increased.  For  example, 
expansion  of  the  interlayer  separation  by  60%  resulis 
m  a  10%  contraaion  in  metal  layer  lateral  dimen¬ 
sions. 

A  total  of  sixteen  low-lying  MO  electron  config¬ 
urations  of  the  cluster  were  studied  at  both  the  cal¬ 
culated  equilibnum  interatomic  disunces  and  at  a 
geometry  corresponding  to  a  60%  expansion  in  the 
interlayer  distance  and  a  10%  contraction  in  lateral 
dimensions.  These  percentages  were  chosen  from  the 
results  of  a  series  of  calculations  at  geometnes  rang¬ 
ing  from  10%  contraction  to  60%  expansion  of  the 
lattice  constants  in  the  vertical  direction  and  20% 
contraction  to  10%  expansion  in  the  lateral  direc¬ 
tion  All  of  these  configurations,  with  the  exception 


of  states  1 3.  IS.  16.  17.  19  and  20.  possessed  two  un¬ 
paired  elearons  and  were  treated  as  pure  spin  triplet 
sutes.  i.c.  the  cluster  wavefunaions  are  spin  eigen¬ 
functions  having  total  spin  quantum  number  5=1. 
Sutes  13.17  and  20  correspond  to  the  same  electron 
configuration  comprised  of  four  open-shell  electrons 
coupled  to  form  fmt  same  applies 

to  sutes  IS. 

Since  the  heavy 'eiealMMig^H|piiivolve  appre¬ 
ciable  relativistic  aflbcta  |^^^■Min-orbit  split¬ 
ting.  a  precise  for  an  MO 

configuration  iaf  j|.^wmfBH|^ble  3)  is  some¬ 
what  ambiguous,  ihie  calctnmlb  reported  here  are 
based  on  L-S  atomic  coupling  in  the  context  of  av¬ 
eraged  relativistic  effeaive  core  potentials  ( REPs ) 
[14].  More  extensive  calculations  that  make  us  of 
REP-based  spin-orbit  operators  and  explicit  elec¬ 
tron  correlation  are  required  to  specify  the  elearonic 
sutes  more  precisely  (14).  However,  the  present 
calculations  provide  an  accurate  representation  of 
cluster  electronic  sutes  in  terms  of  L-S  coupled  sep¬ 
arated  atomic  sutes.  Note  that  for  large  clusters  the 
dependence  of  elearonic  properties  on  specific  total 
orbital  and  spin  angular  momentum  coupling  is  re¬ 
duced  considerably  due  to  the  presence  of  bands  of 
adjacent  sutes  [20]. 

Exciution  energies  are  given  in  ubie  3.  Each  value 
is  the  difference  between  the  total  SCF  energy  of  the 
ground  sute  and  that  of  an  excited  sute.  State  1 1  at 
1.28  eV  is  the  only  closed-shell  sute  lying  within  4 
cV  of  the  ground  sute.  Table  4  reports  eiearon  charge 
per  atom  values  based  on  a  Mulliken  population 
analysis  [  2 1  ]  for  the  lowest-lying  state  at  each  of  the 
aforementioned  geometnes.  The  energy  separations 
for  equilibnum  Bi*!;,-  listed  in  table  3  are  band-like. 
Sutes  I  through  12  exhibit  a  nearly  monotonic  nse 
in  energy  and  lie  considerably  below  the  expenmen- 
lal  band  gap  of  bulk  crystalline  Bilj  at  2.1  eV  [22]. 
Sutes  1 3  through  20  are.  however,  considerably  blue- 
shifted  relative  to  this  band  gap.  This  may  be  attnb- 
uied  to  quantum  confinement  effects,  as  suggested 
by  recent  Pbl>  cluster  studies  [8.23].  A  total  of  65 
single  eiearon  transitions  among  1 90  possible  tran¬ 
sitions  among  the  low-lymg  tnplet  states  were  ex¬ 
amined.  Fig.  2  shows  that  these  transitions  peak  at 
0.3.  2.9  and  3.7  eV.  The  higher-lying  peaks  are  blue- 
shified  by  0.8  and  1.6  eV  with  respect  to  the  energy 
band  gap  of  the  bulk  crystal.  On  the  other  hand,  the 
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Tables 

EIcctnMiiciiMMafBitIft  (eV) 


SuteNo.  Tamiyaibol  MOoaoiiiiintian*’ 


ila,. 

IOsh 

13e, 

14e, 

1 

% 

2 

0 

0 

4 

1 

0 

2 

*A*i 

2 

0 

0 

4 

2 

0  ' 

3 

’A,„*A,^>E,*' 

2 

0 

0 

4 

1 

4 

% 

2 

1 

0 

4 

1 

5 

2 

0 

0 

4 

0 

t'  ■ 

6 

% 

2 

0 

0 

4 

0 

1 

7 

% 

2 

0 

1 

4 

1 

.. 

8 

‘A,. 

2 

1 

0 

4 

0^ 

9 

% 

2 

0 

0 

4 

1 

0 

10 

’A^ 

2 

0 

0 

4 

0 

II 

‘A,. 

2 

0 

0 

4 

0 

0 

12 

’As. 

2 

0 

1 

4 

0 

13 

•E. 

I 

0 

0 

4 

2 

0 

14 

% 

1 

0 

0 

4 

3 

0 

IS 

% 

1 

0 

0 

4 

1 

16 

^A„ 

1 

0 

0 

4 

0 

17 

’Ah 

1 

0 

0 

4 

2 

1 

18 

*An„*AH.’E.** 

2 

0 

0 

3 

3 

19 

’A,. 

1 

0 

0 

4 

1 

20 

‘Ah 

1 

0 

0 

4 

2 

0 

“  OcnipMiBMaf ilw  auMmamriOtwtitow.  *>Seetataadtabie2. 
AvwnBafU— iripHlMlM.  *'L««Ml4yiiigckiM»I<«heUtiaftet(^ 


A£(eqml.)'” 


to  two  or  more  dectron  exciutions 
md  at  lower  eaeiiies  (approximately 


^aiom  values  appearing  in  ubie  4  in¬ 
predominant  interlayer  interactions 
occur  between  the  atoms  comprising  the  Bi  hexagon 
and  the  laifm  I  triangle  (BiB  and  ID  in  fig.  1 )  for 
the  lowest  energy  sute  of  Bi«Ui-  at  both  geometries 
treated.  This  was  also  observed  for  the  other  fifteen 
electron  configurations  of  Bi«If,-  listed  in  table  3. 
No  differences  were  seen  between  charge/aiom  val- 


TaMe  4 

Population  analysis  for  the  lowest  eneify  states  of  BLl},"  *’ 


Atom 

Charge/aiom  ‘ ' 
(equil.  1 

Charge/atom  •’ 
(expanded  Z) 

Bi 

0.89 

0.67 

lA 

-0.26 

-0.18 

1C 

-0.29 

-0.32 

ID 

-0.67 

-0.51 

Transition  Energy  ttV) 


Ref  (21 1.  “'Seefig.1. 


”  Electron  units. 


Fig.  2.  Distribution  of  transition  energies  of  Bislfr .  The  arro» 
marks  the  position  of  the  bulk  band  gap. 
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ues  for  the  ground  sute  and  the  corresponding  val¬ 
ues  for  sutes  2  through  12  at  the  equilibrium  ge¬ 
ometry.  indicating  that  the  electronic  transitions  are 
predominantly  intralayer  in  nature  and  occur  among 
symmetry  related  atoms.  A  more  detailed  analysis  of 
the  Mulliken  populations  revealed  that  these  elec¬ 
tron  exciutions  correspond  to  metal>metal  transi¬ 
tions.  Differences  in  chai^e/atom  values  between 
states  1 3  through  20  of  Bi«l|f  shown  in  uble  3  ( also 
relative  to  the  equilibrium  geometry )  and  those  of 
the  ground  sute  are  as  large  as  0.07.  which  compares 
to  0.02  for  sutes  2- 1 2.  In  these  cases,  the  transitions 
involve  both  Bi  and  I  atoms. 

We  expect  that  a  less  restriaed  geometry  optimi¬ 
zation.  i.e.  one  wnere  tne  Bi  and  1  layers  are  opti¬ 
mized  independently,  will  lead  to  a  cluster  having 
even  shorter  meul-metal  separations.  Such  was  the 
case  for  Pb7li4  .  a  Uyered  semiconductor  very  sim¬ 
ilar  to  Bil).  Calculations  on  this  represenutive  frag¬ 
ment  of  the  bulk  Pbl;  crysul  indicate  a  10%  reduc¬ 
tion  in  metal-meui  separauons  relative  to  crystalline 
Pbl..  while  1-1  disunces  essentially  remained  at  the 
bulk  lattice  consunts  (23). 

Our  findings  regarding  the  reduaion  of  in-plane 
lattice  consunts  when  interlayer  distances  are  ex¬ 
panded  from  their  equilibrium  values  aid  the  un- 
dersunding  of  recent  scanning  tunneling  microscope 
( STM )  inates  of  layered  Bil3  clusters  ( 1 2  ] .  In  that 
stutfy.  atomiarily  resolved  images  showed  structures 
ha>%g  honeycombed  symmetry  after  colloidal  sus¬ 
pensions  of  Bilj  were  evaporated  onto  graphite  sur¬ 
faces.  Since  the  central  bismuth  plane  in  Bil  3  has  ho¬ 
neycombed  symmetry  ( in  contrast  to  the  hexagonal 
s>  mmetry  of  the  iodine  planes ).  the  STM  images  are 
most  likely  individual  bismuth  atoms.  Interestingly, 
the  lattice  consunt  of  these  puuiive  B110I30  clusters 
was  measured  to  be  0.65  nm.  roughly  20%  smaller 
than  the  bulk  in-plane  lattice  consunt.  which  1$  con¬ 
sistent  with  the  present  calculations.  Hence,  these 
STM  images  could  be  revealing  clusters  whose  top¬ 
most  layer  of  iodine  has  been  removed  or  displaced 
trom  us  underlying  bismuth  plane.  Chemical  oxi¬ 
dation  or  cluster-tip  interactions  could  be  among  the 
phenomena  responsible  for  these  large  perturbations. 
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ABSTRACT 

Cesium,  hydrogen  and  oxygen  adsorption  on  beryllium 
restricted  Hartree-Fock  (RHF)  calculations  and  ab  initio  ^latli 
potentials.  The  clusters  are  taken  as  cylindrical  plugs  from  Be  wafers.  Cs*,  H-  and 
0-to~substrate  internuclear  distances  are  optimized.  For  each  system  numerous 
low-lying  electronic  states  are  investigated,  and  the  Mulliken  electron  populations 
are  analyzed.  The  calculations  are  carried  out  in  the  context  of  an  experimental 
study  to  determine  the  effects  of  various  adsorbates  on  the  work  function  of  the 
substrate.  Auger  electron  spectroscopy  and  experimental  work  function  measurements 
indicate  that  H.  does  not  adsorb  on  polycrystalline  Be,  while  photoemission  and  thick 
Cs  overlayer  measurements  show  a  2.3  eV  lowering  in  the  work  function  of  Be  metal 
upon  Cs  adsorption.  The  continuous  oxidation  of  Cs  has  been  studied  using 
ultraviolet  photoelectron  spectroscopy  and  electron  spectroscopy  by  deexcitation  of 
metastable  noble  gas  atoms.  Results  indicate  that  the  work  function  of  Cs  is  lowered 
upon  exposure  of  the  surface  to  small  doses  of  oxygen.  RHF  calculations  show  that  a 
19  atom  Be  cluster,  with  three  layers  of  atoms,  is  too  small  to  adequately  model  the 
Be  surface,  while  the  33  atom  cluster,  a  five-layered  system,  and  the  45  atom 
cluster,  a ,se0Mi-layered  system,  are  more  accurate  representations  of  the  bulk  metal. 
The  emitt^d^iUctron  is  clearly  seen  as  vacating  a  molecular  orbital  which  is 
localiz^  14^111  surface  layer  of  the  cluster,  thereby  giving  further  credence  to  the 
model. 


I 


1 .  INTRODUCTION 


Adsorption  of  cesium  on  a  metal  surface  lowers  the  work  function  of  that  metal 
by  amounts  generally  in  the  neighborhood  of  two  electron  volts (1].  In  addition,  the 
work  function  of  cesium  itself  is  depressed  by  about  0.5  eV  upon  selective  exposure 
to  0.  to  form  CSjjOj  [2).  In  particular,  the  reduction  in  the  work  function  of 
beryllium  arising  from  cesium  adsorption  has  been  the  focus  of  recent  studies(3]. 
Cesium,  hydrogen  and  oxygen  adsorption  on  a  beryllium  surface  are  studied  here  using 
ab  initio  quantum  mechanical  procedures.  Such  theoretical  studies  of  the  adsorption 
process  yield  detailed  information  about  the  metal  surface  which  is  neither 
straightforwardly  nor  accurately  obtained  from  experimental  data.  For  example, 
charge  distributions  and  orbital  structure  result  from  ab  initio  calculations.  Many 
such  studies  have  used  clusters  as  models  of  the  surface  ( 4-7 )  .  A  good  model  of 
adsorption  processes,  however,  requires  the  use  of  clusters  large  enough  to 
accurately  describe  the  surface  interactions  as  well  as  those  involving  substrate, 
while  remaining  small  enough  for  application  of  accurate  theoretical  methods. 


In  this  study,  the  model  systems  are  defined  as  cylindrical  "plugs"  from  Be 
wafers  having  surfaces  corresponding  to  the  (0001)  hep  metal  faces.  A  view  of  this 
wafer  along  the  c-direction  is  illustrated  in  Figure  1.  All  cylinders  of  a  given 
radius  which  may  be  formed  from  a  wafer  of  a  given  thickness  are  defined  in  Table  1. 

Three  such  cylinders  are  treated  here.  The  first  has  a  thickness  of  three  layers 
and  is  comprised  of  19  Be  atoms;  the  second  is  five  layers  thick  and  contains  33  Be 

♦Proceedinas  of  the  SPIE  OE/LASE  '89  Conference  on  Microwave  and  Particle  Beam  Sources 
and  directed  Energy  Concepts,  January  16-20,  1989,  Los  Angeles,  CA. 


Figure  1.  A  view  of  Bulk  Be  metal  along  the  c  direction.  Nuclei  are 
situated  at  the  apices  of  the  hexagons,  as  well  as  at  their  centers. 
Layers  represented  by  dashes  lie  at  c/2  on  the  c  axis  (see  Table  1) . 


Table  1.  Be  Clusters  by  Coordination  Cylinder  and  hep  Layer 


Cylinder 
Height 
(2-coord. ) 

P, 

n 

P 

'1 

X 

R 

L 

R 

3 

•^/ 2 

3c/2 

0 

3 

3 

6 

1 

1 

1 

1 

7 

/ 

7 

c/2 

0 

0 

0 

3 

3 

3 

3 

3 

3 

6 

6 

6 

0  1 

1 

1 

i 

1 

1  1 

1 

1 

7  7 

7 

7 

^  7 

/ 

7 

-  c/2 

0 

0 

0 

3 

3 

3 

3 

3 

3 

6 

6 

6 

-  c 

1 

1 

1 

1 

7 

7 

7 

7 

-3c/2 

0 

3 

3 

6 

—c/2 

No.  At<xns  1 

1 

3 

3 

1  7 

9 

15 

7  13 

27 

31 

7  19 

33 

45 

atesns,  while  the  third  contains  45  atesns  and  seven  layers.  All  three  cylinders  have 
a  radius  R,  which  includes  a  one  unit  cell  step  along  the  a-direction.  Be-Be 
internuclear  distances  are  taken  as  the  same  as  the  bulk  hep  metal  (a-2.29A, 
c»3.58A)(81.  The  resulting  cylinders  have  point -group  syonetry,  and  the 


adsorption  of  atoauc  hydrogen,  one  on  the  top  surface  and  one  on  the  bottoa  surface 
of  all  three  cylinders,  is  modeled  such  that  the  three-fold  symetry  is  preserved. 
The  same  aj^lies  to  the  adsorption  of  atonic  cesium  and  oxygen.  The  Be-adsorbate 
internuclear  dlatances  are  optimized  in  self-consistent  field  (SCF)  calculations. 
The  19-  and  33-atom  cesiated  and  hydrogenated  cylinders  are  shown  in  Figures  2  and  3. 


The  45-ttom  cylinder  is  not  shown,  but  it  is  similar  to  the  one  pictured  in 
Figure  3,  except  that  it  contains  seven  layers  instead  of  five,  with  the  two  surface 
layers  identical  to  the  surface  layers  of  the  19  atom  cluster  (Figure  2) .  Finally, 
the  simultaneous  interaction  of  oxygen  atoms  and  hydrogen  atoms  and  a  cesium 
overlayer  are  modelled  by  the  addition  of  either  H  or  0  to  the  top  and  bottom  layers 
of  Be.r  between  "capping"  Cs  atoms.  That  is,  CsHBe-c,HCs  or  CsOBe.^OCs  (see  Figure 
2) .  ^Results  obtained  using  these  systems  are  compared  to  those  calculated  by 
treating  identical  bare  and  cesiated  Be  cylinders. 

Calculations  were  accomplished  on  a  Cray  X-MP  supercomputer  using  point-group 
symmetry  optimized  programs  based  on  the  "equal  contribution  theorem"  for 
two-electron  integrals ( 9] .  Ab  initio  restricted  closed-shell  and  restricted 
open-shell  Hartree-Fock  calculations,  each  corresponding  to  an  average  energy  of 
configuration,  were  carried  out  on  numerous  low-lying  states  of  each  cluster.  (The 
average  energy  of  configuration  is  defined  as  the  weighted  mean  of  the  energies  of 
all  the  multiplets  for  the  designated  configuration.)  The  largest  clusters  required 
about  two  to  three  hours  of  Cray  X-MP  time  for  each  geometric  orientation.  The 
following  basis  sets  of  contracted  Gaussian-type  functions  were  used  for  beryllium, 
cesium,  hydrogen,  and  oxygen,  respectively:  (3s2p) /<2slp>(10] ,  (SsSp) /<3s2p>lll) , 
(4slp) /<2slp>(12) ,  and  (484p) /<2s2p>(13) .  Ab  initio  effective  potentials  (EP) [14] 
were  used  to  represent  the  Is  core  electrons  in  Be{10]  and  0(13],  and  the  ls-4s, 


2p-4p,  and  3d-4d  core  electrons  in  Cs(lll.  Relativistic  effects,  especially 
important  for  Cs,  were  incorporated  into  the  EP's(14].  Cluster  binding  energies 
calculated  relative  to  the  completely  dissociated  clusters  are  given  in  Table  2. 
Total  adsorption  energies  given  for  the  cesiated  and  hydrogenated  clusters  were 
calculated  relative  to  the  SCF  energies  of  the  naked  clusters  plus  that  of  two 
adsorbate  atoms.  For  0  and  H  plus  Be^^Ca2,  energies  are  relative  to  this  cluster. 


Table  2.  Restricted  Hartree*Fock  Energies  of  Be^^  and  BtyX^j^ttSters 

Occupied  Mb's*  bTF  Adsorf»idftllMw 

Cluster  a^'  a2'  a2*  e'  e"  a^”  (kcal/mole)  (kcal/asA^MH^ 


4  I  342  0  325.25 

4  0  3-4  2  0  394.70 


Be.gCs,  6 

1 

5 

5 

■5 

w 

0 

402.99 

38.87 

[Be,  (jCs.,) 6 

1 

5 

5 

3 

0 

405.37 

4 

1 

3 

4 

2 

1 

393.76 

34.25 

4 

1 

3 

4 

2 

0 

579.66 

Be,,  6 

1 

5 

6 

4 

1 

724.94 

[86,3)'  6 

1 

5 

6 

4 

1 

815.76 

Be,,CsJ’'* 

1 

7 

7 

6 

1 

768.46 

21.76 

(Ee,,c4l  ^,5^ 

t 

4 

-T 

7 

6 

1 

779.74 

"it  2 

1 

c 

6 

4 

1 

727.32 

1.19 

1 

5 

6 

4 

1 

818.70 

8 

7 

8 

6 

2 

1078.74 

8 

2 

7 

n 

6 

1196.77 

10 

C 

9 

9 

7 

2 

1108.95 

15.11 

[Be^^Cs,!*^  10 

2 

9 

9 

7 

2 

1151.78 

8 

2 

/ 

8 

6 

2 

1115.07 

18.16 

[Be^^H,!^  8 

2 

7 

8 

6 

2 

1228.41 

Be_0,Cs,^''^  12 

HJ  ^  ^ 

2 

11 

9 

8 

1214.17 

52.61*^ 

[Be^302Cs,)"^  12 

12 

9 

8 

i. 

1263.79 

[Be^^H.Cs^]^  11 

2 

10 

8 

7 

2 

1058.49 

-50.46 

[Be,,H,Cs,)^^  -H- 
45 

10 

8 

7 

2 

1090.06 

Table  2.  (continued) 

^Underline  denotes  orbital  from  which  the  electron  was  removed. 
Estate  corresponds  to  the  weighted  average  of  configuration. 


^Open-shell  state: 
Open-shell  state: 
xOpen-shell  state: 
Open-shell  state: 
30pen-shell  state: 
•Open-shell  state: 
•Open-shell  state: 
Popen-shell  state: 
^Adsorption  energy 


(a.')Ma,')‘ 

(e‘)2 

(a/)Ma,')> 

(a  ')Me>)  * 

(e»)‘ 

(a,')  ME')' 

(a,')Ma,-)> 

f  0  on  Cs-coated  Be^^ 


Calculated  work  functions  are  reported  in  Table  3.  The  first  values  in  each  row 
were  obtained  as  the  difference  between  the  total  valence  SCF  energy  of  a  neutral 
cluster  and  that  of  the  ion  generated  by  removing  one  electron  from  the  highest 
occupied  orbital  of  that  cluster.  The  values  isnediately  following  are  due  to 
Koopmans'  theorem  and  correspond  to  the  negative  of  the  energy  of  the  molecular 
orbital  (MO)  from  which  the  electron  was  ionized.  Experimental  values  are  also 
shown. 


Table  3.  Work  Functions  of  Be^  and  B^X,  Clusters 

m  m  ^ 


Cluster 

e(ev)  ® 

5.35,5.04 

Be^gCs.^ 

3.41,2.99 

BeigH. 

4.30,3.96 

"'33 

4.41,4.11 

2.92,2.60 

®«33«2 

4.37,4.09 

3.16,2.94 

2.40,2.18 

®®45«2 

3.40,3.15 

1.97,1.89 

®®45H2C®2 

2.96,2.67 

Be  (expt.) 

3.92 

BeCs  (expt.)^ 

1.6 

BeH  (expt.) 

3.92 

Table  3.  (continued) 


( Bej^  j“B€  2^  jCS2  ] 

1.94,2.05 

A»(Be25-Bel9H2l‘^ 

1.05,1.08 

A^[Be^^>Be23CS2 1 ^ 

1.49,1.51 

A#[Be22~Be22B2 1 

0.05,0.01 

A#(Be^5-Be^5Cs2l'^ 

0.76,0.76 

Ae(Be^5-Be^5H2l‘* 

-0.24,-0.21 

Ae(Be^^CS2~Be^^02Cs2 1 

0.43,0.29 

Am(Be^5CS2-Be45H2Cs2l^ 

-0.56,-0.49 

A*fBe-BeCsl%^p^ 

2.3 

A^[Be-B«H]-^^t. 

0 

*Where  two  values  appear,  these  are  Koopmans' 

.theorem  and  AE(SCF)  values,  respectively. 

°At  optimum  Cs  coverage. 

^Values  correspond  to  a  work  function 
flowering. 

•"s*  ,  values  correspond  to  a  work  function  increase. 

Itoaic  net  electron  charge  values  are  given  in  Tables  •),  5,  6  and  7,  as 
calcuiated  using  a  Mulliken  population  analysis[i5] .  Values  of  R  appearing  in  these 
tables  correspond  to  optimum  atom-to-surface  distances.  These  distances  were 
calculated  from  the  minimum  of  a  curve  fit  to  SCF  eneroies  for  three  or  more 
distances. 


Table  4.  Electron  Populations  of  Be^^,  Be^^Cs.,  and  Be^gH2 


Cluster  Z-coord^ 

Atom 

Label 

No.  ^ 
Atoms^ 

Cluster 

Net  Charge 
per  Atom 

(Cluster) 
Net  Charge 
per  Atom 

Total 

Net  Charge^Difference 
([Cluster]  -  Cluster) 

Be  0 

BeO 

1 

1.43 

1.41 

-0.02 

BeA 

6 

-9.08 

-0.0’ 

0.06 

c/2 

BeC 

6 

-0.09 

0.00 

0.54 

0 

BeB 

6 

-0.07 

9.00 

0.42 

Be,.Cs.  0 

BeO 

1 

1.36 

1.35 

-0.01 

c/2 

BeA 

6 

-0.05 

-0.0’ 

-0.12 

c/2 

BeC 

6 

-0.17 

-0.11 

0.37 

0^ 

BeB 

6 

-0.14 

-0.07 

0.43 

R" 

Cs 

L 

0.39 

0.55 

0.33 

Table  4.  (continued) 


0 

BeO 

1 

1.38 

1.33 

-0.05 

c/2 

BeA 

6 

-0.16 

-0.15 

0.06 

c/2 

BeC 

6 

-0.06 

0.07 

0.78 

0. 

BeB 

6 

-0.02 

0.01 

0.18 

H 

2 

0.04 

0.05 

0.02 

:trcn  Fcouiations  cf  Be^,,  Be,,Cs^  and 


Tluster  Z-cccrd^ 

.Atom 

Label 

No. 

Atoms 

Cluster 

Net  Charge 
per  Atom 

[Cluster) 
Net  Charge 
per  Atom 

Total 

Net  Charge  Difference 
((Cluster)  -  Cluster) 

BeO 

< 

0.96 

0.97 

0.01 

'/2 

BeA 

0.55 

0.53 

-0.14 

Dor 

- 

-n _ i_n 

-0.07 

0.14 

»  *1 

BeB 

—  ‘*1  ' 

-0.29 

0.21 

■ 

-■■'.18 

-0.10 

0.89 

■* 

BeD 

0.20 

0.14 

-0.11 

-  _i  ^  ^  J  ^ 

BeO 

■ 

0 . 94 

0.93 

-0.01 

2 

BeA 

0.5~ 

0.52 

-0 .25 

_ 

BeC 

6 

-0.12 

-0.07 

0.27 

BeB 

-0.34 

-0.33 

0.09 

- 

BeH 

12 

-0.28 

-0.22 

0.67 

BeD 

: 

0.44 

0.40 

-0.08 

?  " 

Cs 

- 

.46 

0.60 

0 . 29 

0 

BeO 

1 

1.00 

1.00 

0.00 

'  ■  :  2 

BeA 

0.45 

0.44 

-0.12 

:/2 

BeC 

6 

-0.06 

-':'.04 

0.12 

0 

BeB 

6 

-0.30 

-0.26 

0.24 

/-* 

BeH 

12 

-0 . 23 

-0,16 

0.84 

r 

BeD 

2 

0.87 

0.82 

-0.10 

H 

2 

-0.24 

-0.23 

-0.02 

^r=3 . 58A 

^No.  of  symmetry  equivalent  atoms 
'R=3.7''A  for  Be  to  Cs  distance. 

"P=1.58A  for  Be  to  H  distance  at  local  minimum  (see  text). 


Table  6.  Electron  Populations  of  Be^^  and  Be^jH2 


Atom  Ho.  j 
Cluster  Z-coord*  Labol  ktOKor 


Cluster  (Cluster!  Total 

Net  Charge  Net  Charge  Net  Charge.Difference 
per  Atom  per  Atom  ({Clustorl^r^g^^luater) 


'  —  .rSA 

r 

"'No.  of  symmetry  equivalent  atoms 

"F  =  3.'’0A  for  Be  plane  to  Cs  distance.  (Be  to  Cs  distance  is  3.93  E.) 
*^P=0.86A  for  Be  plane  to  H  distance  at  local  minimum  (see  te::t)  .  Be  to  H 
distance  is  1.  58A. 


Table  7.  Electron  Populations  of  Be^^02Cs2  and  Be^^H2Cs2 


Atom  Mo . 

Cluster  2-coord*  Label  Atoms 


Cluster  {Cluster)  Total 

.  Net  Charge  Net  Charge  Net  Charge .Difference 


per  Atom  per  Atom  ([C^uatejr]  '  Cluster) 


Be,cO-Cs-  0 
c/2 

BeO 

1 

1.01 

1.02 

BeA 

6 

0.48 

0.48 

c/2 

BeC 

6 

-0.19 

-0.18 

0 

BeB 

6 

-0.36 

-0.35 

- 

BeH 

12 

-0.22 

-0.21 

- 

BeO 

1.15 

1.14 

3c.' 2 

BeF 

0.18 

0.16 

•c/i 

BeG 

6 

-0.02 

-0.03 

0 

-0.98 

-1.00 

Cs 

0.34 

0.86 

Bo  u  r=  0 

'  ^  “  C.'2 

BeO 

1 

0.96 

0.93 

BeA 

6 

0.51 

0.49 

:/2 

BeC 

6 

-0.15 

-0.12 

0 

BeB 

6 

-0.33 

-0.31 

BeH 

12 

-0.26 

-0.24 

- 

BeD 

2 

1.10 

1.09 

;  c 2 

BeF 

6 

-0.01 

-0.04 

BeG 

6 

0.02 

0.05 

H 

2 

-0.16 

-0.15 

Cs 

2 

f^.04 

0.33 

o.#r 

0.12 

-0.02 

-0.12 

-0.06 

-0.04 

1.04 

-0.03 

-0.12 

0.18 

0.12 

0.24 

-0.02 

-0.18 

0.18 

0.02 

0.58 


1  ^  *■■■ 

'---■.z'k  % 

'i'..  ;f  s-.Trmetry  equivalent  atoms 

r=l.50A  fcr  Be  plane  to  0  distance.  (Be  to  0  distance  is  2.0A.) 


=  '  .iX 


■^P-0  ,  ?6A 
=2 . 2k 


for  0  to  Cs  distance. 

for  Be  plane  to  H  distance.  (Be  to  H  distance  is  1.58A. 
for  H  to  Cs  distance. 


2.  RESULTS  AND  DISCUSSION 


RHF  calculations  predict  a  net  energy  gain  of  18.2  and  1.2  kcal/mole • atom  upon 
hydrogen  approach  to  the  Be.,  and  Be,,  surfaces,  respectively.  The  approach  of 
hydrogen  on  the  Be.-  surface,  on  the  other  hand,  results  in  an  energy  gain  of  34.3 
kcal/mole • atom.  Cesium  adsorption  on  the  Be.-  and  Be,,  surfaces  leads  to  energy 
gams  of  39.9  kcal/mole • atom  and  21.8  kcal/mole • atom,  respectively.  Therefore,  while 
cesium  IS  predicted  to  adsorb  on  both  Be.-  and  Be, 3  (i.e.  at  a  hole  and  at  a  head-on 
site),  at  the  SCF  level  of  theory,  hydrogen  adsorbs  strongly  on  Be,  -  and  appear.^  to 
adsorb  moderately  on  (i.e.,  at  a  hole  site)  but  only  negligibly  on  Bejj  (i.e., 
at  a  head-on  site) . 


Table  3  contains  calculated  values  for  the  ionization  potentials  of  all  the 
clusters  studied,  as  well  as  values  for  the  shift  in  ionization  potential  (A^) 
resulting  from  Cs  or  H  adsorption.  Experimentally  determined  values  of  the  work 
functions  due  to  Cs  or  H,  adsorption  on  Be  metal  are  also  included(3] . 

Results  for  the  45  atom  system  indicate  a  slight  rise  by  0.2  eV  in  the  ionization 
potential  of  the  cluster  due  to  H  adsorption,  while  results  for  the  Be,,  system 
indicate  a  lowering  by  1.5  eV  in  the  ionization  potential  of  the due  to  Cs 
adsorption,  and  no  lowering  due  to  H  adsorption.  The  experiMiita^*^ii^ngs  are  2.3 
eV  and  0.0  eV  lowering  in  the  work  function  of  Be  metal  due  to  .9i|M9(p||tion  of  Cs  and 
of  H. ,  respectively.  Adsorption  of  Cs  and  of  H  on  Be..,  h^l^pr,  results  in 
ionization  potential  lowerings  of  approximatley  2.0  eV  ana^l.O  iv,  respectively. 
.Although  the  first  value  is  in  close  agreement  with  experiment,  the  second  is  not. 

These  discrepancies  are  explained  by  referring  to  the  electron  populations  listed 
'.n  Tables  4  and  5.  Total  net  charge  differences,  in  the  case  of  Be,„,  indicate 
participation  in  rharge  redistribution  by  both  the  surface  (BeC,  z-coora.=c/2)  and 
'Tiiddle  (BeB,  z-cccrd.=6)  layers  upon  ionization  of  the  cluster.  In  Be.,Q,  ionization 
:f  the  cluster  results  in  a  contribution  of  0.4  electron  from  Be  atomg^(BeB)  in  the 
middle  layer  compared  to  0.5  electron  for  Be  atoms  on  the  surface.  Thus,  the  middle 
and  surface  layers  cf  the  cluster  are  involved  almost  equally  in  the  electron 
ionization  process.  For  Bej^.Cs.,,  both  the  surface  (BeC)  and  middle  (BeB)  layers 
."ntribute  '^.4  electron,  while’Cs^donates  0.3  electron.  Since  electron  emission  is 
from  the  surface  of  the  bulk  metal,  it  is  concluded  that  Be^^  is  too  small  a  cluster 
to  model  the  Be  metal  surface  involved  in  Cs  adsorption. 

?n  the  other  hand,  the  middle  layer  of  Be,(sH.  contributes  0.2  electron,  with  the 
surface  layers  (Bed  being  predominantly  involved  in  ionization.  H  donates  only  a 
minimal  amount  (0.02  electron)  of  charge.  Part  of  the  charge  contributed  by  the 
middle  (BeB)  layer  ic  shifted  toward  the  BeC  group  cf  atoms  in  Be.gH..  Since  the 
!r=>atest  *harge  perturbation  occurs  on  the  surface,  this  model  at  first  Appears  to  be 
adeauate  for  the  description  of  H  adsorption.  However,  further  analysis  reveals  that 
"hic  IS  not  the  oase.  .^s  was  stated  previously,  the  middle  layer  'BeB)  of  Be,  ^ 

partioipates  appreciably  (donating  ".4  electron)  in  the  ionization  proces^I 
indicatina  that  this  cluster  does  not  reasonably  model  the  bulk  surface.  Adsorption 
:f  '5  or  H  onto  this  surface,  therefore,  yields  information  about  the  interaction  of 
these  -wo  species  with  a  finite  19  atom  Be  cluster,  but  this  information  cannot  be 
interpreted  as  also  applying  to  the  process  whereby  a  Cs  or  H  atom  adsorbs  onto  the 
buDc  3e  surface.  Such  an  e.xtension  may  be  made  only  in  cases  where  the  bare  cluster 
IS  a  good  model  of  the  bulk  metal. 

,  3  five-layer  cluster,  and  Be,,,  a  seven-layer  cluster,  appear  to  model  the 

Be  metal  surface  more  appropriately.  loHization  of  these  clusters  indicates  that  the 
surface  layers  (BeH  for  Be,,  and  BeG  for  Be,,;)  donate  the  emitted  electron.  Although 
the  middle  layers  also  contribute,  the  involvement  of  the  surface  layers  is  three  to 
four  times  as  great.  The  adsorption  process,  then,  is  represented  more  reasonably  by 
usina  Be.,  or  Be,,  clusters  rather  than  Be, ^  to  model  the  metal  surface. 

■)  j  1  r 

Table  6  shows  that  the  surface  (BeG,  z-coord.=3c/2)  layers  of  Be,^  and  Be^^H^ 
are  most  affected  during  the  ionization  process.  Although  the  inner  layefs  of 
are  also  involved  in  ionization,  their  contribution  is  minimal  compared  to  the  charge 
donated  by  the  surface  (BeG)  layers  (0.2  electron  vs  0.6).  As  with  the  Be.g  systems, 
some  of  the  contributed  charge  is  shifted  toward  several  groups  of  atoms  ift  the  Be., 
clusters;  for  example,  BeA  (z-coord.=c/2) .  However,  the  largest  percentage  of  tne 
charge  involved  in  ionization  is  donated  by  the  BeG  atoms  in  the  surface  layer.  It 
IS  also  noted  that  the  net  charge  differences  for  Be^^H^  are  nearly  identical  to 


those  for  Be.,/  the  only  difference  being  0.06  for  the  BeA  layer.  The  approach  of  H 
to  the  Be.(.  ^surface  does  not  affect  the  ionization  process.  The  same  conlusion  is 
reached  frora’an  analysis  of  the  electron  populations  of  Be, 3  and  Be,,H.,  which  appear 
in  Table  5 . 

The  values  of  appearing  in  Table  3  ar'  another  indication  that  H  has  no  effect 
on  the  Be,,  and  Be.,  surfaces.  The  greatest  shift  in  ionization  potential  occurs  for 
the  latter  surface^and  is  0.2.  AO  for  Be,  3  is  practically  ^ 

predicted  to  chemisorb  on  the  Be  surface  and,  in  fact^  in  the  case 

of  the  head-on  approach,  as  modelled  by  Be,,  (the  *"-«ii^HH|^Kenergy  is  1 
kcal /mol  -  atom)  .  Adsorption  of  H  on  Be..,  as  shown  in  Taa>le.^^p|3|||»ored  by  18.2 
kcal/mol •  atom.  The  optimum  H-to-Be-plafte  distance  is  0^. corresponding 
■•alues  for  adsorption  onto  a  Be,,  cluster  reported  by  Bagna  et  al.  are  42.3  kcal/mol 
for  the  adsorption  energy  anc°  0.95A  for  the  distance  from  H  to  the  surface.  The 
adsorption  energy  of  H  on  Be,.,  although  greater  (18  kcal/mol • atom) ,  is  attributed  to 
a  Itoal  minimum  due  to  state  orossina  and  is  akin  to  physisorption  since  both  the 
copulation  analysis  and  ao  values  indicate  no  effect  on  the  Be,^  surface  due  to 
adsorption  of  atomic  hydrogen. 

Hydrogen  approaches  the  Be. ,  surface  above  a  Be  atom  (directly  overhead  site)  in 
-he  present  calculations.  There  are  three  other  possible  sites  for  adsorption. 
These  are  the  Be-Be  midpoint,  eclipsed  and  open  sites.  The  adsorbate  approaches  the 
?e  surface  directly  between  two  Be  atoms  located  on  the  surface  in  the  first  case, 
’he  second  involves  adsorption  over  a  triangle  of  Be  atoms  located  on  the  surface, 
the  renter  of  which  is  directly  above  a  Be  atom  situated  in  a  layer  next  to  the 
surface  (which  is  the  case  for  Be^^) .  The  open  site  has  the  adsorbate  approaching 

-  h<=  surface  directly  above  the'center  of  a  triangle  of  Be  atoms  located  on  the 

surrace.  However,  in  this  case,  there  is  no  Be  atom  directly  below.  According  to  a 
study  by  Bagus  et  al.,  based  on  SCF  calculations  using  clusters  containing  two  to 

-  hree  layers,  the  directly  overhead  site  is  the  least  stable  for  the  adsorption  of  H 
nt:  ,  j  ■'  Be  surface!  1?).  Be.  is  a  five-laver  cluster  and  thus  should  yield  a 

r-ascnaijie  value  for  the  adsorption  energy  associated  with  the  directlv  overhead 
joprcach  of  an  atom  -r  -.he  surface,  while  adsorption  on  an  eclipsed  site  is 
vdeauately  mo#sled  by  the  seven-laver  Be.c  cvlinder. 

It  is  clear  that  the  surface  layers  (BeH)  of  Be,,Cs-  are  predominantly  involved 
in  -lectron  removal.  These  layers  contribute  0.7  electron,  whereas  the  greatest 
charae  -ontributed  by  the  inner  layers  is  0.3  (see  Table  5).  The  tctal  net  charge 
differences  of  Be,,Cs.  differ  from  those  of  Be,,,  indicating  that  Cs,  unlike  H,  is 
affecting  electron  removal.  In  fact,  Cs  donates  0.3  electron  to  this  process.  The 
effects  :f  charge  redistribution  resulting  from  ionization  are  greater  for  the  inner 
layers  (relative  to  the  surface  layers!  of  Be,,Cs.  than  for  those  of  Be,,H..  The  BeC 
group  cf  atoms  contribute  0.3  electron,  while  the  BeA  atoms  accept  0.3  electron  in 
Be.Cs  .  The  sum  of  these  amounts  still  corresponds  to  less  than  half  cf  the  charge 
contributed  by  the  surface  layers.  In  Ee.,H  ,  the  areatest  contribution  by  the  inner 
layers  is  approximately  one  fourth  of  the  charge  contributed  by  the  surface.  The 
surface  layers  of  Be^^Cs,  are  also  predominantly  involved  in  the  lonicaticn  process. 
However,  the  atoms  in'the  inner  layers,  particularly  those  on  the  outer  edges  of  the 
cylinder,  are  involved  significantly  (<;ee  Table  6)  .  In  this  case,  as  opposec  to 
Be,y  the  cesium  interacts  more  closely  with  the  Be.,  surface,  thus  making  edge 
effects  more  apparent.  In  the  hydrogenated  case,  these  edge  effects  were  not  pre.sent 
since  H  IS  much  smaller  than  Cesium  (0.5  vs.  1.67  for  the  ionic  radius  of  H  and  Cs, 
respectively) .  The  small  but  non-negligible  contribution  of  the  inner  layers  of 
Be.,Cs-  and  Be^^Cs-  to  the  ionization  process  indicates  that  the  discrepency  between 
the  calculated  ana  measured  workfunction  lowering  (Table  5)  may  be  improved  by 
-seating  a  cluster  corresponding  to  a  model  based  on  cylinders  having  the  same  height 


as  Be, ,  or  Be-c  but  a  greater  radius  (see  Table  1  and  Figure  1) .  In  addition  to 
improvements  in^^A*^,  adsorption  of  Cs  on  an  eclipsed  site  is  predicted  to  be  more 
stable,  providing  that  the  cluster  model  has  a  surface  large  enough  to  minimize  the 
edge  effects.  Unlike  the  head-on  site,  this  position  allows  the  Cs  to  interact  more 
closely  with  the  surface  triangles  of  Be  since  there  is  no  center  atom  present. 

Both  oxygenation  and  cesiation  of  Be-c  results  in  a  greater  work  function 
lowering  than  does  cesiation  alone,  while  hydrogenation  of  Be^c^s.  results  in  a  work 
function  increase  (see  Table  3) .  The  electron  population  anaiysip  for  Be.jOo^®? 
indicates  that  the  emnitted  electron  is  predominantly  vacating  a  Cfltli  orbital.  In 
Be,cH.,Cs.,,  however,  t.he  Be  cylinder  is  more  invovled  in  ion^Mlfrion,  donating 
apptoximltely  40%  of  the  emitted  charge  (see  Table  7). 

Experimentally  Cs  adsorption  on  Be  results  in  work  a  function  lowering  of  2.3 
“V,  whereas  exposure  of  a  Be  surface  to  H.  yielded  no  observable  change(3].  That  is, 
the  Be  substrate  Auoer  sianal  was  not  attenuated  in  the  case  of  H.  and  there  was  no 
ciscernable  work  function  shift. 


3 .  CONCLUSIONS 


.Analysis  of  the  Mulliken  populations  for  Be.g  indicates  that  this  cluster  is  too 
small  to  model  the  bulk  surface.  However,  the  population  analyses  for  Be,,  and  Be.c 
inaioate  that  these  clusters  are  adequate  models  of  the  Be  metal  surface  involved  in 
•byarcaen  adsorption.  A  wider  cylinder  is  needed  for  the  adequate  modelling  of  cesium 
adsorption.  The  adsorption  energy  of  Cs  on  the  Be,j  cluster  is  21.8  kcal/mol-atom. 
The  decrease  in  the  ionization  potential  of  Be,,Cs.  relative  to  Be,,  is  calculated  to 
1.5  eV.  The  e.xperimental  work  function  lowering  resulting  from  Cs  adsorption  is 
measured  as  2.3  eV.(3)  The  adsorption  energies  of  H  on  Be,_,  and  Be,,  3re  1  and  18 

kcal/ mol • atom,  respectively.  The  first  value  is  negligible,  while  the  latter  is 
attributed  to  a  local  minimum  since  the  population  analysis  and  AO  values  for  both 
Be  and  Be,^  indicate  no  effect  on  these  clusters  due  to  H  adsorption.  Experiment 
'^'•-ais  no  adsorption  ^f  H  on  Be  metal,  thus  ruling  out  dissociative  chemisorption. 
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Layered  aaieoirapic  eeeeiooadiictar  ciuctert  of  the  type  Pb>I]<  arc  etadied  in  ab 
Fock  calealatiooe  that  indade  relathrietic  effects.  Total  energies,  intennelsar  dM 
atoode  charges  are  calcalated  Large  increases  in  the  equilibriun  iodiaadaai  inn 
resalt  in  as  much  at  a  10%  redaction  of  the  in-plane  intemudear  separetiaM  anil 
eiactronic  spectra,  which  are  attributed  to  quantum  tire  effects  observed  in  reeai 
studies. 


lattto  Hartree- 


The  properties  oi  dusters  are  attracting  increased  at¬ 
tention  because  they  represent  an  intermediate  physical 
regime  where  neithw  atomic  nor  solid-state  descriptions 
seem  toully  adequate.'*^  Besides  interest  in  fundamen¬ 
tal  questions  such  as  how  electronic  structure  evolves  as  a 
function  of  size,  there  is  also  the  poonbility  of  using 
finite-size  effecu  to  create  new  matadais  having  unique 
electronic  and  optical  properties.*'*  In  modem  semicon¬ 
ductor  research,  for  eaampie,  it  is  commonplace  to  fabri¬ 
cate  struaures  whoae  physical  dimensions  are  consider¬ 
ably  smaller  than  carrier  wavelengths,  which  leads  to 
physical  phenomena  that  cannot  be  ofaeerved  in  bulk 
matter.*  Ousters  containing  from  about  10  to  i(XX) 
atoms  are  partienh^  ooteworthy  hecauae  their  proper¬ 
ties  set  the  HWnM  Bmitt  on  the  use  of  physical  size  to 
alter  the  hehavlor  immatter. 

Many  ^xperlnliimi  techniques  are  now  available  to 
syntbeiiai  ehmiM^Inth  in  gas  and  condensed  phases,^'** 
but  as  ydt  thde  is  no  generally  available  experimenul 
procedure  to  determine  the  structure  and  spectra  of  the 
clusters  produced.  Hence,  theoretical  insight  into  cluster 
structure  as  a  fwMiion  of  cluster  size  takes  on  a  major 
role.  The  computational  difficulty  inherent  in  treating 
large  systems  comprised  of  heavy  elements,  however,  has 
precluded  rigorous  ab  initio  studies  on  such  clusters.  Up 
to  now.  theoretical  treatments  on  systems  of  this  type  have 
largely  neglected  relativistic  effects,  as  well  as  certain  ex¬ 
change  interactions  among  electrons.  Consequently,  the 
majority  of  ab  initio  work  in  semiconductor  systems  has 
appear«l  on  small  clusters  of  light  elements  such  as  sil¬ 
icon. '*'*  In  such  strongly  covalent  systems,  it  is  found 
that  small  isolated  clusters  tend  to  reconstruct  and  to  form 
tightly  bound  symmetrical  structures  with  a  high  degree 
of  coordination.  Moreover,  this  tendency  to  minimize  the 
number  of  exposed  dangling  bonds  has  the  consequence 
that  the  bulk  silicon  structure  is  recovered  only  for  very 
large  clusters. 

In  this  Rapid  Communication,  we  report  the  results  of 
rigorous  ab  initio  quantum-mechanical  calculations  on 
semiconduaor  clusters  containing  as  many  as  21  heavy 
atoms.  All  of  the  theoretical  resulu  presented  here  were 
calculated  using  fully  ab  initio  procedures:  that  is,  no  ex- 
penmental  parameters  were  used  and  all  of  the  terms  in 


the  Hamiltonian  were  reuined.  Since  no  integrals  were 
apfuoximated  or  disregarded,  all  Coulomb  and  exchange 
interactions  among  electrons  were  explicitly  evaluated  for 
each  cluster.  Furthermore,  relativistic  effects,  which  arc 
significant  in  these  heavy^lement  systems.  '*  were  also  in¬ 
cluded  in  the  calculations. 

The  structures  of  the  layered  semiconductors  investigat¬ 
ed  here  exhibit  a  radically  different  dependence  on  cluster 
size  than  do  those  of  isotropic  silicon.  Unlike  clusters  of 
Si.  layered  semiconductors  of  the  type  Pb^Ia  essentially 
retain  the  structure  of  the  bulk  crystalline  lattice  in  the  la¬ 
teral  direction  for  the  range  of  geometries  studied.  We 
find,  however,  that  there  can  be  significant  electronic  cou- 
(rfing  between  interlayer  perturbations  and  intralayer  lat¬ 
tice  constantt  in  these  clusters.  Indeed,  for  very  large  in¬ 
creases  in  the  Pb-I  interlayer  spacings,  we  see  in-plane  lat¬ 
tice  constants  decrease  by  as  much  as  20%.  With  these  re¬ 
sults.  it  might  be  possible  to  explain  recent  compelling 
scanning  tunneling  images  of  structurally  related  Bil; 
clusters.'^ 

Two  cluster  models  are  studied,  one  of  Pb6li2  and  two 
of  Pb7li4  stoichiometry.  These  are  shown  in  Fig.  1.  Each 
has  the  layer  symmetry  of  the  bulk  crystal  and  is  com¬ 
prised  of  a  plane  of  lead  atoms  above  which  are  situated 
two  triangular  layers  of  iodine  atoms.  Pb«Ii2  has  D}h 
point-group  symmetry  (hexagonal  close  packing),  while 
both  Pb7li4  clusters  have  Du  point-group  symmetry  (cu¬ 
bic  close  packing).  All  three  represent  fragments  of  the 
bulk  crystalline  solid  (one  unit  cell  thick)  having  iodine 
layers  slightly  perturbed  from  the  bulk  arrangement  to 
preserve  the  symmetry  of  each  cluster.  Pb-Pb,  I-l,  and 


FIG.  1.  Pbilii  clusters  (jr  *6.7). 
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Pb-I  intemuciear  distances  were  initially  determined  from 
the  lattice  constants  of  crystalline  Pblj.  All  intemuciear 
distances  were  then  optimized  in  ab  initio  Hartree-Fock 
calculations,  consistent  with  the  variation  principle  and 
within  the  symmetry  constraints.  Self-consistent  field 
(SCF)  energies  for  the  lowest-lying  molecular-orbital 
electron  configurations  were  calculated  for  each  geometric 
orientation  of  the  clusters.  SCF  calculations  on  heavy- 
element  systems  based  on  relativistic  effective  core  poten¬ 
tials  (REP’s)  have  been  shown  to  yield  reliable  geometric 
parameters.  In  addition,  this  method  yields  fully  opti¬ 
mized  molecular  orbitals  for  each  electronic  sute  treated, 
in  contrast  to  many  sundard  methods  based  on  band 
theory.  REP’s  (Ref.  19)  were  used  to  represent  the  l5-4r 
2p-4p,  and  3d -Ad  con  electrons  in  I  (Ref.  20(a)]  and  the 
Ir-Ss.  2p-5p,  3d-Sd.  and  4/  core  electrons  in  Pb.^**** 
Relativistic  effects  cannot  be  ignored  for  systems  contain¬ 
ing  heavy  elements  such  as  Pb  and  I.  Basis  sets  of  con¬ 
tracted  Gaussian-type  functions  were  used  to  represent 
the  valence  orbitals  of  Pb  and  I.^  Calculations  were  ac¬ 
complished  on  a  Cray  X-MP  supercomputer  using  pro¬ 
grams  based  on  the  “equal  contribution  theorem”  for 
two-electron  symmetry  orbital  integrals.^'  These  ab  initio 
restricted  closed-shell  and  restricted  open-shell  Hartree- 
Fock  linear  combination  of  atomic-orbital-molecular- 
orbitai  (LCAO-MO)  calculations  consumed  approxi¬ 
mately  20  min  of  Cray  X-MP  time  for  each  nuclear 
configuration.  The  geometry  optimizations  required  a 
minimum  of  ten  orienutions  per  cluster. 

Energy-optimized  interlayer  and  iniralayer  distances  of 
Pbtlii  were  found  to  be  1 1%  expanded  and  1%  contracted, 
respectively,  relative  to  bulk  Pbl].  This  swollen  interlayer 
distance  is  not  unexpected  for  a  Pb«Ii2  cluster  having  Dji, 
point-group  symmetry  because  the  electron-rich  I  atoms 
of  the  top  layer  are  arranged  in  “head-on”  positions  rela¬ 
tive  to  those  in  the  bottom  layer.  This  is  in  contrast  to 
bulk  Pbh,  where  each  iodine  atom  in  a  top  layer  is  posi¬ 
tioned  directly  above  the  midpoint  between  two  iodine 
atoms  located  in  a  bottom  layer.  This  system  tends  to 
contract  laterally  as  the  interlayer  separations  expand. 
For  example,  expansion  of  the  interlayer  separation  by 
40%  results  in  a  7%  contraction  in  metal-layer  lateral  di¬ 
mensions. 

A  total  of  five  MO  electron  configurations  were  studied 
for  neutral  Pbsli:-  Four  of  these  states  were  nearly  degen¬ 
erate  and  lay  1.6  eV  higher  in  energy  than  the  ground 
state.  The  electronic  excitations  involved  orbitals  associ¬ 
ated  with  both  iodine  and  lead  atoms.  Furthermore,  five 
of  the  six  lowest  unoccupied  MO's  (LUMO’s)  for  each 
electron  configuration  have  energies  lower  than  —  2.7  eV, 
indicating  that  the  neutral  system  is  electron  deficient. 
This  conclusion  is  supported  by  means  of  Koopmans' 
theorem.-* 

If  the  cluster  wave  function  'K2N).  a  2A’-clectron  an- 
tisvmmetrized  product  of  MO’s  IpmK  is  stationary  and  Om 
IS  an  eigenfunction  of  the  Fock  operator,  then  Y(2jV  +  1 ). 
wherebv  an  electron  is  added  to  is  also  stationary  with 
respect  to  any  further  vanations  in  the  orbitals.  This 
means  that  the  orbital  energy  corresponds  to  the  elec¬ 
tron  affinity  of  the  neutral  system,  providing  that  Cm  is 
jnoccuDied  Therefore,  systems  possessing  very  low-lying 
LIMO'^  .ire  electron  deficient:  i.e..  the  addition  of  one  or 


more  elearons  to  these  orbitals  will  increase  the  overall 
stability  of  the  system. 

To  satisfy  this  electron  deficiency,  we  studied  a  doubly 
charged  anionic  Pbsiij  cluster.  We  also  believe  that  the 
anion  is  the  most  natural  cluster  to  consider  from  the  ex¬ 
perimental  point  of  view  since  clusters  of  this  type  are 
prepared  in  solution  as  colloidal  particles  in  the  presence 
of  excess  iodide  ions.^^  The  total  energy  of  Pb«Ii2^~  is 
more  than  6  eV  lower  than  the  corresponding  neutral 
species.  In  this  case,  nearly  all  of  the  LUMO’s  have  ener¬ 
gies  greater  than  zero  (^2.7  eV),  indicating  a  quenching 
of  the  system’s  electrophilidty.  ' 

A  geometry  optimization  ww'Hdio  carried  out  for 
Pb6li2^~-  The  equilibrium  inteiiqwr  and  intralayer  dis¬ 
tances  were  found  to  be  expanded  by  14%  and  1%,  respec¬ 
tively,  relative  to  bulk  Pbl2.  This  cluster  is  also  slightly 
swollen  with  respect  to  neutral  Pb6li2,  but  this  is  expected 
due  to  the  increased  ionicity.  As  with  the  neutral  cluster, 
expansion  of  the  interlayer  separation  results  in  a  contrac¬ 
tion  of  intralayer  distances,  although  not  as  pronounced  as 
for  the  neutral  system:  that  is,  a  40%  expansion  in  inter¬ 
layer  separation  results  in  a  3%  contraction  in  lateral  di¬ 
mensions  relative  to  the  bulk  crystal.  As  previously  men¬ 
tioned,  results  of  SCF  calculations  on  heavy-element  sys¬ 
tems  have  been  shown  to  be  reliable,  with  the  accuracy  of 
the  method  increasing  as  the  number  of  electrons  treated 
explicitly  increases. Therefore,  SCF  results  for  Pb*- 
Ii2*“,  a  1 1 0-valence-electron  system,  are  expected  to  be 
reliable,  and  any  structural  differences  between  our  mod¬ 
els  and  those  of  bulk  Pbh  are  attributable  to  cluster 
effects. 

A  total  of  ten  low-lying  electron  configurations  of  the 
cluster  were  studied  at  both  the  calculated  equilibrium  in¬ 
teratomic  distances  and  at  a  geometry  corresponding  to  a 
60%  expansion  in  the  interlayer  distance  and  a  10%  con¬ 
traction  in  lateral  dimensions.  These  percentages  were 
chosen  from  the  results  of  a  series  of  calculations  at 
geometries  ranging  from  30%  contraction  to  100%  expan¬ 
sion  of  the  lattice  constants  in  the  vertical  direction  and 
30%  contraction  to  20%  expansion  in  the  lateral  direction. 
All  of  the  open-shell  configurations  were  treated  as  pure 
spin  triplet  states:  i.e..  the  cluster  wave  functions  are  spin 
eigenfunctions  having  total  spin  quantum  number  S'*!. 
Excitation  energies  relative  to  the  lowest  state  found  are 
given  in  Table  I.  while  Table  II  reports  electron 
charge/atom  values  calculated  based  on  a  Muliiken  popu¬ 
lation  analysis*^  for  the  lowest-lying  state  at  each  of  the 
aforementioned  two  geometries. 

The  energy  separations  for  equilibrium  Pbeli:'”  listed 
in  Table  1  are  clearly  bandlike,  with  the  two  lowest-lying 
slates  nearly  degenerate  in  energy  and  slates  3  through  8 
lying  at  or  near  the  band  gap  of  bulk  crystalline  Pbl:  at 
2.5  eV  The  last  state  corresponding  to  Pbsli:'”  listed 
in  Table  1  is  considerably  blueshifted  relative  to  the  band 
gap  of  the  bulk.  The  same  type  of  energy  separations 
occur  at  an  expanded  interlayer  distance,  the  only  excep¬ 
tion  being  for  state  4,  which  has  an  energy  value  between 
three  essentially  degenerate  low-lying  states  and  four 
states  grouped  near  the  band  gap  of  the  bulk  solid.  A  to¬ 
tal  of  45  electron  transitions  among  the  ten  lowest-lying 
states  found  were  examined.  The  density  of  states  shown 
in  Fic  2  indici.ifc  th.it  these  transitions  are  peaked  at  0  2 
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TABLE  1.  Electronic  tUto  of  Pb.Ib  cluiten  (eV). 


Cluster 

State  No. 

A£ 

(EquiL)* 

A£ 

(Expanded  Z)* 

Pbsl.j'-" 

1 

0.00 

0.00 

2 

0.03 

0.14 

3 

114* 

0.47 

4 

121 

1.19* 

5 

125 

1.31 

6 

136 

105 

7 

141 

112 

8 

149 

139 

9 

153 

147 

10 

4.84 

3.06 

Pb,I,4*-" 

1 

0.00 

0.00 

2 

0.26 

0.35 

(five  layers) 

3 

0.28 

0.76 

4 

0.35 

0.95 

5 

0.59 

1.11 

6 

0.63 

114 

7 

0.75 

171 

8 

0.76 

190 

9 

0.99 

3.29 

10 

1.01 

3.49 

II 

1.26 

3.99 

Pbrlu-’-' 

1 

0.00 

0.00 

2 

0.68 

1.71 

(three  layers) 

3 

1.54 

1.99 

4 

1.80 

139 

5 

1.83 

3.17 

. 

^  6 

119 

4.29 

7 

139 

4.67 

1.  - 

■4  8 

3.27 

4.75 

9 

3.31 

10 

4.97 

‘See  text.  "See  Fig.  Kb). 

"SeeFif.  !(•).  'See  Fig.  I  (c). 

'Lowest-lying  closed-thcli  singlet  state. 


2.5,  and  4.9  eV.  with  the  second  peak  lying  near  the  exper¬ 
imental  bulk  band  gap  of  Pbl:  and  the  third  blueshifted 
relative  to  the  band  gap  by  2.3  eV. 

Charge/atom  values  appearing  in  Table  11  indicate  that 
the  predominant  interlayer  interactions  occur  between  the 
atoms  comprising  the  smaller  Pb  and  1  triangles  lPb.t  and 
1 4  in  Fig.  I  (a)}  for  the  lowest-energy  state  of  Pbelii"'  at 
both  geometries  treated.  This  is  also  the  case  for  the  other 
eight-electron  configurations  of  Pb«,l,i-"  listed  in  Table  I. 
with  the  exception  of  the  highest-energy  state  (state  10). 
Comparison  between  charge/atom  values  for  this  state  at 
the  equilibrium  geometry  and  at  an  expanded  interlayer 
separation  shows  that  the  interlayer  interactions  predom¬ 
inantly  involve  both  tvpes  of  Pb  and  1  atoms  [Fig.  1(a)]. 
No  differences  arc  seen  between  charge/atom  values  for 
the  ground  state  and  t.ie  corresponamg  values  for  each  of 
the  higher-energy  states  (sutes  2-9  in  Table  1)  at  both 
geometries  treated  indicating  that  the  elearonic  transi¬ 
tions  are  predominantly  intralayer  and  occur  among 
symmeirv-relaied  atoms.  A  more  detailed  analysis  of  the 
Mulliken  populations  reveals  that  these  electron  excita¬ 


tions  correspond  to  meui-meul  transitions.  Differences 
in  charge/atom  values  between  the  highest-energy  state  of 
Pb6li2^~  (state  10)  listed  in  Table  1  and  the  ground  state 
are  as  high  as  0.2  (compared  to  ~0.02  for  sutes  2-9).  In 
this  case,  the  transition  involves  both  Pb  and  I  atoms. 

As  noted  above,  two  cluster  orienutions  were  studied 
for  Pbrlu  (Figs.  I  (b)  and  1  (c)l.  Since  the  neutral  system 
was  found  to  be  electron  deficient,  analogous  to  Pb6li2. 
two  electrons  were  added  to  the  cluster.  A  geometry  and 
electronic  suu  optimization  of  the  first  (Fig.  1(b))  result¬ 
ed  in  a  20%  contraction  and  10%  expansion  of  interlayer 
and  intralayer  distances,  respectively,  relative  to  bulk  Pb- 
l2-  The  lowest  1 1  sutes  (the  energy  for  each  of  which  cor¬ 
responds  to  the  weighted  mean  of  the  energies  of  all  the 
pure  spin  multiplets  for  that  particular  electron 
configuration)  are  nearly  monotonically  distributed  for 
the  five-layered  cluster,  as  shown  in  Table  I.  On  the  other 
hand,  bandlike  behavior  is  observed  when  the  interlayer 
disunces  are  expanded  by  60%  and  lateral  dimensions  are 
contracted  by  10%  relative  to  the  bulk  crystal.  Energy 
gaps  of  2.7  eV  and  higher  may  be  attributed  to  quantum 


TABLE  II.  Population  analysis  for  the  lowest-energy  states 
of  Pb.lj,  clusters  (Ref.  22). 


Cluster 

Atom 

Charge/atom' 

(Equil.) 

Charge/atom' 
(Expanded  Z) 

PW„'-  •> 

Pb., 

0.81 

0.57 

Pbr 

069 

0.67 

I  4 

-0  44 

-0.33 

Ir 

-0.64 

-0.62 

Pbrlu'-' 

Pbn 

I.OI 

-0.29 

(five  layers) 

Pb» 

0.80 

0.55 

I4 

-0.41 

-0.15 

Ir 

-0.59 

-0.51 

lo 

-0.93 

-0.53 

Pbrlu-”" 

Pbo 

0.08 

-0.49 

(three  layers) 

Pb. 

0.80 

0.67 

I4 

-0.53 

-0,45 

Ir 

-0.53 

-045 

Id 

-0.26 

—0.05 

'Electron  units 
°See  Fig  l(a> 

'See  Fig.  Kb) 
"See  Fig.  t (c). 
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size  effects,  which  are  responsible  for  blueshifts  in  the 
spectrum.  The  predominant  interlayer  interactions  occur 
among  the  Pb  layer.  and  Id  [Table  II  and  Fig.  Kb)], 
and  a  comparison  of  charge/atom  values  for  the  ground 
state  and  the  corresponding  values  for  each  of  the  excited 
states  at  both  geometries  indicates  that  all  of  the  transi¬ 
tions  involve  both  Pb  and  I  orbitals. 

The  second  Pb7li4  cluster  studied,  which  contains  three 
layers  [Fig.  Kc)],  was  also  found  to  be  electron  deficient. 
Consequently,  a  doubly  charged  anion  was  investigated. 
Initial  intralayer  atomic  distances  for  this  cluster  corre¬ 
sponded  to  those  in  the  bulk  crystal.  although  interlayer 
atomic  distances  differed  from  those  of  the  bulk  due  to  the 
positioning  of  the  hexagonal  layers.  Upon  optimization  it 
was  found  that  the  lowest-energy  geometry  of  Pb7lu^~ 
corresponds  to  a  2.2%  lateral  contraction  and  a  38.7% 
vertical  expansion  relative  to  the  bulk  lattice  constants. 
The  cluster  was  also  studied  at  a  geometry  corresponding 
to  a  60%  vertical  elongation  and  a  20%  lateral  contrac¬ 
tion.  These  percentages  were  chosen  from  a  series  of  cal¬ 
culations  analogous  to  those  on  Pb«Ii2^~  A  total  of  ten 
electron  configurations  were  studied  for  the  anion  at  both 
the  equilibrium  geometry  and  at  an  expanded  interlayer 
separation.  The  energies  of  these  states  appear  in  Table  1 
(dashed  lines  indicate  energies  greater  than  S  eV).  As  for 
the  five-layered  cluster,  the  states  correspond  to  the 
weighted  mean  of  all  the  multiplets  for  a  particular  elec¬ 
tron  configuration.  Clearly,  the  system  exhibits  a  band¬ 
like  behavior  (Table  1).  The  first  excited  state  at  0.68  eV 
corresponds  to  a  single  electron  transition.  The  next  three 
states  (1. 34- 1.83  eV)  are  double  excitations  from  the 
ground  sute.  while  sutes  6  and  7  (2.19-2.39)  correspond 
to  single  excitations.  The  remaining  sutes  (3.29  eV  and 
higher),  which  are  higher  lying  than  the  bulk-energy  band 
gap.  correspond  to  double  excitations  from  the  ground 
state.  It  is  noteworthy  that  elongated  Pb7li4^~  exhibits 
energy  gaps  of  3.17,  4.67.  and  4.75  eV,  which  may  be  at¬ 
tributed  to  quantum  size  effects,  as  suggested  by  Sandroff. 
Hwang,  and  Chung  (3.42,  3.95,  and  4.80  eV).^^ 
Charge/atom  values  for  all  Pb7li4  clusters  studied 
(Table  ID  indicate  that  the  interlayer  interactions  involve 
both  types  of  Pb  and  I  atoms  in  each  layer  [Figs.  1  (b)  and 
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1(c)].  The  center  lead  (Pbo)  and  the  center  iodines  (lo) 
are  especially  affected  by  interlayer  perturbations  since 
the  charge  on  these  shifts  by  as  much  as  1.3  elearons  (for 
Pbo)  in  the  case  of  five-layered  Pbtlu*".  The  largest  per¬ 
turbations  are  expected  for  these  atoms,  however,  since 
they  lie  along  a  symmetry  axis  and  therefore  experience 
the  greatest  number  of  nearest-neighbor  interaaions. 
Pbo,  in  fact,  is  positioned  at  the  center  of  symmetry.  Al¬ 
though  not  tabulated,  there  are  also  significant  differences 
(as  large  as  0.6)  between  chaifs/atom  values  for  the 
lowest-energy  state  and  the  rni  mipwiilH  values  for  each 
of  the  higher-energy  states  the  Phfli4  anionic  cluster  at 
the  two  geometries  listed  in  Table  farther  indicating 
that  the  electron  exdutions  in  these  dusters  involve  both 
the  Pb  and  I  layers. 

Our  findings  regarding  the  reduction  of  in-plane  lattice 
constants  when  interlayer  distances  are  expanded  from 
their  equilibrium  values  aid  in  the  understanding  of  recent 
scanning  tunneling  microscope  (STM)  images  of  layered 
Bil)  clusters.'^  (These  layered  semiconductors  are  very 
similar  chemically  to  those  of  PbK  and  both  systems  have 
been  studied  experimentally. " )  In  that  study,  atomically 
resolved  images  showed  structures  having  honeycombed 
symmetry  after  colloidal  suspensions  of  Bil)  were  evap¬ 
orated  onto  graphite  surfaces.  Since  the  central  bismuth 
plane  in  Bil)  has  honeycombed  symmetry  (in  contrast  to 
the  hexagonal  symmetry  of  the  iodine  planes),  the  STM 
images  are  most  likely  individual  bismuth  atoms.  In¬ 
terestingly,  the  lattice  constant  of  these  putative  Biiplso 
clusters  was  measured  to  be  0.6S  nm.  roughly  20%  smaller 
than  the  bulk  in-plane  lattice  consunt.  Hence,  these 
STM  images  could  be  revealing  clusters  whose  topmost 
layer  of  iodine  has  been  removed  or  displaced  from  its  un¬ 
derlying  bismuth  plane.  Chemical  oxidation  or  cluster-tip 
interactions  could  be  responsible  for  these  large  perturba¬ 
tions. 
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IhibaomamdlahmuaiUnnerocktakdttdomtoiuidytiiaon^alqeaaammetmeami 
I U  li|Wii  FMm  ebmm*.  As  tks  ctasur  gaameuy  it  energy  opamized  in  several  steps,  the  ttraoerc  distends  into 
l^sndipkeMisOwkinnsBadwiGinngslKstfnMliendlByer.  At  tin  eqndibrinin  geometry  for  the  loweeteiectiaiic 
SMumt  are  oomracted  hy  11.3%  oomperad  to  the  expernnentti  bulk  Pblj  vsinct.  Opticsl 
to  be  blne-ebifted  by  0.3-2. 1  eV  relative  to  the  energy  band  gap  of  crytuUine  Pbl,. 


Thare  ia  iperoMint  bnarest  in  esploittng  dw  nniqae  properties 
attribuuMe  to  qiuuuwn  tiae  or  anaatnn  confmement  effecu  in 
the  design  of  optoaleawoic  devioes.''*  In  fact,  it  is  currently 
within  the  bmitt  of  technology  to  tymbaaine  aeaiioandiictor  duster 
systems  that  have  vahunea  and  sul^boe  areas  that  are  of  the  same 
Older  of  magnitnde  as  carrier  wivahwglhi^-*  These  ade-dependem 
phenomena  typically  yidd  bhie-ahified  optical  spectra.^  exdton 
Gonfinement.*^  dfpreased  ionization  energies,*-^  and  nonbulk 
fragment  stroanr^  compared  to  bulk  propertiet. 

In  this  article,  we  address  the  origin  of  quantum  oonfinemem 
phenomena  in  semiconductor  clusters  from  the  perspective  of 
rigorous  ab  initio  molecalar  orbiul  theory. 

In  reoem  yams,  dnaur  sradim  involving  geometry  opdmizatiam 
and  electronic  structure  determinations  have  been  made  possible 
by  the  avsilability  of  superoompuiers.^' '  Yet  the  intermediate 
reghne  between  motaenlat  and  aolids.  where  certain  quantum  size 
effecu  emerge,  it  not  fully  nnderatood.  However,  that  threshold 
may  require  duaunooMaiaing  as  many  as  1000  atoms. where 
rigorous  ab  initio  cakulationt  are  still  not  feasible.  Even  for 
smaller  clutters,  geometry  optimizations  employing  gradient 
methods  require  enormous  oompoutional  effort.*  However,  the 
advent  of  scanning  tunneling  micwacopy  (STM)  has  provided 
useful  information  about  duster  geometry'^'*  although  tip-surface 
interactions  may  perturb  the  surface  and  absorbed  clusters. '* 
Transmission  electron  micraaoopy  (TEM)  is  also  used  to  directly 
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TAILEI:  Eawiy-OptMaad  latm 


(A) 


geometry* 


pair 

0* 

a 

b 

c 

d 

Pb-Pb 

4.56 

4  56 

446 

4.04 

4.04 

U-Ia 

4.56 

4  56 

446 

461 

4.61 

Ia-1o 

4.56 

4.56 

4.46 

4.61 

464 

»^a-1a 

3.22 

300 

346 

3  44 

3.44 

Pho-lo 

3.22 

1.85 

2.57 

2.56 

3.07 

Pb-P 

1.85 

1.85 

2.57 

2.56 

'See  Figure  1 

*Bulk  geometry  (see  ref  21) 

'  Interlayer  dist 

measure  the  size  and  geometry  of  larger  dusters.^' ^  "  Information 
derived  by  using  these  approaches  may  be  used  as  a  surting  point 
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(bi  (c)  (d) 

ngvtl.  Pb,lMGlMtr|HMifta;  (o)i(initMMifWkFbl>(a)ilM 
iiuual  fMMiciry.  (k)  apiiMiMin  ■  nw  4iaNM  if  6«iiMk  (c)  ipu- 
mixauon  in  ihrw  dqrni  if  fmdMB,  (4)  i^tiaiaiin  if  x  ooatiiMwi 
of  lo  nuMiM  in  nddniM  M  (c). 


6  2  0  0  10  1 

7  2  0  2  0  0  0 

8  0  0  0  3  10 

9  2  I  loo  0 

10  0  «  J  r2  0 

*NMBbnr  of  ilnaroii  hi  iti— Q  *  Imdwibie  rep- 

rcaeautioM  of  MO'i  in  foim  fiwotiiiilMry.  'Nambcr  of  eiec- 
trow  in  inner  *ni«Mi  Mtfk. 

Ab  initio  resthctsd  cloned'iheil  end  restricted  open-shell 
Hnrtree-Fock  linenr  oombinntion  of  atomic  orbiuls  molecular 


forthecomputttioaaldiMnRiMiiMafelMMrpaaMry.  Recent 
STM  studies  on  metal  indhw  anaieaiidMMr  dnHn  sanatted 
the  existence  of  fflicr«iMMn  ha«ia|  a  iaywiri  mraetwc,  which 
appear  to  be  cylindrical  ftadaMMa  of  tka  bak.'* 

been  investipted  in  racaM  yaan  hava  haaa  aallad  aapcnnatcc- 
ules><^  and  supmion"  (Ml  naadak).  TlaiBp— pjanileiare 
spherKal  or  cytindiieal  Wk  ftapMatt  iv  wWah  dnamr  prapeniei 
were  derived  thraagliakiMowlBalBliBai.aiybtlwiapiraiiiwii 
are  based  on  jaUianaNMiadllipMaalbaBadaiiaa  Howwar. 
in  either  case.  sarCMe  aid  hManai  faBaMlric  faoaamractinns  due 
to  the  size  of  the  micracladlan  awaa  aa>  treated,  latberegiaw 
of  microciustin  af  ipal  to  ioMnaadiBli  M  (tin  to  sevaral 
hundrad  atooH),  adlte  ciaractarkiiai  MayaBwne. 

.  ^ 


A  raceatlnM  ttkIpdB  Mj  micredaaien  Wicaied  the  size 
distribBtioaiafsiMbaapliyaiiorbid  daitenpMlMdat  12,  IS,  and 
29  A>'*  intheptaaaaiabiaiiiostBdy  acyliadriealttraGnirefor 
s  gas-phase  PbrlHMMT  is  ananid.  A  aomhiBaiion  of  hexagonal 
lead  and  iodine  plaaaa  h  finad  at  balk  lattioe  positions  and  used 
u  the  initial  dastarpaaMryfFifafa  la).  A  poRion  of  cryaiattine 
lead  iodide  appears  in  Figara  lo  (e  ■  4.3d  A,  c  ■  6.96  A).^' 
Interlayer  disianoe  and  in-pinnB  iaieraianBC  disuaoes  woa  initialiy 
set  at  their  bulk  values,  I.IS  and  4.56  A,  respectively  (Table  I, 
column  ol  The  central  atoms  on  the  hexafonal  layers  were  in 
head-on  positions,  leading  to  Pb^-lA  and  Fbo-io  atomic  distanoes 
that  are  shortened  in  comparison  to  the  bolk  (Table  I.  column 
0)  The  diameter  of  this  cylindrical  duster  corresponds  to  the 
lowest  peak  in  the  experimenul  size  distribution.^  The  coordi¬ 
nation  numbers  of  the  cluster  are  2  (Pbo),  4  (PbA).  I  do)-  and 
2  (Ia).  while  in  bulk  Pbl;  they  are  6  and  3  for  lead  and  iodine 
atoms,  respectively.  Ab  initio  studieB  on  noiMoichiometnc  lattice 
matched  clusiera  can  be  found  elsewhere.^-^^ 


orbitals  (LCAO-MO)  calculations  were  carried  out  for  Pbrlu 
dusters  that  explicitly  induded  the  vaknce  elections  of  each  atom 
(6s  and  6p  for  Pb,  Ss  and  5p  for  I).  Basis  seu  of  contracted 
Gauasian-type  functions  were  used  to  represent  valence  orbiuls 
of  Pb^  and  I.^  The  core  deetrons  ( ls-3s,  2p-Sp.  3d-Sd.  and 
4f  for  Fb,**  Ie-4a,  2p-4p.  and  3d-4d  for  l“)  were  represented 
by  ab  initio  efTectne  core  potentials  that  indude  relativistic  effects, 
which  are  not  negligible  for  heavy  dentenu  such  as  lead  and 
iodine.** 

AO  of  the  eneigitt  of  the  highest  occupied  mdecular  orbiuls 
(HOMO’S)  and  lowest  unoccupied  molecular  orbitals  (LUMO's) 
wete  obeemd  to  be  negative  for  the  low-lying  sutes  of  neutral 
Pbrlu  dusten.  To  compensate  for  this  electron  defldency,  two 
extra  elearans  were  added  to  the  valence  MO's  (Pbrl|4**) 
throughout  the  celmlations.  This  procedure  is  consistent  with 
Koopmans*  theeram**  and  was  used  in  earlier  studks  on  Pbrli],’ 
and  in  which  geometries  were  frozen  at  near-bulk 

valian. 

Geometry  optimizations  based  on  toul  valence  energy  calcu¬ 
lations  required  60  independent  geometries.  The  computer  pro¬ 
grams  employed  here  uke  advanuge  of  the  cluswr  symmetry  to 
significantly  decrease  the  number  of  two-electron  integrals  that 
must  be  calculated.**  The  cluster  was  treated  by  using  its  full 
Du  point  group  symmetry.  Toul  valence  energies  and  partial 
and  grou  atomic  populations**  for  1 1  singlet  and  1 1  triplet  sutes 
induding  the  ground  sute  were  calculated  for  each  geometry 
(Table  li).  From  these  sutes.  55  singlet  and  SS  triplet  transition 
energies  were  obuined  as  differences  in  toul  energies  of  initial 
and  final  sutes.  This  procedure  is  more  rigorous  than  conventional 
HOMO-LUMO  gap  approximations  for  transition  energies. 
LUMO's  are  merely  virtual  orbiuls  due  to  the  self-consistent-rield 
(SCF)  procedure  and  are  deficient  in  the  sense  that  ceruin 
two-electron  self-interactions  are  present.***  In  addition,  if  they 
are  used  as  a  final  sute  of  a  transition,  cluster  orbiul  relaxations 
that  accompany  the  transition  are  ignored.*'  These  orbiul  re¬ 
laxations.  as  well  as  the  proper  inclusion  of  all  two-electron  in- 
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Quantum  Conrinemem  Effects  in  (Pbi}>7 

_ geomeiry* _ 

tnni*  tpin* _ a _ b _  c _ d _ 

1-2  T  -0.65  0.05  -0.17  0.23 

S  -0.11  0.06  0.00  0.01 

1-3  T  0.01  0.73  1.73  1.95 

S  0.42  1.00  1.19  1.99 

1-4  T  2.70  2.46  3.07  3.21 

S  2.11  2.07  2.69  2.82 

1- 6  T  4.53  2.68  4.25  4.36 

S  3.98  2.28  3.88  3.97 

2- 5  T  -0.38  1.60  2J1  3.59 

S  -0.97  1.69  2.22  3.88 

2-8  T  -0.48  1.93  3.17  4.45 

S  -1.12  1  95  3.03  4.67 

‘Trinsiiion  i-j:  i  and  >  are  the  initial  and  final  ttaua.  rcipectively 
(see  Table  II).  *Spin  nwltipiicity;  T  «  triplet  and  S  ■  tingiet.  'See 
Figure  I 

teractions,  were  explicitly  included  in  the  preaent  calculations. 

The  transitions  are  chanKterized  by  an  analyiit  of  the  dominant 
electron  population  shifts.  Comparison  of  atomic  orbital  popu- 
lations  for  initial  and  rinal  nates  ^  a  transition  repremnt  changes 
due  to  the  transition  (population  shifts).  Electronic  charge  is 
considered  to  have  exited  from  the  orfaitais  with  the  mon  negative 
population  shifu  and  entered  those  orfaitais  with  the  most  positive 
shifts.  These  characteristia  may  be  denoted  by,  for  example, 
“Pb-Pb  in-plane  transition*,  “I-Pb  interlayer  transition”,  etc. 

Results  and  Diaeuaaiaa 

The  geometry  was  optimized  adiabatically  (viz.,  the  lowest 
clearonic  sute  was  fon^  at  each  geometry)  with  respect  to  the 
toul  valence  energy  at  three  leveb  (Table  ID.  First  two  pa¬ 
rameters,  the  interlayer  and  intralayer  atomic  dmainces  (height 
and  diameter  of  the  ^indrical  cluster),  were  varied  (Figure  lb). 
The  in-plane  imcreti—ir  disunoes  in  the  lead  and  i^ne  layers 
were  assumed  to  be  egual  at  this  level.  This  process  resulted  in 
an  ettimau  of  apprteiimate  cluster  size.  The  in-piane  atomic 
disttnoe  on  oKh  layer  dsrrwssed  by  2.2%  relative  to  the  bulk  lattioe 
consunt  (Tgble  I,  ealMiin  b).  Tlie  Pb^-lA  atomic  disunce  and 
Pb-I  interlayer  datanee  were  larger  than  the  bulk  values  by  7.3% 
and  38.7%.  respectiveiy.  This  vertical  expansion  of  the  cluster 
is  attributed  to  the  choice  of  a  compressed  initial  geometry. 

Single-electron  trauiitiont  from  kiw-lymg  sutes  (sutes  I  and 
2).  chosen  from  among  1 10  passible  transitions,  are  ubulated  in 
column  b  of  Table  III  for  this  case.  Compared  to  the  sutes  for 
the  initially  chosen  geometry,  many  sute  croasinp  were  observed 
(Table  III.  column  a).  Sutes  I  and  2  were  found  nearly  de¬ 
generate.  Transition  energies  from  sutes  I  to  4  and  to  6  (tran¬ 
sitions  1-4  and  1-6).  2.07-2.68  eV.  are  comparable  to  the  ex- 
penmenul  energy  band  gap  in  crystalline  Pbl;  of  137  eV.^’  Since 
the  population  shifu  in  both  transitions  were  small  in  comparison 
to  all  the  others,  it  is  concluded  that  the  wave  functions  of  the 
initial  and  final  sutes  were  very  similar.  As  a  result,  it  is  expected 
that  their  transition  probabilities  are  higher.  The  other  transitions 
(0.05-1  95  eV)  were  found  to  lie  below  the  band  gap.  At  this 
suge.  the  in-ptane  atomic  distanoes  and  optical  properties  are  seen 
to  still  be  bulklike. 

Second,  three  geometric  parameters  were  vaned:  the  interlayer 
disunce.  diameter  of  the  middle  layer  (lead  layer),  and  the  di¬ 
ameter  of  the  top  and  bottom  layers  (iodine  layers)  (Figure  Ic). 
The  in-plane  interatomic  disunces  in  the  lead  layer  contracted 
by  1 1 .3%  while  the  intralayer  iodine  disunces  suyed  essentially 
the  »me  ( 1 .2%  elonpted)  compared  to  the  initial  geometry  (bulk 
lattice).  The  interlayer  disunce  was  found  to  be  nearly  equal  to 
that  of  geometry  b  of  Figure  I  (contracted  by  0.01  A).  This  result 
showing  a  contraction  of  the  lead  layer  is  consistent  with  the 
experimental  STM  studies  on  surface  meui  clusters  of  Bi'*  and 
of  Au  and  Ag  ”  In  all  of  these  cases,  the  meul-meul  interatomic 
distances  were  observed  to  be  contraaed  relative  to  bulk  values 

;:i  Hirbene  G  .  Touaiii.  E  KCA  Hn  IV75  So  4i> 
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FignK  1  Electraeic  sutes  of  Pbrlu’*  clusieis  at  each  level  of  geometry 
optimization.  T  and  S  indicate  triplet  and  singlet  suits,  respectively. 
Nnmbers  correspond  to  Table  II. 


Figure  2  shows  that  this  geometry  deformation  resulted  in  a 
large  spreading  of  energy  sutes.  Toul  energies  were  lowered  by 
1 . 14-1 .36  eV  in  the  low-lying  sutes  I  and  2.  However,  energies 
were  lowered  by  only  0.12-0.68  eV  in  sutes  3,  4,  5,  and  8  and 
even  increased  in  the  remaining  sutes  by  0.28-1 .46  eV  relative 
to  geometry  b  of  Figure  I .  As  suggested  the  spreading  of  suta 
shown  in  Figure  2.  vibronic  properties  of  the  clusters  appear  to 
be  emerging  at  this  level  of  geometry  optimization.  That  is,  each 
sute  has  its  unique  equilibrium  geometry  with  a  deep  and  narrow 
potential  well.  However,  consistent  with  the  Franck-Condon 
principle,  vertical  electronic  transitions  with  no  simuluneous 
geometry  shifts  are  expened.  As  a  result,  the  energies  of  singlet 
electron  transitions  1-4.  1-6,  and  2-8  arc  blue-shifted  from  0.12 
to  1.68  eV  relative  to  the  energy  band  gap  (Table  III.  column  c). 

The  emergence  of  a  blue-shifted  ultraviolet  (UV)  spectrum  for 
geometry  c.  but  not  for  geometry  b,  indicates  that  quantum 
confinement  effecu  in  cluster  systems  cannot  be  discussed  un¬ 
ambiguously  without  first  determining  the  optimized  geometnes. 
Since  both  systems  had  an  identical  ground-sute  electron  con¬ 
figuration.  and  were  nearly  of  equal-size  nuclear  frameworks,  the 
blue-shifted  optical  spearum  must  be  due  mainly  to  geometric 
reconstruction  in  the  lead  layer. 

This  reconstruction  leads  to  difficulties  in  using  a-particle-in- 
a-box  model  with  accompanying  effective  mass  approximation 
in  microcluster  theory.^  In  this  model  quantum  size  effects  arise 
from  the  finite  boundary  conditions.  However,  this  method  may 
not  be  able  to  correctly  estimate  effecu  due  to  internal  recon¬ 
struction,  which  are  observed  to  cause  large  perturbations  in 
transition  energies  in  the  present  calculations.  For  example,  the 
following  equauon  for  the  energy  rise  due  to  oonfinefnent  in  layered 
Pblj  clusters  within  an  effective  mass  approximation  has  been 
used:^ 

£•  (eV)  .  235/Z.,.' +  26.85/Z.,^  (1) 

L„  IS  the  diameter  and  1,  the  height  (both  in  anptroms)  of  the 
cylindneal  cluster.  The  dimensions  of  geometry  b  of  Figure  I 
can  be  estimated  as  *  1 3.38  A  and  L.  *  9.60  A  by  assuming 
;ne  in-piane  interatomic  disunce  as  the  atomic  diameter  and  the 

C'.rcumscTic'ine  cvhndf  as  ;hr 


ThoM  of  feoneiry  c  in  Fignre  I  are  >  13.26  A  (average 
diameter)  and  L  •  9.73  A.  Substitution  of  tbcK  values  aito  eq 
I  gives  energy  increases  for  geometries  b  and  c  of  Fipire  1, 1.60 
and  1.62  eV,  respsetiveiy.  Tte  increment  in  transition  energies 
due  to  variaiian  fiom  geometry  b  to  c  results  in  oniy  0.02  eV.  while 
the  average  tncremem  in  single  clectmn  transitiao  energies  between 
geometries  b  and  c  is  OJOeV  in  the  rigorous  ab  initio  calculations 
(Table  IIIV 

Finally,  geomenyc  of  Fignre  I  was  fmen.  ooept  for  the  central 
iodine  atoms  on  the  tog  and  bottom  layen  (to  atoms).  The  z 
coordinates  of  these  atoms  alone  were  otninuasd  (FigHre  Id.  Table 
I.  column  d).  The  equMihriuni  geometry  in  tta  echeme  corre- 
spoi^  to  the  lowest  total  vnisnee  ener^  all  Pblj  dusters 
studied.  (The  system  was  not  fully  ogtimind  in  all  four  degrees 
of  freedom  beause  of  oomgutatioital  limitations.)  The  Final 
optimum  geometry  is  guasHaghavical.  Variation  from  geomethes 
c  to  d  of  Figure  I  resulted  in  sute  crossings,  blue^ifted  UV 
speara  (Table  Hi.  cdanui  d)  and  the  further  energy  spreading 
shown  in  Figure  2d.  Tmaailiens  1-3. 1-4.  and  1-6  (1.95-4.36 
eV)  are  seen  to  be  predominantly  Pbo-PbA  in*piane  transitions. 
The  transition  energies  incraaisd^  only  0)09-0.22  eV  frornchmer 
c  to  cluster  d  (TaM  III,  columns  c  and  d),  because  these  tran¬ 
sitions  are  largely  Pb-Pb  inplanc  transitions  and  ate  unperturbed 
by  the  changed  positiam  of  the  Ip  atoms.  On  the  other  land,  large 
increases  in  transition  energies  in  transitions  2-5  and  2-8 
( 1 .28-1 .66  eV)  are  attributable  to  the  Iq  locations.  In  these  cases, 
lo-PbA  electron  transfer  is  predominant  The  tiaaaitioo  energies 
2.82-4.67  eV  in  transitions  1-4, 1-6.  2-5.  and  2-8  (Table  III) 
agree  well  with  the  experimentally  observed  optioal  peido  at  3.42 


3.95.  and  4.80  eV'O  or  3.1,  3.6.  and  4  eV.» 

CiaclueisBi 

Geometric  and  optical  propertia  of  layered  Pb^lu  ciusten. 
reported  here  at  several  levds  of  geometry  opthniation.  indicate 
an  evolution  of  quantum  oonfinetnent  effeett.  Surfsoe  and  internal 
geometric  reconstructioru  arc  important  oomponeau  in  under¬ 
standing  the  nature  of  such  effects.  A  carefully  optimized  ge¬ 
ometry  for  layered  Plhlu  indicated  a  drastic  contract ioo  in  the 
Pb-Pb  in-plane  inieratoroic  distances  and  biueihifted  UV  spectra 
compared  to  crystalline  Pbl}.  Bulk  fragmentt  treated  by  usuig 
ab  initio  procedures  without  full  jBoaMVy  optimiatioifa^  or 
a-partide-in-a-boa  modali  (uddi  «|j|Mlui8PBas  appraKimatiom)^ 
are  not  sufficieatly  tigaioui  iy  tjg  t|pi|jiBiguous  treatment  of 
quantum  size  effeett  in  aaruiamiAgllidwpMunrs.  Surhueand 
internal  geometric  rnmnsirufitiOBS  l— ^j^imsie  in  transition 
energies  resulting  in  bhw-ahllhiB  absorption  spectra. 

Cengiete  geometry  ophnuBtieB  is  aagattid  10  lesnh  in  even  larger 
blue  shifu. 
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Ab  initio  study  of  the  effects  of  cesium,  hydrogen 
and  oxygen  adsorption  on  the  work  function  of  beryllium 

M.M.  Mahno  : 

Depanmeni  of  Chemairv  and  Chemuai  Enginetnng.  Sievent  Insiitute  ofTechnotogy.  Hoboken,  NJ  OTOMl 

and 

W.C.  Ennler  ‘ 

Depanmeni  of  Chemtsin.  The  Ohio  Slate  Universiiv.  Columbus.  OH  43210.  USA 
Rtceivcb  1 2  Oaober  1990 


Hydroten.  oiyten  and  oenum  adaorption  on  beryllium  metal  were  modelled  using  a  Be,,  cluster  conuining  seven  layers.  The 
sb  initio  Hsi  irec-Pock  cakulations  employ  efleaive  core  potentials  and  full  D^i  point-group  symmetry.  Cesiation  of  the  surface 
IS  found  10 1  jwrr  the  work  function  by  2.2  eV.  m  good  agreement  with  the  experimental  value  of  2.3  eV.  Adsorption  of  oxygen  on 
ihe  cesMkUl  surface  lowers  the  work  funcuon  by  an  additional  0.3  eV  over  the  2.2  eV  lowering  due  to  cesiation  alone,  whereas 
hydropenatioa  raises  it  by  0.4S  eV. 


It  1$  well  known  that  the  work  function  of  a  meui 
can  be  lowered  by  the  adsorption  of  cesium  ( I  -4  ] . 
Co-adsorption  of  Oj  has  been  found  to  lower  the  work 
function  even  further  (S).  Studies  of  this  type  are 
important  because  of  the  many  technical  applica¬ 
tions  of  low  work  funaion  surfaces,  one  example  of 
which  IS  in  the  generation  of  H~  and  D~  sources  in 
thermonuclear  fusion  research  [6].  For  many  years, 
the  lowering  of  the  work  function  due  to  alkali-meui 
adsorption  has  been  explained  by  way  of  a  mecha¬ 
nism  involving  charge  transfer  from  adsorbate  to 
substrate  (7).  Recently,  however,  a  more  covalent 
interaction  between  the  two  has  been  proposed  ( 8.9  ] . 

In  this  work,  the  effects  of  cesium,  hydrogen  and 
oxygen  adsorption  on  the  work  funaion  of  beryllium 
are  investigated  using  a  45-atom  cluster  to  model  the 
Be  ( 0001 )  surface.  The  use  of  clusters  as  models  for 
surfaces  is  common  [  1 0. 1 1  ] .  Beryllium  was  chosen 
as  the  metal  substrate  for  several  reasons.  First.  Be 
has  a  closed-shell  ground  electronic  state  and  is  the 

Prnnxneni  address:  E)epanmetii  of  Chemisiry  and  Chemical 

Engineering.  Stevens  Institute  of  Technology,  Hobohen.  NJ 

O’OJO.  USA 


smallest  such  atom  to  exhibit  meullic  charaaer.  Sec¬ 
ond.  the  theoretical  results  may  be  compared  to  ex- 
penmenial  work  involving  Cs  adsorption  on  Be  meuI 
[4].  Finally,  the  Be,,  cluster  used  in  this  study  has 
previously  been  determined  to  provide  an  accurate 
model  of  the  bulk  surface  involved  in  adsorption 
( 1 2. 1 3  ].  All  of  the  calculations  are  earned  out  at  the 
Hanree-Fock  level  of  theory. 

The  model  systems  are  taken  as  cylindncal  “plugs" 
from  Be  wafers  having  surfaces  corresponding  to  the 
( 000 1 )  hep  meui  face  [13].  The  cylinder  chosen  for 
this  study  conuins  45  atoms  and  seven  layers  and 
has  a  radius  which  includes  a  one-unit  cell  step  along 
the  a-direction.  Be-Be  intemuclear  distances  are 
equal  to  those  in  the  bulk  hep  meui  (a=2.29  A. 
c=3.58  A)  (14).  This  cylinder  has  Dih  point-group 
symmetry,  and  the  adsorption  of  cesium  and  oxy¬ 
gen/hydrogen.  one  on  the  top  surface  and  one  on  the 
bottom  surface  is  modeled  such  that  the  threefold 
symmetry  and  honzonul  mirror  plane  are  pre¬ 
served.  In  the  case  where  only  one  type  of  adsorbate 
IS  involved,  the  Be-adsorbate  intemuclear  disunces 
have  been  optimized  in  self-consistent-field  (SCF) 
calculations.  The  45-atom  cylinder-plus-adsorbates 
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it  sbowa  in  fig.  I .  Resulu  obtained  using  these  sys- 
tetntaie  compared  to  those  calculated  by  treaung  an 
idaacMilbare  Be  cylinder.  For  cases  where  two  types 
of  admbates  are  involved  (i.e..  Bea^OjCs}  and 
),  no  geometry  optimization  was  per¬ 
formed  due  to  the  compuutional  expense.  Instead. 


Be-adsorbate  and  adsorbate-adsorbate  iniemuclear 
separations  correspond  to  van  der  Waals  radii.  Re- 
sulu  for  these  systems  are  compared  to  those  ob¬ 
tained  by  treaung  identical  bare  Be  cylinders  which 
have  been  cesiated. 

Calculationa  were  canM||||y|p  a  Cray  Y-MP  su¬ 
percomputer  wiaf  on  the  "equal- 

contribution  l||g^ipB**m[^^^pl|Bctron  integrals 
[IS].  Ab  initio and restnaed 
open-shell  HMirrir  each  corre- 

orbital  (MO)  e&^R^nntion,  wefCtarhed  out  on  nu¬ 
merous  low-lying  states  of  each  cluster.  (The  average 
energy  of  configuration  is  defined  as  the  weighted 
mean  of  the  energies  of  all  the  multiplets  for  the  con¬ 
figuration  [16].)  The  calculations  for  the  largest 
clusters  consumed  about  one  hour  of  Cray  Y-MP  time 
for  each  geometric  onenution.  The  following  basts 
sets  of  contracted  Gaussian-type  funcuons  were  used 
for  berylliuim  cesium,  oxygen  and  hydrogen,  respec¬ 
tively:  (3t2p)/[2slp]  [17].  (5s5p)/(3s2p]  [18]. 
(4s4p)/[2s2pl  [19]  and  (4slp)/[2slp]  [20].  Ab 
initio  effective  potentials  (EP)  (21  ]  were  used  to 
represent  the  IscoreelearonsinBe  [17]  andO  [19] 
and  the  ls-4s.  2p-4p.  and  3d-4d  core  electrons  in 
Cs  [  1 8 ].  Relativistic  effecu.  especially  important  in 
Cs.  were  incorporated  into  the  EPs  [21 ). 

Binding  energies  (BE)  calculated  relative  to  the 
completely  dissociated  clusters  are  given  in  table  I . 
(Valence  SCF  energies  for  atomic  Be.  Be*.  Cs.  Cs*. 
O.  Hand  H- are  -0.95083.  -0.65457.  -19.84225. 


Tabic  I 

Rcttncicd  Hanree-Fock  eneriie*  of  Be<i  and  Be4iX,  clunen 


Cluster 

Occupied  MOs*’ 

BE 

(kcal/mol  i 

Adsorption  energ> 

( kcal/mol  atom  i 

a; 

a: 

*7 

c 

C 

a; 

Be.,»' 

8 

2 

7 

7 

6 

1120.00 

(Be.,r 

8 

■1 

7 

7 

6 

1203.66 

Be„C$; 

10 

9 

9 

7 

1131.02 

5.51 

(Be.,Cs,r'" 

10 

■A 

9 

9 

7 

2 

1266.52 

Be„0,C$.  •*’ 

i: 

*> 

1  1 

9 

8 

1378.68 

123.83  " 

(Be.,0,C5,r  •' 

i: 

•> 

10 

9 

8 

1520,99 

(Be.,HjCijl 

10 

9 

9 

7 

1272.94 

70  96  ' ' 

(Be„H,C5:r 

10 

■A 

9 

9 

7 

1397.62 

■ '  Underiine  denotes  orbital  from  which  the  eiearon  was  removed.  *’  Sute  corresponds  to  the  weishted  averafe  of  connsuraiion. 
‘ ’  Open-shell sute:  ( e  )'  Open-shell  state:  (a,  I '(aH' 

'  ’  Adsorption  energy  of  O  on  Cs-coated  Be.,.  ' '  Adsorption  energy  of  H  on  Cs-coaied  Be,,. 
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-19  72876.  -15.64121.  -0.49928  and  -0.44815 
hanree.  respectively,  for  the  basil  sen  used  in  this 
study. )  The  total  adaoiption  energies  given  in  ubie 
I  for  the  cesiated  duster  was  calraiiaied  relative  to 
the  SCF  energies  of  the  naked  cluster  plus  that  of  two 
adsorbate  atoms.  For  BcesOCt)  aad  Be45HiCs2, 
energies  are  relative  to  H  and  O  adsorbed  on  Be45Cs2. 

Theoretical  work  function  are  given  in  table  2.  The 
first  values  in  each  row  were  calculated  as  the  dif¬ 
ference  between  the  total  valence  SCF  energy  of  a 
neutral  cluster  and  that  of  the  ion  generated  by  re¬ 
moving  one  electron  from  the  highest  occupied  or¬ 
bital  of  that  cluster.  The  values  immediately  follow¬ 
ing  are  due  to  iCoopmans'  theorem  and  correspond 
to  the  negative  of  the  energy  of  the  molecular  orbiul 
( MO)  from  which  the  electron  was  ionised. 

Atomic  net  elearon-charge  values  are  given  in  ta¬ 
ble  3  as  calculated  itaiag  a  Mulliken  population  anal- 
vsis  [  22 ).  Value  of  R  appearing  in  these  ubies  cor¬ 
respond  to  optimum  atom-to-sur&ce  disiances  These 
distances  were  calculated  from  a  curve  fitted  to  SCF 
energies  for  at  least  three  surface-adsorbate  sepa¬ 
rations  corresponding  to  the  lowest-lying  electron 
configuration. 

The  toul  adsorption  etier,  for  C9|i«ed  Be45  ap- 
pean  in  ubie  1  na  5.3  kcal  .nol  atom.  As  previously 
mentioned,  ihii  airgy  was  calcutaied  raiaAive  to  the 
toul  SCF  ensr^  of  the  naked  duaier  plus  that  of 
the  two  adsorbaia  atoms  in  their  ground  sutes.  and 
the  value  corrcspuads  to  an  average  energy  of  con- 


Table  : 

Work  functioni of  Pc«t  and  Be^Xj  cluiien 


Cluiier 

<P  (eVi  •' 

Be.. 

4  62.  4  43 

Be.,(.i: 

;.4i,  2.19 

Be..O:Cs: 

181.  1.89 

Be,.H..CS: 

2.95.2.67 

Be  le»pt  ) 

392 

BeCs  ie»pt  )  *' 

16 

.ii>fBe4,-Be.,Ci,)  " 

:.:i.  2  24 

.i0(Be4,C5:-Be4,O.C$j|  " 

0.60.  0.30 

J0(Be.,Cs:-Be.,H,Cs;)  ■" 

-0  54.  -0  48 

.l<J(0e-BeCs|expi 

2  3 

■ '  Whf  rt  iixo  values  appear,  ihese  are  koopmans  theorem  and 
A£(SCF'  values,  respectiveh 
Ref  |4i 

'  Values  correspond  to  a  work  function  lowenng. 

Values  correspond  to  a  work  function  increase 


figuration.  Therefore,  cesium  is  predicted  to  adsorb 
weakly  on  Bcas  ( i  c.  at  a  hole  site  on  the  Be  ( 000 1 ) 
surface )  at  the  SCF  level  of  theory.  Table  2  contains 
calculated  values  for  the  ionization  potentials  of  Beas. 
as  well  as  values  for  the  shift  in  ionization  potential 
(d^)  resulting  from  Cs.  H  and/or  O adsorption.  Ex¬ 
perimentally  determined  veiiioiof’ihe  work  func¬ 
tions  due  to  Cs  adsorption  on  %mpdllo  also  in¬ 
cluded  { 4  ] .  Results  indieoif  a  eV  in 

the  ionization  poiantiai  of  the  <j|aCTrlpa^  Cs  and 
H  adsorption.  The  exponmentm  finding  is  a  2.3  eV 
lowering  in  the  work  function  of  Be  metal  due  to  ad¬ 
sorption  of  Cs. 

It  is  the  surface  layer  of  Bca,  which  is  predomi¬ 
nantly  involved  in  ionization.  This  layer  ( BeG  in  fig. 
I )  contributes  0.3  elearon  to  the  process  ( table  3 ). 
with  Cs  contributmg  0.4  electron.  The  greatest  change 
contnbuted  by  the  inner  layers  is  0.2  electron  ( BeH 
in  fig.  1 ).  However,  comparison  of  the  net  charge 
differences  between  Be45Cs2  and  Be4s  indicates  that 
the  group  of  atoms  most  affected  by  Cs  adsorption 
are  on  the  surface  of  the  cluster  ( BeG.  fig.  1 ).  Cs  ad- 
sorpuon  shifts  the  electron  contributions  of  this  layer 
by  about  0.2  eV. 

Combined  oxygenation  and  cesiation  of  Be4s  re¬ 
sult  in  a  greater  work  function  lowenng  than  does 
cesiation  alone.  Oxygenation  of  the  cesiated  cluster 
further  lowers  the  ionization  potential  by  0.3  eV  for 
a  total  ionization  potential  lowenng  of  2  5  eV  rela¬ 
tive  to  bare  Be45  (table  2).  O  adsorbs  strongly  on 
cesiated  Be45.  the  adsorption  energy  being  123.8  kcal  / 
mol  atom  (table  1).  Net  charge  differences  for 
Be4502Cs7  (table  3)  indicate  that  cesium  contnb- 
utes  significantly  to  the  ionization  process.  In  fact, 
cesium  IS  donating  1 .0  electron  to  the  process,  while 
the  role  of  the  surface  layers  is  subsiantialiy  reduced 
relative  to  cesiated  Be,*.  The  BeG  and  BcF  group  of 
atoms  on  the  surface  of  Be«5Cs,  have  a  charge  of  0.30. 
while  the  charge  on  the  same  layer  in  Be450:Cs7  is 
-0.06. 

Hydrogenation  and  cesiation  of  Be4v  on  the  other 
hand,  result  in  a  nsc  of  about  0.5  eV  in  the  ioniza¬ 
tion  potential  of  Be4«Cs:  (table  2).  Since  cesiation 
of  Bc.v  lowers  the  ionization  potential  of  the  cluster 
by  2.2  eV.  there  is  still  a  total  lowering  of  1.7  eV  in 
the  ionization  potential  of  Be44H7Cs:  relative  to  bare 
Be4v  The  adsorption  energy  of  H  on  cesiated  Be4»  is 
■^0.96  kcal/ mol  atom  (table  1)  Therefore.  H  ad- 
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Table  3 

Atomic  populations  of  Bets.  BcsiCsj.  BetsOjCS}  and  BetiHjCS} 


Cluster 

Z-cootd.  •’ 

Atom 

label 

No. 

atoms  *' 

Cluster 
netcbarfc 
per  atom 

( Cluster  )- 
net  charge 
peruom 

Total 

net  charge  diflerence 
(1  duster)"  -cluster) 

Bc„ 

0 

BcO 

1 

I.OI 

1.00 

^aoi 

c/2 

BeA 

6 

0.42 

adO 

■  4-0.12 

c/2 

BeC 

6 

-0.20 

-0.9  .. 

0 

BeB 

6 

-0.33 

-0.31 

c 

BeH 

12 

-0.19 

-0.17 

c 

BeO 

1.36 

1.34 

ic/2 

BeF 

6 

-0.14 

-0.14 

0.00 

3c/2 

BeG 

6 

0.01 

0.09 

0.48 

Be„Csj 

0 

BeO 

1 

0.97 

0.97 

0.00 

c/2 

BeA 

6 

0.46 

0.45 

-0.06 

c/2 

BeC 

6 

-0.19 

-0.17 

0.12 

0 

BeB 

6 

-0.37 

-0.34 

0.18 

C 

BeH 

12 

-0.23 

-0.21 

0.24 

c 

BeO 

A 

1.22 

1.21 

-oo.-* 

3c/ 2 

BeF 

6 

-0.06 

-0.09 

-0.18 

3c/2 

BeG 

6 

-0.05 

0.00 

0.30 

R  •’ 

Cs 

•) 

0.29 

0.51 

0.44 

0 

BeO 

1 

1.01 

1.02 

0.01 

c/2 

BeA 

6 

0.48 

0.48 

0.00 

c/2 

BeC 

6 

-0.19 

-0.18 

0.06 

0 

BeB 

6 

-0.36 

-0.35 

0.06 

c 

BeH 

12 

-0.22 

-0.21 

0.12 

c 

BeD 

2 

1.15 

1.14 

-0.02 

3c/2 

BeF 

6 

0.18 

0.16 

-0.12 

3c/2 

BeG 

6 

-0.02 

-0.03 

-0.06 

R  •' 

O 

-0  98 

-1.00 

-0.04 

R-’ 

Cs 

•» 

0.34 

0.86 

1,04 

0 

BeO 

1 

0.93 

0  94 

0  01 

t7  2 

BeA 

6 

0.46 

0.45 

-0.07 

c/2 

BeC 

6 

-0.15 

-0.12 

0.13 

0 

BeB 

6 

-0.38 

-  0.36 

0.17 

c 

BeH 

12 

-0.19 

-0  18 

0.20 

r 

BeD 

n 

0,14 

0.13 

-0,04 

ic/2 

BeF 

6 

-0  08 

-0.1 1 

-0.17 

ic/2 

BeC 

6 

-0  04 

O 

o 

\  J 

0,38 

R" 

H 

-0.11 

-0.12 

-0.02 

R  t> 

Cs 

•) 

0.23 

0  44 

0.42 

*'  £■»  3.58  A  *'  No.  of  symmetry  equivaieni  atoms.  /?  =  3.70  A  for  Be  plane  to  Cs  distance,  i  Be  to  Cs  distance  is  3.93  .A. ) 

7?  =  1.50  A  for  Be  plane  to  O  distance.  ( Be  toOdtsunce  is  2.0  A.  1  "  R  =3.1  A  forOioCsdisunce. 

' '  R  =  0.86  A  for  Be  plane  to  H  distance.  ( Be  to  H  distance  is  1.58  A. )  •’  R  =2.2  A  for  H  to  Cs  distance 


sorbs  more  than  twice  as  strongly  on  cesiated  Bc.s 
than  on  bare  Bctv  Net  charge  differences  appearing 
in  table  3  indicate  that  hydrogenation  of  BctsCs;  has 
little  effea  on  the  ionization  process.  Tfie  net  charges 
on  hydrogenated  Bcss  are  nearly  the  same  as  those 


on  cesiated  Be.s,  with  the  greatest  change  occurring 
for  the  BeG  group  of  atoms  on  the  surface  (0.30 
electron  versus  0.38  for  Be^yCsy  versus  Be<,H:Cs:. 
respectively ). 

NeKharge  per  atom  values  shown  in  table  3  for 
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each  Cs  of  neutral  Be45Cs2  (Be4502Cs2)  indicate  that 
the  interaction  between  these  atoms  and  the  naked 
( oxygenated )  surface  is  mote  covalent  in  nature,  with 
the  charge  polarized  inward  toward  the  surftce.  This 
IS  in  agreement  with  previous  theoretical  ( 8  ]  and  re¬ 
cent  experimental  (9]  data  for  alkali  adsorption  on 
meui  surfaces,  lliis  polarization  is  more  pro¬ 
nounced  for  the  interaction  of  Cs  with  the  oxygen¬ 
ated  surface.  The  same  phenomenon  prevails  for  the 
charged  systems.  However,  in  these  cases  it  is  more 
pronounced,  especially  for  the  oxygenated  surface, 
where  the  net  charge  per  atom  value  of  0.9  indicates 
an  ionic  interaction  between  the  Cs  and  the  oxygen¬ 
ated  surface.  As  expected,  the  oxygen  atoms  on  both 
the  neutral  and  charged  species  are  highly  electro¬ 
negative.  The  smallest  dipole  between  the  cesium  at¬ 
oms  and  the  surface  of  both  the  neutral  and  ionized 
system  exists  for  the  hydrogenated  surface,  although 
the  difference  is  qutte  small.  It  is  noteworthy  that  the 
dipole  interamon  between  Cs  and  the  surface  is 
greatest  for  the  system  having  the  lowest  work  func¬ 
tion.  namely  the  604502  surface,  and  smallest  for  the 
system  having  the  highest  work  function,  t.e.  the  hy¬ 
drogenated  6045  surface. 
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Gold  colloidal  suspensions  have  been  deposited  onto  highly  oriented  pyrrolitic  graphite  and  analyzed  by  using  a  scanning 
tunneling  microscope  (STM).  The  size  distribution  of  individual  particles  measured  with  the  STM  ( 16  Sb  1 .5  nml  agreed 
well  with  measurements  made  by  transmission  electron  microscopy.  STM  images  of  colloidal  aggregates  containing  up  to 
I S  individual  particles  could  be  obuined  under  a  wide  variety  of  experimental  conditions.  The  quality  of  these  images  was 
extremely  sensitive  to  the  type  of  aggregating  agent  used,  with  ionic  salts  giving  better  lesulu  than  molmlar  species.  Under 
our  colloidal  deposition  conditions  we  often  observed  highly  ordered  hexagonal  patterns  with  a  periodicity  of  7S-IS0  nm. 
extending  over  areas  as  large  as  I  SO  X  ISO  nm’  Preliminary  evidence  suggests  that  this  hexagonal  lattice  is  associated  with 
the  graphite  substrate  itself,  and  not  with  an  ordered,  twoslimensional  array  of  colloidal  particles. 


Perhaps  the  most  significant  advanuge  of  the  scanning  tun¬ 
neling  microscope  (STM)  over  other  surface  science  tools  is  its 
ability  to  scan  large  areas  of  a  surface  and  to  probe  local,  non¬ 
periodic  structures  with  atomic  scale  resolution.  This  feature  ts 
of  great  utility  in  searching  for  adsorbed  chemical  speaes  ranging 
from  molecules'  to  metallic  and  semiconductor  clusters. 
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Meullic  clusters  adsorbed  on  graphite  have  had  a  special  appeal 
for  STM  workers’^  because  they  are  easily  produced,  are  relatively 
robust,  and  are  readily  differentiated  from  the  graphite  substrate 
through  their  tunneling  charaaehstics. 

Encouraged  by  the  STM  resulu  on  small  meul  clusters,  we 
began  to  explore  the  properties  of  larger  metallic  systems— 
colloidal  particles — deposited  on  the  surface  of  highly  ordered 
pyrolytic  graphite  (HOPG).  In  principle,  many  important  features 
of  colloids  could  be  understood  by  exploiting  the  atomic  resolution 
achievable  with  the  STM.  For  example,  local  variation  of  the 
electrical  double  layer  surrounding  colloidal  particles,  and  the 
geometry  of  adsorb^  molecular  species  on  colloidal  surfaces,  are 
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Figarc  I  STM  imafe  of  itoloted  ooHoidal  goW  particles  on  a  HOPC 
surface.  The  averape  hk  of  the  ooUoidt  is  I S  nm.  The  spiky  noise  in 
the  image  is  bdieved  to  come  from  salts  which  precipitate  on  the  graphite 
surface  during  the  depoMtion  of  the  colloids. 


particularly  relevant  in  the  science  and  technology  of  dispersed 

systems. 

in  this  Letter  we  repon  some  initial  STM  remits  from  colloidal 
gold  deposited  onto  HOPG.  For  a  number  of  reasons,  the  gold 
system  is  an  ideal  suning  point  for  determining  the  utility  of  the 
STM  in  colloid  science.  Not  only  do  simple  chemical  syntheses 
produce  small  particles  with  a  narrow  site  distribution.’  but  both 
the  aggregation  behavior  and  surface  chemistry  of  the  gold 
particles  have  been  explored  in  considerable  detail.'*'*’ 

Our  aqueous  gold  colloids  were  piepared  through  the  chemical 
reduction  of  sodium  tetrachloroaurate  using  trisodium  citrate 
according  to  a  modified  synthesis*  first  reported  by  Hillier  et  al.’ 
This  recipe  results  in  the  formation  of  suble  gold  particles  with 
a  mean  diimeter  of  **'  I S  nm.  These  colioidat  suspensions  remain 
dispersed  for  many  months,  presumably  stabilized  by  citrate  ions 
strongly  lAemitorbed  to  the  colloid  surface. 

Because  the  aqueous  gold  suspension  does  not  wet  the  graphite, 
we  found  it  necessary  to  modify  the  HOPG  surface  in  order  to 
obtain  a  sample  surface  which  was  uniformly  coated  with  isolated, 
unaggregated  colloidal  particles.  This  modification  was  achieved 
by  piacmg  a  small  drop  of  I -propanol  on  the  graphite  surface  pnor 
10  depositing  the  gold  sol.  We  found  that  surfaces  treated  with 
the  alcohol  could  be  partially  dried  under  vacuum  with  no  apparent 
loss  in  wetting  efficiency 

.'Ml  of  our  tunneling  images  were  obtained  with  an  STM  op¬ 
erating  in  air  at  room  temperature."  To  isolate  the  system  from 
low-frequency  vibrations,  the  sample  head  was  placed  on  a  SO-kg 
anvil  suspended  by  rubber  cords  stretched  to  approximately  twice 
their  equilibnum  length.  We  used  Pl/Ir  STM  tips,  and  operated 
the  microscope  at  biases  of  -500  mV  to  -hSOO  mV  and  tunneling 
currents  of  ~  I  nA  in  the  constant  height  mode.  To  ensure  that 
the  STM  was  operating  optimally,  we  would  first  scan  the  graphite 
surface  under  conditions  which  allowed  us  to  resolve  atomic 
features  on  the  graphite  surface.  We  would  then  scan  over  the 
surface  in  search  of  gold  particles. 

Figure  I  displays  the  typical  appearance  of  a  colloid-covered 
HOPG  surface  scanned  by  the  STM  under  these  conditions.  The 
scan  reveals  a  high  density  of  particles  on  a  flat  HOPG  back¬ 
ground  The  stability  and  brightness  of  the  STM  images  indicate 
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the  presence  of  conductive  metallic  panicles  strongly  adsorbed 
on  the  HOPG  substrate. 

To  confirm  that  these  images  were  of  gold  colloids  and  not 
random  adsorbed  impurities  or  graphite  anifaces.  we  measured 
the  dimensions  of  the  panicles  in  the  xy  direaion  (parallel  to  the 
graphite  substrate)  with  the  STM.  This  measurement  gave  an 
average  panicle  diameter  of  1 6  ±  1.5  nm.  Considering  that 
tip-cluster  interactions  can  change  the  apparent  size  of  measured 
features.^  this  diameter  is  excellent  agreement  with  the  14.5-nm 
colloid  diameter  measured  by  TEM  in  this  and  in  previous  work.' 

Though  xy  scans  with  the  STM  aocarately  reflected  gold 
panicle  diameters,  we  fouad  tet  :i>§J|iight  mode  was  not 
quantitatively  reliable.  GenezaHy.  save  panicle  heights 

in  the  2  direaion  varying  bctweea^(|Si  2.0  nm.  a  dimension 
considerable  smaller  than  the  I  ^^teRicie  diameter.  We 
believe  that  the  contracted  dimenaiai  gBM  direction  is  a  result 
of  operating  the  microaoope  in  the  couaiM  bdght  mode  rather 
than  a  flattenmg  of  the  gold  ooUoid  paniGie  on  the  graphite  surface. 
Hence  dimensions  measured  by  tunneling  are  more  reliable  in  the 
xy  than  in  the  z  direction,  and  we  conclude  that  the  STM  can 
be  used  as  a  reliable  and  semiquantitative  tool  for  measuring  the 
shape  and  dimensions  of  colloidal  metallic  particles  deposited  on 
conducting  substrates. 

After  our  initial  characterization  of  single,  isolated  colloidal 
particles  with  the  STM.  we  began  to  explore  the  properties  of 
colloidal  aggregates.  The  particles  were  aggregate  by  adding 
1-3  drops  of  I  M  NaCI  or  I  M  NaCIOj  to  5  cm^  of  the  colloidal 
suspension.  After  the  addition  of  the  salt  solutions,  the  wine-red 
Au  suspension  turned  blue,  indicating  that  the  small  aggregates 
of  gold  were  beginning  to  form.’  At  this  stage  several  drops  of 
ihe  suspension  were  deposited  onto  a  freshly  cleaved  HOPG 
substrate  and  evacuated  to  dryness.  The  water  was  removed  at 
a  reduced  pressure  of  ICT'  Tort  without  allowing  the  sol  to  solidify 
and  sublime.  As  described  earlier,  we  modified  the  graphite 
surface  with  I  -propanol  and  deposited  the  aggregates  both  at 
ambient  temperatures  and  at  1 10  *C.  Elevated  temperatures  did 
not  change  the  struaure  or  distribution  of  aggregates  but  did  speed 
up  the  removal  of  solvent  and  propanol. 

Figure  2  shows  typical  STM  and  TEM  images  of  colloidal  gold 
aggregates  flocculated  under  these  conditions.  Samples  for  TEM 
analysis  were  obtained  by  depositing  the  colloids  on  carbon-coated 
copper  grids.  The  STM  images  of  the  flocculated  colloid  are  very 
similar  to  the  appearance  of  the  aggregates  by  TEM:  the  indi¬ 
vidual  panicles  appear  to  be  physicallv  touching  but  rarely  fused 
together.  As  with  the  isolated  panicles  the  aggregates  were  stable 
for  periods  of  several  hours,  and  they  could  be  imaged  under  a 
wide  vanety  of  bias  conditions  with  little  effea  on  their  appearance. 

One  of  the  most  important  issues  in  colloid  science  concerns 
the  nature  of  the  electric  double  layer  which  provides  the  repulsive 
Coulombic  barrier  that  keeps  the  panicles  in  suspension.  We 
hoped  to  use  the  STM  to  study  the  local  charaaer  of  the  double 
layer  by  aggregating  the  colloids  with  different  flocculating  agents. 
For  example,  lomc  salts  are  known  to  induce  aggregation  by 
causing  the  collapse  of  the  diffuse  double  layer  surrounding  each 
colloidal  panicle,  while  molecular  adsorbates  are  thought  to  cause 
aggr^ation  by  displacing  more  weakly  bound  charged  speaes  from 
the  colloid  surface  These  different  agents  should  produce  colloidal 
surfaces  with  significantly  different  electrical  and  molecular 
characteristics. 

Though  flocculating  with  ionic  salts  produced  aggregates  with 
excellent  tunneling  charaaenstics.  we  were  not  as  successful  with 
the  molecular  flocculating  agent,  teirathiafulvalene  (TTF)  ”  In 
fact  we  never  obtained  stable  STM  images  from  colloids  floccu¬ 
lated  with  TTF.  even  though  TEM  micrographs  show  that  these 
aggregates  are  identical  with  those  aggregated  with  salt.  We 
speculate  that  adsorbed  TTF  forms  an  insulating  bamer  around 
the  colloid  surface  that  prevents  tunneling  from  occurring.  Thus, 
to  obtain  the  best  STM  images  of  flocculated  colloids,  it  seems 
necessary  to  aggregate  with  ionic  agents  which  screen  the  Cou- 
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Hfure  2  STM  ind  TEM  imaics  of  colloidal  aold  aggregates  la,  tor* 
^TM  image  on  graphite',  ib.  boitomi  TEM  image  on  caroon-coateo 
copper  grid  The  colloid  Jizes  were  lound  to  be  ~  I  '  nm  bv  both  tech 

•'iqucc 

ombic  interaction  between  particles  without  hinderine  the  lun- 
'’CiMC  r'.'cess 

U  hiie  explorine  the  properties  oi  colloidal  aggregates  we  often 
■’oticea  trim  the  STM  images  contameO  sharp  spihes  Believing 
tnese  were  oue  to  the  presence  oi  codeposited  salts,  we  treated 
:nc  aggreeatea  suspensions  with  ion-e»change  resin  ‘  to  reduce 
the  concentration  ol  ions  in  solution  higure  '  displass  a  re 
"i.iric.]fiic  “sTM  image  oDiained  aiic  me  aggregated  conoids  were 
•^eaico  with  lon-excnanee  resm  T'ts  image  consists  oi  rough's 
so  itoms  which  have  avsumed  .i  nishis  ordered  ncxagonai  arr.is 
:nj:  measures  '  'fl  nm-  *  '  S(i  nm-  Ljch  ol  the  "atoms’  ir  the 
.litice  n.is  a  oiameter  oi  ^  '4  nm  .mo  ,i  iwivaimensinnai  Fourie' 
■  '  ,nsi,s'"s  .nf  'atiice  C'-'',;.':','  t  ■  nc  '4  nm. 

.oincc  prcsious  wor»  had  demonstrated  that  charged  colloidu. 
p.irticics  deposited  Irom  soiuiuir'  onti'  smooin  surijces  c.tr 
sometimes  crvstallize  in  hexagonal  lattices  '  we  originalis  asso- 
.  .tied  these  STM  images  with  our  sots  with  each  ".t'om  'dentilin.' 

,  0,.  g  ii’.  r.t't.^.e  IL.'wese-  ^c'.cr.i,  ^ n.ir.iv tcmstiss  ,c 

.i",!ss  .rgue  against  their  being  lattices  ot  colloid.’.'  gout 
iM'iiv.c.  !  "st,  ;r,c  individual  n,r'..^.es  .irc  extremcis  umiorm 
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Figat*  3  STM  image  of  anomalous  long-range  order  on  graphite  The 
pertodicitv  of  the  hexagonal  arrav  was  ~  .s  nm  and  it  extended  over  ar 
area  ol  '  I  50  x  |  50  nnr  T  ne  arravs  were  aiwavs  oDservea  near 
terrace  like  the  one  seen  in  the  lower  part  of  the  ftgure 

in  size,  much  more  uniform  than  isolated  particles  appear  wnen 
seen  by  either  TEM  or  STM  .Also  the  particles  in  the  lattice  oo 
not  conduct  as  much  current  as  isolated  gold  particles  do  under 
identical  expenmenial  conditions  Finally,  images  quaiitaiiveiv 
similar  to  those  in  Figure  3  were  obtained  alter  the  deposition 
of  colloidal  dispersions  which  were  not  treated  with  lon-excnange 
resin  Lnder  these  circumstances  the  tonic  strength  of  the  soi  is 
high,  and  the  long-ranged  Coulombic  forces  necessary  to  induce 
two-dimensional  crystallization  should  not  be  present 

We  speculate  that  this  highly  regular  lattice  is  associated  with 
the  graphite  substrate  itself  Indeed  similar  images  nave  oeen 
obtained  from  graphite  surfaces  after  the  deposition  ol  organic 
salts  from  solutions  which  contained  no  dispersed  pariicies 
Moreover,  verv  recent  studies  of  graphite  surfaces  bv  STM  have 
appeared,  and  almost  identical  images  were  obtained  wnere  nc 
deliberate  adsorption  of  material  occurred  '*  Kuwaoara  ei  a:  ' 
argue  that  the  anomalous  long-range  order  seen  on  graphite  car 
be  attributed  to  Moire  interference  between  two  siiehtiv  mtsv'- 
'icnteo  lasers  oi  crapnue  However,  is  possiDie  'r.i'.  severe 
mecnanisms  are  operating  simuiiancousiv  to  give  rise  to  me 
long-range  order  on  graphite  indeed,  prcliminarv  compute- 
simulations  snow  mat  anomalous  uint-range  order  car.  dcveior 
when  a  grapniie  llaxe  picKco  up  d\  the  nr  slides  acros,s  the  grap.nr.e 
substrate 

In  conclusion,  we  have  snown  mat  ;he  STM  can  oc  used  wnr 
some  conlioence  m  trie  studs  o!  trie  properties  oi  meta,  ..oiioios 
The  phvsicai  dimensions  ol  colloidal  panicles  can  be  measured 
with  rcasonaole  accuraev  m  me  plane  oi  the  substrate,  anc  laree 
aggregates  can  Pc  rcadilv  imaees  so  long  as  eieancailv  msuiannc 
■noicci'ics  arc  not  adsorbed  I’n  the  colloid  sur'avC  I  naer  some 
conditions  It  IS  possible  to  observe  hiehiv  ordered  hcxaconai  ,i-rjvv 
■n  crapnitc  sunaccs  Mihoucn  me  (’’■ipt  oi  ihis  arrav  :v  no:  we 
unocrsiiXH)  we  believe  n  i-  a  pncnomcnon  asscK-'iaico  w:ih 
crapti'.tc  suriacf  ■.isc'.l  >,ncc  me  la'.nce  p.mame’.cr  oi  the  orae-tc 
orav  comparable  to  me  size  ol  colloidal  nanicies  i;  :s  imrx’-.a-' 

cxcruisc  ^j:c  wnen  usinki  me  STM  to  studs  vO'ioiU.s  or  oir.e- 
njcromoiecuiar  species  adsorbed  on  grapniie  surlasCs 
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A  model  for  the  scanning  tunneling  microacope  (STM)  experuneni  is  presented  which  explicitly  treats  the  STM  tip  as  a 
polyatomic  crystalline  surface.  The  net  tunnehng  cuirent  is  the  sum  of  ail  tunneling  occurrences  from  each  tip  atom  to  every  substrate 
atom.  With  the  aid  of  computer  stmulauons.  the  atomic  resolution  of  an  hexagonally  close-packed  substrate  with  a  nearest  neighbor 
distance  of  0.24S  nm  is  reproduced  using  an  hexagonally  close-packed  up  with  a  nearest  neighbor  distance  of  0.277  nm.  The  resulting 
penodiaty  of  the  image  is  independent  of  tip  sue.  The  recent  observauons  of  anomalous  long  range  penodiaty  is  explained  by 
employing  the  model  in  a  sunulauon  of  a  2.7  nm  super  laiuce. 


Since  its  introduction  [1.2].  scanning  tunneling 
microscopy  (STM)  has  found  immediate  accep¬ 
tance  as  a  powerful  tool  for  surface  analysis.  Cur¬ 
rent  and  future  applications  of  this  teclmology  is 
highly  dependent  on  our  understanding  of  the 
measurement  technique.  This  Letter  presents  a 
new  model  for  describing  the  STM  experiment 
based  on  mutual  scanning  between  the  tip  and 
surface.  This  model  of  the  STM  experiment  is 
shown  to  reproduce  atomic  resolution  of  ordered 
surfaces,  as  well  as  provide  a  general  explanation 
of  anomalous  long  range  perio^city  that  has  been 
recently  observed  {3-5). 

As  a  first  approximation,  it  is  generally  be¬ 
lieved  |6)  that  the  STM  current  image  results  from 
an  electron  of  a  lone  tip  atom  tunneling  to  a  single 
surface  atom  when  an  external  bias  potenual  of 
appropriate  sign  is  applied.  The  exponential  de¬ 
pendence  of  tunnehng  current  on  both  interatomic 
distance  and  barrier  height  of  the  tunneling  elec¬ 
tron  selects  for  a  single  dominant  tunneling  path¬ 
way.  This  results  in  the  abihtv  to  map  out  the 
electron  orbital  charge  density  of  the  surface  as  a 
function  of  tunneling  current  when  a  tip  is  scanned 
parallel  to  the  surface.  This  mechanism,  although 
It  correctly  explains  many  phenomena  assoaated 


with  STM  images,  does  not  predict  the  long  range 
periodic  structures  that  are  often  observed  (3-5). 

We  propose  a  new  interpretation  of  the  STM 
experiment  for  tunnehng  to  an  ordered  surface 
such  as  highly  oriented  pyrolytic  graphite  (HOPG) 
from  a  metal  tip  having  a  close-packed  structure. 
The  proposed  mechanism  is  predicated  on  the 
following  STM  experimental  model.  Firstly,  that 
microscopically,  a  pure  metal  tip  is  crystalhne 
with  the  (111)  surface  dominating  in  many  cases. 
Secondly,  that  it  is  possible  that  luimehng  can 
occur  from  every  tip  atom  to  every  surface  atom. 
The  net  tunnehng  current  is  initially  assumed  to 
be  a  simple  summation  of  all  individual  tuimeling 
occurrences  which  are  assumed  to  be  independent 
and  behave  according  to  the  Tcrsoff  and  Hamann 
equations  (7-9).  Finally,  defects  and  step  edges 
that  posses  lower  barrier  heights  than  their  corre¬ 
sponding  bulk  surfaces  will  appear  as  adatoms 
and  contribute  umquely  to  the  resulting  image. 

The  modeling  of  a  sharp  metal  tip  as  an  ordered 
crystal  surface  is  a  logical  assumption  based  on 
the  following  developments,  lijima  et  al.  (10)  have 
observed  large  metal  clusters  containing  over  1000 
atoms  which  mamtain  single  crystal  structure  with 
a  tendency  to  maximize  the  (111)  surface  area 
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under  TEM  experimental  conditions  (1.3  x  105 
electrons  per  in^trdm  per  second).  Clusters  of  as 
few  as  500  atoms  experienced  frequent  rearrange¬ 
ment  to  new  single  crystal  structures.  This  work 
suggests  that  even  under  highly  energetic  condi¬ 
tions  large  close-packed  met^  clusters  are  crystal¬ 
line.  This  work  has  been  confirmed  by  Malm  et  al. 
[11].  Trevor  et  al.  [12]  observed  that  metal  surfaces 
which  have  undergone  electrochemical  oxidation/ 
reduction  to  tome  STM  tip  preparations 

tend  to  minimize  defects,  holes,  and  step  edges 
preferring  to  maxiimie  crystal  plane  areas.  Per¬ 
haps  the  most  compelling  evidence  for  a  planer  tip 
comes  from  work  done  by  Colton  et  al.  [13].  This 
group  showed  that  graphite  pencil  lead  or  graphite 
coated  tungsten  could  be  used  as  an  STM  tip  to 
reprodudbiy  image  HOPG  with  atomic  resolution. 
Finally,  the  STM  experiment  itself  has  failed  to 
show  that  pure  metal  surfaces  contain  stable 
adatoms  under  STM  experimental  conditions.  This 
conclusion  is  based  on  the  lack  of  pubhshed  data 
supporung  the  presence  of  adatoms  on  pure  metal 
surfaces  under  STM  conditions.  It  is  therefore 
logical  that  the  tip  be  treated  on  the  same  footing 
as  any  other  close-packed  metal.  In  fact  the  dis- 
imction  between  tip  and  substrate  is  purely  con¬ 
figurational  having  no  spedal  significance  when 
considering  the  relative  motion  and  interactions 
bMween  the  two  surfaces. 

^.The  net  tunneling  current  from  the  (1 11)  surface 
of' a  metal  tip  to  a  coplaner  (111)  surface  of  a 
metal  substrate  is  given  by 

y).  (1) 

•  J 

1,  (1.  j)’^A  exp|-fl[(jc,^-x,*)‘-t-(y,,- v,*)‘ 

•»!  1/21 

]■.  (2) 

where  eq.  (1)  is  a  simple  summation  over  all 
tip/ surface  tunneling  currents  and  (2)  is  Tersoff 
and  Hansma’s  familiar  result  [7,9].  As  the  up  is 
moved  to  a  new  position  this  current  vanes.  Since 
both  surfaces,  which  can  iniually  be  considered 
infuiitely  extended  in  parallel  planes,  inherently 
posses  translauonal  symmetry  in  both  the  x  and  \ 
direcuons.  the  net  current  is  penodic  havmg  a 
penod  related  to  the  smaller  atomic  spacings. 


This  can  be  exemplified  using  the  following 
gedanken  experiment.  Consider  two  coplaner  in¬ 
finite  lattices  that  are  superimposed  on  other 
with  one  lattice  translated  with  respect  to  the 
other  in  the  lattice  plane.  An  observer  could  then 
record,  at  some  finite  interval  in  the  plane,  a 
unique  set  of  images  for  each  of  the  configurations 
as  one  lattice  is  tianslatoi  itilliMhe  omiespond- 
ing  to  the  period  of  the  spacing. 

Due  to  the  translational  qpMUpjF  of  the  two 
lattices,  each  subsequent  m>vt  Images  would  be 
identical  to  the  initial  set  la  fact  since  the  direc¬ 
tion  of  the  translatioo  is  arbitrary,  only  the  first 
half  of  the  images  of  the  initial  set  will  be  unique. 

Consider  next  the  STM  experiment.  It  is  known 
that  as  a  crystal  surface  of  a  tip  is  scanned  across 
a  crystal  surface  of  a  substrate,  the  periodicity 
corresponding  to  the  smaller  atomic  periodicity 
(substrate  or  tip)  will  be  imaged.  Furthermore,  it 
is  known  that  HOPG  contains  two  inequivalent 
carbon  sites  [14]  resulting  in  the  STM  image 
dominated  by  the  tunneling  contribution  of  one  of 
these  atoms.  Fig.  1  reports  the  results  of  a  series  of 
two  dimensicmal  simulations  where  a  multi-atom 
tip  (with  the  nearest  neighbor  distance  0.2T7  nm 
corresponding  to  platinum)  has  scanned  a  25  atom 
substrate  (with  nearest  neighbor  distance  of  0.245 
am  corresponding  to  the  beta-beta  carbon  dis¬ 
tances  in  HOPG).  Within  the  constraints  of  the 
independent  multiple  tunnelmg  assumpuon.  de¬ 
scribed  above,  it  is  evident  that  multiple  atom  tips 
do  not  necessarily  yield  multiple  images.  It  is  clear 
from  these  two  dimensional  cross  secuons  that  the 
multi-atom  tip,  with  a  larger  rnteratomic  distance, 
correctly  reproduces  the  simulated  current  image 
of  HOPG.  The  simulauon  was  simplified  usmg  a 
straight  forward  reduction  of  vanables  [16]  and 
performed  in  a  constant  height  mode  in  keeping 
with  recent  experimental  observauons  [5],  Fig.  1 
shows  that  neither  the  atomic  periodicity  nor  the 
corrugauon  amphtude  is  effected  bv  tip  size. 

If  the  surface,  with  the  smallest  atomic  spacmg. 
contams  a  defect  which  results  in  a  lower  bamer 
height  than  that  of  the  bulk  surface  [IS],  a  long 
range  hexagonally  close-packed  image  with  a  pen- 
odiaty  much  larger  than  the  substrate  is  observed. 
This  results  from  the  surface-to-defect  tunneling 
path  of  the  large  atomically  spaced  surface  to  the 
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Fig.  1 .  A  (w<Miuiiensioaal  crou  mcuoo  o(  a  siroulated  STM  consuat  faaghi  unage  as  a  funcuon  of  up  size.  The  number  of  Up  auans. 
which  are  scparaind  by  0.2T7  am.  u  indicaied  next  lo  each  curve.  The  subttfate  is  repreaeaied  by  25  atoms  separated  by  0.245  mn. 
Note  that  a  multiple  atom  up  does  not  necessarily  produce  a  muluple  image  under  this  independent  lunnelmg  assumpuoa.  The  edge 

disioruons  are  caused  by  over  scanning  the  substrate  surface. 


surface  defect  of  the  smaller  atomically  spaced 
surface.  In  essence,  the  defect  scans  the  larger 
latuce  of  the  surface  as  the  tip  translates  across 
the  substrate.  Therefore,  the  image  of  the  tip  is 
added  to  the  net  current  image  as  produced  by  the 
defect.  Since  the  defect  must  result  in  a  lower 
barrier  height  in  order  to  act  as  a  pseudo-tip.  it 
contributes  significantly  to  the  net  current  image. 
Thus,  the  resulting  STM  image  is  the  superposition 
of  both  the  tip  and  substrate  images.  Fig.  2  ex- 


hibiu  a  simulation  of  the  image  due  to  a  3S  x  IS 
atom  hexagonally  close-packed  tip  emnning  a  3S 
X  IS  hexagonally  close-packed  substrate  which 
contains  a  smgle  defect. 

The  superposition  of  the  tip  and  substrate  image 
can  result  in  a  variety  of  large,  long  range  peri¬ 
odicities.  This  is  because  the  xv-plane  of  the  tip 
will  be  onented  at  an  angle  v  (0<v<  30® )  rela¬ 
tive  to  the  xy-plane  of  the  surface.  The  intensity  of 
the  resulting  image  vanes  as  a  function  of  the 
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Lgi«i»l|B  pcnodidty  liiiwlilil  by  the  MpcnxMitiati  of  •  (III)  (cc  lurface  with  ■  nearest  neighbor  distance  of  0.277  nm 
which  n  oa^HMK  and  rotated  by  a  -  0*  with  reipeCT  to  an  hexagonai  latuce  «nih  a  nearest  neighbor  distance  of  0.24S  nm.  The 
image  m  MS  by  I.SS  nm.  and  dearly  shows  a  long  range  penodiatv  of  approumaielv  2.87  nm.  (The  apparent  compression  of  the 
hmgonai  stnicuire  m  the  image  is  due  lo  the  graphical  output  device  and  u  not  a  feature  of  the  dau. ) 


phase  of  the  two  images,  which  is  a  function  of 
their  atomic  surface  periodicities  as  well  as  the 
angle  rj  between  them. 

Long  range  periodic  structures  on  graphite  have 
been  previously  reported  (3.4].  In  many  cases, 
these  structures  have  been  observed  to  be  associ¬ 
ated  with  graphite  defects  [3-Sj.  In  our  experience 
ihev  usually  can  be  assoaated  with  basal  plane 
step  edges  (5).  The  carbon  atoms  lying  at  step 
edges  are  unique  with  respect  to  bulk  graphite 
carbon  atoms  due  to  their  unfilled  (open-shell) 
atomic  orbital.  If  those  step  edge  atoms  have 


significantly  different  ionization  potentials  in 
companson  to  the  defect  free  graphite,  they  should 
logically  constitute  low  barrier  height  defects. 
Hence,  what  appears  in  some  STM  images  to  be 
long  range  periodicity  can  be  produced  by  point 
defects  located  at  basal  plane  step  edges. 

Although  we  have  chosen  to  investigate  the 
speaai  case  of  a  coplanar  tip  and  substrate,  we  do 
not  wish  to  imply  that  this  is  representative  of  the 
STM  expenment  in  general.  Clearly  the  effects  of 
non-coplanantv  need  to  be  investigated  in  a  more 
thorough  studv.  However,  the  reader  is  reminded 
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that  the  aDomalous  long  nuife  oocnigauons  such 
as  those  reported  to  be  modaied  with  HOPG 
|3.S]  are  also  speciai  cue  ieMgM  which  are  not 
represenutive  <k  the  typical  graphite  surface.  Fur¬ 
thermore.  the  lastrictive  coaditioe  that  the  STM 
tip  contain  a  1bi|c  planar  surface  can  be  lessened 
by  recognizing  that  STM  lips  have  been  shown  to 
adsorb  large  HOPG  flakes  [17]  which  are  inher¬ 
ently  planar. 

We  do  not  ilisagerr  that  many  STM  expen- 
menu  can  be  Modeled  by  assumiiig  a  single  atom 
tip.  This  is  reaeonable  since  a  aiicroecopically 
crystalline  tip  wiO  have  edges  and  ooners  awl  it  is 
unlikely  that  the  substrate  will  be  anpt****  with  a 
crysul  face  of  the  tip.  However,  sinoe  most  tips 
are  made  of  bulk  nutterial  (tangBian  or  platinum 
wire,  for  example)  and  upon  initial  currant  flow 
there  exisu  sufficient  energy,  possibly  due  to  resis¬ 
tive  heating,  for  the  tip  to  anneal  to  its  lowest 
energy  interatoinic  orientation.  It  is  still  unclear  if 
a  planar  tip  will  be  unable  to  reproduce  sharp  step 
images  as  some  have  suggested.  It  is  also  unclear 
to  what  extent  the  ooplanar  oonfignratkm  can  be 
relaxed  such  that  long  range  periodicity  can  still 
be  produced.  However,  we  have  shoim  that  a 
multi-atom  crystal  domain  STM  tip  can  explain 
both  atomic  pmiodicity  and  the  anomalous  long 
range  order  often  observed  on  graphite. 

At 

This  naaaich  is  supported  in  part  through  the 
Air  force  OflUce  of  Sdrotific  Research  and  by  the 
National  ScienPf  Foundation  under  grant  No. 
CHE-8912674.  We  are  also  grateful  to  Professor 
D.M.  Kalyon  and  Dr.  J.  W'ang  for  access  to  the 
graphical  facilities  used  in  preparauon  of  this 
manuscript. 


Referenees 


(11  G.  Biaius.  H.  Rohrer.  Ch  Gerber  and  E  Weibel.  Phvs 
Rev.  Lett  49  (1982)  37 

|2]  G.  Buiait.  H.  Rohrer.  Ch.  Gerber  and  E  Weibel.  AppI 
Phya.  Lett  40  (1982)  178 

(3)  M.  Kuwabwa.  D.R.  Clarke  and  D.A.  Smith.  AppI.  Phvs 


Lett  36  (1990)  2396. 

(4)  J.W.  Lydias. 

Brockaabraush.  J.I 
TechnoL  A  6  (1918)163. 

(31  J.F.  Womatadorf.  ftSCl 
iitbed. 

|6)  RJ.  Haaen.  Aaan.  hcv.  PhyalSMlL’^i  11989)  331. 

[7]  3.  TeraoTf  and  D.R.\HaaaBa.  Asa.  Mm.  Lau  30  (1983) 
1998 


S.  Skala.  R 
J.  Vac.  Sci 


oapub- 


(8)  J.  Tenoff  and  D.R.  Hamann.  Phvs.  Rev.  B  31  (1985)  803 
(91  N.D.  Lang,  Pbyi.  Rev.  B.  34  (1986)  3947 

(10]  S.  lijtma.  in:  Micraduften.  Eds.  S.  Sugano.  Y  Nishina 
and  S.  Ohnishi  (Sprager.  Berlin.  1987)  p  186 

(111  J.-O.  Mabn  and  J.-O.  Bovin.  Surf.  Sci  200  (1988)  67 

(12)  DJ.  Trevor.  C.E.D.  Chidsev  and  D  N  Loiacono.  Phvs 
Rev.  Leu.  62  (1989)  929 

(13|  RJ.  Colton.  S.M.  Baker.  RJ.  Dnscoll.  M.G  Youngquisi 
J.D.  BaldeiCbwifler  and  W.J  Rauer.  J  Vac.  Sci  Technol 
A  6  (1988)  349. 

(14)  G.  Biania.  H.  Fuchs,  Ch.  Gerber.  H.  Rohrer.  E  Stoll  and 
E  Toaatit.  Europbys.  Lett.  1  (1986)  31 

(13)  B.  Kxahl-Urbaa.  H.  Wapier  and  R.  Buiz.  Surf  Sci  93 
(1980)  423. 

(16]  The  reducuon  of  variables  was  accomplished  as  follows 
The  duiance  term  tn  the  exponent  was  expressed  m  Cane- 
Stan  coordinaiea.  Since  the  simulauon  was  performed  in 
constant  hei|ht  mode,  the  constants  in  the  exponential 
(haghL  work  function,  etc  i  were  grouped  lofether  and 
set  arbitranly  equal  to  one  as  was  the  pre.cxponeniial 
factor.  This  simphficauon  was  etnploved  since  the  lou 
bias  lunnebng  rdationahips.  previously  denved  bv  Tersoff 
and  Hamann  iref.  |8.10|)  are  assumed  to  be  correct  It  ivav 
found  that  assigning  to  the  exponential  constant  values 
between  0  and  100  simpiv  scales  the  c<omponeni  ( i  c 
ampUiude)  of  the  image  W’e  found  this  simplification 
acceptable  unce  corrugation  amplitude  is  not  the  locus  o( 
this  siudv 

(17)  T  Tiedie.  J  Varon  and  J  Siokes.  J  tac  Sci  Technol  ^  ^ 
(1988)  37: 


Ab  Zaitio  study  of  tho  Bloetronio  Spootrua 
of  a  ^*1x^3  Clustor 

Ping  Wang 

Dopartaont  of  Chaaistry  and  Cheaical 
Stavans  Instltuta  of 

Hobokan,  Naw  Jersey  07030 
and 

Walter  C.  Ermler 
Dapartaent  of  Chemistry 
The  Ohio  State  University 
Coltiabus,  Ohio  43210 

ABSTRACT 

Caoaatries  and  electronic  states  of  the  cesium  suboxide 
cluster  are  investigated  using  ab  initio  Hartree-Fock 

calculations.  Relativistic  effective  core  potentials  are 
employed  to  represent  core  electrons  in  the  atoms  while 
Gaussian-type  orbital  basis  sets  are  used  to  describe  the 
valence  electrons.  Interatomic  distances  are  optimized  for 
the  ground  electronic  state.  The  workfunction  of  is 

calculated  as  0.65  eV,  in  excellent  agreement  with  an 
experimentally  determined  value  of  0.7  Ev  for  bulk  cesium 
suboxide.  Mulliken  populations  of  negatively  charged  and 
neutral  clusters  have  been  calculated  and  indicate  that 

upon  ionization  the  emitted  electron  exits  a  molecular  orbital 
encompassing  the  outermost  cesium  layer  of  the  cluster. 


♦Permanent  address  of  Walter  C.  Ermler 


•Ittctron  sp«otrosc(^>y  by  de-excitation  of  aetastable  noble 
gas  atoas  (MDS)  to  follow  the  continuous  oxidation  of  cesiun 
filas.®  HU.S  group  applied  UPS  for  a  finite  inforaation 
depth  and  MDS  for  the  properties  of  the  outeraoM  atonic 
layer.  Ihe  esqperinental  results  show  that  —all  of  O2 


CS11O3  is  decreased  in  craparison  to  pure  Cs 


cause  the  fomation  of  a  "monoplane"  of  CS33O3 
ions  are  incorporated  below  the  surface,  sfeOLle  tiM  surface 
itself  consists  solely  of  Cs  atoms.  The  workfunction  of 

ital.  A  theory 

of  the  vorkfunction  of  cesium  suboxide  and  a  model  of  its 
electronic  structure  %ras  proposed  by  Burt  and  Heine^  based 
on  Simon's  studies  of  the  crystallography  of  cesium 
suboxide.  They  determined  that  a  substantial  drop  in 
morkfunetion  is  to  be  expected  in  going  from  Cs  to  08^^303 
$  dms^w  changes  in  boundary  conditions  on  the  conduction 
^  eXedIrons. 


The  present  work  investigates  a  cesium  suboxide  cluster 
CS33^03  using  ab  initio  calculations.  The  geometry  of  the 
cluster  in  low-lying  electronic  states  is  studied.  The 
workfunction  of  the  cluster  is  calculated  and  Mulliken 
populations  are  analyzed. 


II.  CALCOZATIOMS 

The  average  relativistic  effective  potentials  (AREP) 
for  oxygen  were  obtained  from  Pacios  and  Christiansen.^^ 
They  used  a  refined  version  of  the  Dirac-Fock  based 
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relativistic  effective  core  potential  procediure  of  Lee, 
Ennler  and  Pitzer^^  to  conpute  AREP's  of  Li  through  Ar.  They 
also  reported  snail  optimized  Gaussian-type  orbital  (6T0) 
basis  sets  for  the  valence  electrons  of  Li  through  Ar.  The 
AREP  and  Gaussian-type  basis  set  used  for  cesium  were  taken 
from  Ross  et  al.^^ 

2 -electron  core  and  6-electron  valence  spaces  for 
oxygen  and  46-electron  core  and  9-electron  valence  spaces 
for  cesium  were  chosen.  In  cesium,  the  five  s-type  and  five 
p-type  GTO's  are  contracted  to  three  s-  and  two  p-typm 
functions,  respectively.  f 

In  the  RHF  calculations,  one-  and  two-electron  malbrix 
elements  involving  symmetry-adapted  linear  combinations  of 
GTO's  were  calculated  using  "equal  contribution"  integral 
program  of  Pitzer  and  Hsu.^^  The  wavef\inctions  and  energies 
were  calculated  using  a  restricted-Hartree-Fock  self- 
consistent-  field  program.^® 

Geometry  optimizations  were  performed  beginning  with 
the  crystal  structure  of  CSj^j03^®  using  a  one-dimensional 
"breathing"  model.  The  equilibrium  geometry  was  then 
obtainmi^lifnii  a  plot  of  energy  versus  breathing  radius.  The 
geosetrii^^^mtmmeters  associated  with  the  minimum  energy  for 
this  curve  were  then  used  in  RHF  calculations  to  obtain  the 
corresponding  energy. 

Electronic  state  optimizations  were  performed  for 
anionic,  neutral  and  cationic  clusters  and  the  corresponding 
workf unctions  were  extracted. 
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III.  RESULTS 

Table  I  shows  the  coordination  numbers  of  the  Cs^^03 
cluster.  The  structure  of  the  CS]^303  cluster  can  be 
described  as  being  formed  by  layers  of  Cs  and  O  atoms  (see 
Fig.l).  The  upper  layer  consists  3  CsB  atoms  possessing  a  C3 
rotation  axis.  The  middle  layer  consists  of  3  OA  and  3  CsA 
atoms.  The  lower  layer  is  structurally  similar  to  the  upper 
layer  and  spaced  from  the  middle  layer  by  the  sauae  distance 
as  the  upper  layer.  The  two  CsC  atoms  are  nested  somewhat 
closer  to  the  middle  layer  than  the  upper  and  lower  layers 
and  lie  along  the  C3  axis. 

The  abovementioned  AREP's  and  contracted  GTO  basis  sets 
lead  to  Hartree>Fock  energies  of  -19.842X8  and  -15.64121 
hartrees  for  the  cesium  and  oxygen  atoms,  respectivelyj^^ 
Table  II  shows  the  structural  pareuaeters  of  CS3]^03 
obtained  experimentally  by  Simon  and  Westerbeck. The 
cluster  geometry  was  optimized  and  the  effected  of  geometric 
changes  on  total  energies  are  shown  in  Table  III.  The 
geometric  parameters  associated  with  the  minimum  energy  for 
this  co^rMpond  to  a  7%  expansion  of  the  crystal  structure 
of  Table  2. 

Tables  IV-VI  show  total  valence  energies  of  several 
low-lying  electronic  states  for  anionic,  neutral  and 
cationic  clusters,  respectively.  All  calculations  correspond 
to  the  average  energy  of  molecular  orbital  configurations.^^ 
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We  chose  to  treat  such  averaged  spin  states  because  the  size 
of  the  cluster  (11  Cs,  30  and  117  valence  electrons) 
results  in  severe  computational  complexity.  All  electronic 
state  optimizations  were  performed  at  the  7%  swelled 
equilibrium  geometry  of  TaUsle  III  and  the  lowest  states  were 
found  for  the  negative,  neutral  and  positive  clusters.  It  is 
noteworthy  that  the  negatively  charged  cluster  is  found  to 
be  the  most  stable. 

Mulliken  populations  on  each  of  the  atoms  in  the  ground 
states  of  the  negative,  neutral  and  positive  clusters  are 
listed  on  Table  Vll.  By  comparing  the  electron  populations 
on  each  type  of  atom  in  the  negative  and  neutral  clusters, 
it  is  clearly  seen  that  the  ionized  electron  comes  from  the 
outermost  cesium  layer. 

The  energies  of  the  highest  occupied  molecular  orbitals 
(HOMO)  and  lowest  unoccupied  MO's  (LUMO)  are  listed  in  Table 
VIII.  The  energies  of  the  HOMO'S  of  all  the  clusters  are 
negative.  However,  the  energies  of  the  LUMO's  for  the 
positive  cluster  are  all  negative,  those  for  neutral  cluster 
being  partially  negative  and  partially  positive,  and  those 
for  ne^^illpiii^  cluster  all  positive.  These  indicate  that  both 
positisiiilil|||||^  neutral  clusters  are  electron  deficient  and  the 
negative  cluster  is  the  most  stable.  This  conclusion  is 
reinforced  by  comparing  their  total  energies.  The 
workfunction  of  the  cluster,  obtained  by  subtracting  the 
total  energy  of  the  negative  cluster  from  that  of  neutral 
cluster,  is  0.651  eV.  This  may  be  compared  to  the 
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corresponding  value  of  0.700  eV  calculated  using  Koopnans^ 
theorem. 


IV.  DISCUSSION 

The  Cs^]^03  cluster  model  studied  is  shown  in  Fig.  l.  It 
has  the  layer  symmetry  of  the  bulk  crystal  and  has  four 
unique  atoms.  The  cluster  has  031^  point  group  symmetry  and 
represents  a  fragment  of  the  bulk  crystalline  solid,  '^e  use 
of  the  03j^  point  group  symmetry  and  averaged  relativistic 
effective  core  potential  to  replace  the  core  electrons  for 
CS11O3  led  to  substantial  reductions  of  the  numbers  of  two- 
electron  integrals  that  must  be  computed  and  processed  for 
each  SCF  calculation. 

The  contractions  of  the  basis  sets  for  cesium  and 
oxygen  atoms  leave  the  outermost  (smallest  exponent) 
primitive  functions  uncontracted,  making  single  contracted 
functions  of  the  innermost  (largest  exponent)  primitives. 
This  permits  the  linear  SCF-MO  coefficient  of  the  outer 
function,  having  maximum  amplitude  in  the  interatomic 
region r respond  to  the  changes  which  occur  upon  cluster 
formatlHiii; " 

The  interatomic  distances  were  initially  determined 
from  the  lattice  constants  of  crystalline  Csj^203^^  (see 
Table  II)  and  the  geometry  optimization  was  carried  out 
using  a  breathing  model.  The  results  appearing  in  Table  3 
indicate  that  the  total  energy  of  the  cluster  reaches 
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minlMia  at  7%  expansion  of  the  crystalline  size  and  the 
naked  cluster  is  stable  in  the  gas  phase.  The  total  energy 
of  the  cluster  is  lower  by  81.5  kcal/nole  than  that  due  to 
the  summation  of  the  energies  of  the  separated  atoms. 

The  electronic  state  optimizations  were  carried  out  for 
negative,  neutral  and  positive  clusters.  The  procedures  are 
similar  to  those  following  by  Ermler  et  al.  in  their  studies 
of  the  properties  of  large  Be  clusters. The  results 
appearing  in  T2d)le  IV-VI  show  that  the  most  stable 
electronic  state  has  an  open  shell  configuration  and  that 
the  negative  CS]^]^02~  cluster  is  the  most  stable.  The 
behavior  of  the  lowest  electronic  state  is  similar  to  t^t 

i  1 

of  Be^3^^,  and  is  indicative  the  onset  of  the  bulk 
characteristics,  i.e.,  the  presence  of  a  partially  filled 
band  in  which  spin  couplings  are  not  important  energetically 
and  electrons  are  highly  delocalized.  This  provides  some 
justification  for  carrying  out  electronic  state 
optimizations  based  on  averaged  energy  of  configuration 
states  for  open  shell  cases. 

Tii%jDonclus ion  that  the  most  stable  cluster  is 

supported  by  HOMO  and  LUMO  analyses. 

As  VIII,  the  six  LUMO's  for  the  lowest-energy 

configuration  of  the  positive  cluster  Csj^2^03'’’  are  all 
negative  and  the  majority  of  the  LUMO's  for  the  neutral 
cluster  are  negative,  indicating  that  the  positive 

and  neutral  system  are  electron  deficient  for  the  low-lying 
states.  The  singly  charged  anionic  Cs^3^03~  was  then  studied 
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in  an  attempt  to  satisfy  this  electron  deficiency.  The 
results  appearing  in  the  Table  VIII  for  show  that 

its  six  HOMO'S  are  all  negative  and  all  of  its  LUMO's  are 
positive,  indicating  that  the  electrophilicity  of  the  system 
has  been  quenched  and  addition  of  one  electron  to  the 
neutral  cluster  increases  its  overall  stability.  These 
findings  are  similar  to  those  reported  by  Marino  et  al.  on 
their  studies  of  the  layered  semiconductor  cluster  Biglj^g.^^ 

Charge  per  atom  values  appearing  in  Table  VII  indicate 
that  the  predominant  interlayer  interactions  occur  between 
the  CSC  and  OA  atos  (Fig.  1)  for  the  lowest-energy  state  of 
CSii03~.  This  was  also  observed  for  the  other  five  low^lylng* 
states  of  CS2^]^03'.  However,  those  for  083^3^03  and  ^83^3^03'** 
indicate  that  predominant  interlayer  interactions  occur 
between  two  layer  six  CsB  atoms  and  three  OA  atoms  (see 
Fig.l).  By  comparing  charge  per  atom  values  in  the  negative 
and  neutral  clusters,  we  also  found  that  the  values  for 
oxygen  are  constant  at  -1.59  whereas  those  for  the  CsB  and 
CsC  atoms  are  changed  by  +0.34e  and  -0. 41e,  respectively.  A 
detailed  analysis  of  Mulliken  populations  reveals  the 
followim^t  (a)  In  ionization  from  the  negative  cluster 

the  neutral  cluster  Cs3^3^03,  the  electron  exits 
from  the  outermost  cesium  layer,  and  (b)  the  predominant 
interlayer  interactions  change  from  CsC-OA  to  CsB-OA. 

The  workfunction  of  the  cesium  suboxide  is  one  of  its 
most  important  properties.  Our  calculated  results  show  that 
CS11O3”  has  very  low  workfunction,  consistent  with  that 
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experinentally  for  cesium  suboxides. ^  Two  theoretical  values 
correspond  to  a  calculation  of  the  workfunction  as  the 
difference  between  total  eneirgy  of  CS|^2.^3~ 

Koopmans'  theorem.  In  the  latter  case,  two  inherent 
approximations  can  result  in  a  more  reasonable  value  of  the 
workfunction  for  systems  with  large  numbers  of  electrons: 

(a)  The  ** frozen  orbital"  approximation.  This  assumes  that 
the  spin  orbitals  for  the  N-electron  states  are  identical 
with  those  for  the  (N+1) -electron  states.  This  approximation 
neglects  relaxation  of  the  spin  orbitals  in  the  (N+1)- 
electron  state  to  produce  a  workfunction  that  is  too 
positive,  (b)  "constant  differential  correlation  effects" 
approximation.  The  single  determinant  vravef unction  leads  to 
well-know  errors  in  description  of  instantaneous  electron 
correlation  effects  and  one  must  go  beyond  the  Hartree-Fock 
approximation  in  order  to  recover  such  corrections  to 
Koopmans'  theorem.  In  particular,  correlation  energies  are 
larger  for  the  system  havig  more  electrons.  Therefore, 
correlation  effects  tends  to  cancel  the  relaxation  error  for 
the  workfunction.  Our  model  cluster  contains  118 

valence  electrons  and  the  calculated  workfunction  value 
using  Kocq^mans'  theorem  is  0.7  eV,  which  is  comparable  with 
that  of  0.65  eV  using  the  energy  difference  between  two 
Hartree-Fock  calculations-  By  comparing  our  calculated 
workfunction  values  with  experimental  value,  it  is  seen  that 
Koopmans'  theorem  gives  a  value  in  close  agreement  with 
experiment.  This  close  agreement  between  SCF,  Koopmans' 
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theorea  and  experimental  values  indicates  that  since  CS2^2.^3~ 
is  a  system  with  large  number  of  electrons,  Koopmans' 
theorem  is  applicable  and  may  be  expected  lead  to  an 
accurate  workf unction. 

Gregory  et  al.^  reported  the  workfunction  of  cesium 
suboxide  to  be  0.7  eV.  To  obtain  this  value,  they  studied 
the  oxidation  of  cesiiim  using  ultraviolet  photoemission 
spectroscopy.  They  carried  out  the  experiment  in  a 
stainless-steel  ultrahigh-vacuum  system  with  a  base  prtBSSure 
of  about  3xl0~^^  Torr  and  measured  the  distribution  curves 
(DEC)  using  a  standard  retarding  potential  system. They 
found  that  upon  exposure  of  a  fresh  cesium  film  to  oxygen, 
very  narrow  peak  which  grew  with  increasing  €uq;>osure 
appeared  in  the  DEC'S  about  2.6  eV  below  the  Fermi  level 
(Ef) .  This  peak  is  associate  with  oxygen  ions  dissolved 
below  the  surface  of  the  cesium  metal.  After  3xl0“®  Torr/ sec 
of  exposure,  additional  structure  began  to  appear  and  the 
structure  associated  with  the  oxides  changed  with  increasing 
oxygen  exposure.  At  500  L  oxygen  exposure  the  workf unction 
reached  a  sharp  minimum  of  0.7  eV,  which  was  obtained  from 
the  DE^f^4|g  subtracting  the  width  of  the  EDC  from  the 
photon  used  to  produce  the  DEC'S.  The  structure  of 

present  model  cluster  is  consistent  with  this  experimental 
observation,  viz.  oxygen  ions  dissolved  in  the  cesium  metal 
below  the  surface. 

Burt  and  Heine^  studied  cesium  suboxide  by  calculating 
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the  band  structure  using  a  layer  method.  Their  model  solid 
is  an  FCC  lattice  of  hard  spheres  with  a  constant  potential 
between,  a  lattice  constant  of  30  a.u  corresponding  to  the 
suboxide  was  used.  The  model  of  the  electronic 

structure  of  the  system  proposed  by  Burt  and  Heine  is 
composed  of  one  region  approximating  cesium  metal  with  other 
regions,  occupied  by  oxygen  ions;  the  latter  highly 
repulsive  to  the  conduction  electrons.  Their  calculated 
reduction  in  workf unction  compared  with  cesium  metal  for 
CSi^o^  is  0.9  eV.  This  result  is  in  general  agreement  with 
the  drop  in  workf unction  upon  oxidation  of  thick  cesium 
films  observed  by  Gregory  and  with  the  present  calculations. 
In  their  studies,  the  lower  workfunction  of  the  cesium 
suboxide  monolayer  film  compared  with  that  of  bulk  cesium 
was  explained  as  a  quantxim  size  effect.  That  is, 
confinement  of  an  electron  to  a  monolayer  raises  its  kinetic 
energy  and  hence  the  Fermi  energy,  leading  to  an  equivalent 
decrease  in  workfunction.  This  claim  is  further  verified  by 
our  calculations. 

Woratschek  et  al.®  presented  evidence  for  a  quantum 
size  efiillHbt  for  the  conduction  electrons  during  oxidation 
Cs.  ThllPgfoup  stT'.died  the  continuous  oxidation  of  a  Cs  film 
using  ultraviolet  photoelectron  spectroscopy  and  electron 
spectroscopy  by  deexcitation  of  metastable  noble  gas  atoms. 
They  found  that  the  transition  from  metallic  Cs  to  Csj^j^03 
leads  to  a  twofold  increase  of  the  intensity  of  electron 
emission  from  the  valence  band,  even  though  the 
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concentration  of  conduction  electrons  is  reduced  by  about 
50%.  This  is  attributed  to  a  quantiun  size  effect. They 
also  found  that  a  thin  layer  of  082^^03  was  created  by 
exposing  a  clean  Cs  surface  to  small  amount  of  03.  The 
resulting  worXftinctions  were  2.0  eV  for  Cs-poly  and  1.5  eV 
for  €833^03 -poly.  Although  the  workf unction  values  are 
different  from  those  observed  by  Gregory  et  al.,^  Burt  et 
al.,^  and  with  the  present  work,  the  decrease  in 
workfunction  compared  with  Cs-poly  for  CS33^03-poly  is  still 
in  general  agreement  with  the  present  calculations  and  other 
experiments. 

A  very  recent  study  by  Limberger  and  Martin^  ^  reported  ^ 
that  with  increasing  oxidation  of  the  cesium  clusters  the  IP 
decreases  and  reaches  a  minimum  at  Cs  (€820)^3.  They 
investigated  cesium  and  cesium  oxide  clusters  by  ionization 
in  a  one-photon  experiment  by  means  of  a  tunable  CW  dye 
laser  and  a  mass  spectrometric  detection  system.  The 
interpretation  of  the  ionization  threshold  was  given  in 
terms  of  a  displaced  harmonic  oscillator  model  which  leads 
to  a  reltfiKidn  between  adiabatic  and  vertical  ionization 
potentSH^inimir  findings  support  our  results  as  in  the 
folloiM|^HBPhiMPS  (i)  stability  of  the  negatively  charged 


cluster.  In  the  bulk,  they  found  that  the  oxygen  appears  as 
an  oxide  0^“,  a  peroxide  03^”,  or  a  superoxide  anion  02”. 
These  observations  indicate  that  the  negative  cluster 
031103“  most  stable,  in  agreement  with  our 

calculated  results.  (2)  Influence  of  the  number  (n)  of 
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cesitui  atoBs  per  cluster.  Their  experimental  results  show 
that  CSj^Ojg  cluster  with  n  an  odd  number  has  a  lower  IP  than 
those  with  n  even  because  those  with  odd  n  contain  an 
unpaired  electron.  This  finding  is  also  consistent  with  our 
model  which  includes  11  cesium  atoms  and  an  open-shell 
ground  state.  (3)  Influence  of  the  number  (m)  of  oxygen  atom 
per  cluster.  A  doped  oxide  cluster  C8(Cs20)j^  containing  one 
excess  cesiiui  has  the  lowest  IP.  Our  model  could  be  %nritten 
as  Cs5(Cs20)3  and  is  in  general  agreement  with  Cs(Cs20)|^ 
when  the  number  of  cesium  atoms  in  the  cluster  remains  odd. 
(4)  The  cluster  has  several  isomeric  forms,  each  with  its  | 
o%m  threshold.  This  may  explain  the  workfunction  diffsOrsnce 
between  our  calculated  value  for  Cs^^^^S  that  observed  by 
Woratschek  et  al.,  that  is,  €032^03  could  have  stable  isomers 
possessing  different  workf unctions. 
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V.  COHCLU8ZON8 

Our  findings  regarding  the  workfunction  of  Cs-^2P2~ 
provide  a  theoretical  basis  for  experimental  observations  of 
very  low  workfunctions  for  cesium  stiboxides.  The  consistency 
between  calculated  and  experimental  values  indicates  that 
our  model  provieds  a  good  representation  of  the  cesium 
suboxide  surface.  Mulliken  population  analysis  reveals  that 
the  ionized  electron  comes  from  the  outermost  cesium  layer. 
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Table  I.  Coordination  Numbers  of  the  Cs^^^O^  Cluster.^ 


Atom 

Number 

• 

CSA*^ 

3 

CsB° 

6 

CsC^^ 

2 

OA® 

3 

• 

^See  Fig.  1. 

^CsA  =  Csl;  2;  9. 

°CsB  »  Cs5;  6;  7;  8;  10;  11. 
^CsC  =  Cs3;  4. 

®OA  =01;  2;  3. 
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Table  II.  Geometric  Structure  of  the  Cluster 


Atom  Pair  Interatomic  Distance  (A) 


CS3 

-  CS4 

3.653(2) 

CS3 

-  Csl;  2;  9 

3.786(3);  3.783(3); 

3.766(2) 

CS4 

-  Csl;  2;  9 

3.754(3);  3.741(3); 

3.769(3) 

CS3 

-  Cs5;  7;  10 

4.165(3);  4.153(3); 

4.123(3) 

Cs4 

-  Cs6;  8;  11 

4.153(3) ;  4.156(3) ; 

4.185(3) 

CSl 

-  Cs5;  6;  10;  11 

4.244(4);  4.304(3); 

4.250(4)  ; 

4.344(3) 

CS2 

-  cs7;  8;  10;  11 

4.310(3) ;  4.327(3) ; 

4.311(4) ; 

4.285(4) 

CS9 

-  Cs5;  6;  7;  8 

4.275(4);  4.370(3); 

4.388(3)  ; 

4.310(3) 

Cs5 

-  Cs6 

4.022(3) 

Cs7 

-  Cs8 

4.103(4) 

CslO 

-  Csll 

4.063(4) 

01  - 

Csl;  2;  3;  4;  10;  11 

2.940(14) ;  2.914(15) 

;  2.970(14); 

2.957(13) ;  2.688(14) 

;  2.760(15) 

02  - 

3;  4;  7;  8;  9 

2.967(13);  2.981(14) 

;  2.964(14); 

7.731(14) ;  2.731(14) 

;  2.918(13) 

03  - 

csl;  3 ;  4 ;  5 ;  6 ;  9 

2.916(14);  2.982(14) 

;  2.921(14); 

2.684(15) ;  2.763(14) 

;  2.948(13) 

oi  - 

02;  3 

4.061(20) ;  4.021(18) 

02  - 

03 

4.037(20) 
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Teible  III.  Geometry  optimization  of  CSj^2^03 

Expansion  (percent  of  Total  Energy 

the  lattice  parameters)^  E  (a.u) 


50% 

-265.0296602077 

25% 

-265.2343008441 

15% 

-265.2982444057 

10% 

-265.3162851913 

7% 

-265.3200065394 

5% 

-265.3181039860 

4% 

-265.3164955496 

^See  Table  2. 
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TaUale  IV.  Electronic  State  Optimization  for  035^2^03”. 


State^  MO  Configuration^  Total  Energy 


10a 

4a2' 

If 

832 

12e' 

8C'' 

It 

2ai 

E  (a.g.) 

1 

2 

0 

1 

3 

0 

0 

-265.320006 

2 

2 

0 

2 

2 

0 

0 

-265.313010 

3 

0 

0 

2 

4 

0 

0 

-265.259282 

4 

1 

0 

2 

3 

0 

0 

-265.29367.2. 

5 

2 

0 

0 

4 

0 

0 

-265.309983 

^Average  energy  of  configuration.  Ref.  17. 
^Occupations  of  the  outermost  MOs'  are  shown. 
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Table  V.  Electronic  State  Optimization  for  CS]^^03. 

state^  MO  Confiqfuration^  Total  Energy 

lOaj^'  4a2'  Saj"  12e'  Se"  2aij^"  E  (a.u.J 

0  -26S. 296097 

0  -265.286243 

0  -265.287395 

0  -265.278357 

0  -265.264475 


^Average  energy  of  configuration.  Ref.  17. 
^Occupations  of  the  outermost  MOs'  are  sho%m. 
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TzUale  VI.  Electronic  State  optimization  of 

State^  MO  Configuration^  Total  Energy 

lOa^'  4a2'  Vaj”  8a2"  12e'  Be**  2ai"  E  (a.u.J 

2020200  -265.200801 

2021100  -265.209878 

2022000  -265.199121 

1  0  2  1  2  0  0  -265.198115. 

2010300  -265.130255 

^Average  energy  of  configuration.  Ref.  17. 

^Occupations  of  the  outermost  MOs'  are  shown. 


1 

2 

3 

4 
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Table  VII.  Population  Analysis  for  Cluster  Groiand  States. 


s-type  p-type  Total 


CS1103" 

OA 

5.825 

16.934 

22.759 

7.59 

-1.59 

CSA 

6.916 

18.764 

25.679 

8.56 

0.44 

CSB 

15.950 

36.856 

52.806 

8.80 

0.20 

CSC 

3.635 

13.120 

16.755 

8.38 

0.62 

CS1103 

OA 

5.825 

16.944 

22.769 

7.59 

-1.59 

CsA 

7.205 

18.671 

25.875 

8.62 

0.37 

CsB 

14.274 

36.503 

50.778 

8.46 

0.54 

CsC 

3.681 

12.896 

17.577 

8.79 

0.21 

CSiiOs'^ 

OA 

5.825 

16.949 

22.774 

7.59 

-1.59 

CsA 

7.214 

18.563 

25.777 

8.59 

0.41 

CsB 

13.528 

36.288 

49.815 

8.30 

0.70 

CSC 

4.849 

12.784 

17.633 

8.82 

0.18 

^Clustec, 

.■1 

^Type  of  atom  (See  Table  l  and  Fig.  I.) 
^Gross  atomic  population. 

‘^Charge  per  atom. 
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Table  VIII.  HOMO  and  LOMO  Energies  of  CS]^]^03,  CS11O3  and 

in  Their  Lowest-Energy  States 


Cluster  Orbital 

HOMO  (a.U.) 

UJMO  (a 

CS11O3  a^' 

-0.0436 

0.0492 

^2' 

-0.1285 

0.0977 

^2” 

-0.0377 

0.0140 

e' 

-0.0257 

0.0341 

e' 

-0.1339 

0.0566 

^1" 

-0.6163 

0.1236 

-0.1003 

-0.0103 

^2' 

-0.1964 

0.0334 

^2" 

-0.0989 

-0.0454 

e' 

-0.0877 

-0.0470 

e" 

-0.2022 

-0.0009 

^2" 

-0.6870 

0.0590 

-0.1587 

-0.0704 

■JtT^  ■ 

■4.  -  *■■ 

-0.2667 

-0.1615 

-0.0321 

-0.1057 

B* 

-0.1517 

-0.1294 

&” 

-0.2732 

-0.0598 

SB  ** 

-0.7590 

-0.0070 
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Figure  Captions 


Figure  1.  The  ^3^3^03  cluster  from  several  perspectives. 


Figure  2.  Total  valence  energy  (a.u.)  vs.  geometry  for  a 
CSii03“  cluster. 
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Ab  Initio  Stndy  of  Hydrogen  Adsorption  on  Be  (0001) 

M.  M.  Marino 

Department  of  Chemistry  and  Chemical  Engineering* 

Stevens  Institute  of  Technology 
Hoboten,  New  Jersey  07030 

and 

W.  C.  Ermler 
Department  of  Chemistry 
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C*"  Ci  ^  Ay  5.  V  ^  ^ ^ 

ABSTRACT 

Hydrogen  adsorption  on  the  (0001)  beryllium  surface  is  modelled  using 
a  Be^5  cluster  containing  seven  layers.  The  ab  initio  Hartree-Foct 
calculation*  «aploy  effective  core  potentials  and  full  D^^  point  group 
symmetry.  lowlying  electron  configurations  are  investigated  and  the 

effect  of  etbctron  spin  coupling  on  the  adsorption  process  is  discussed. 

The  adsorption  of  H  on  Be  is  calculated  to  be  stable  by  20  kcal/mol  and 
lowers  the  work  function  by  1.3  eV. 
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I .  INTRODUCTION 


The  use  of  metal  clusters  to  represent  bulk  surfaces  involved  in 
adsorption  phenomena  is  widespread.  A  good  model  of  the  adsorption 
process,  however,  requires  the  use  of  clusters  large  enough  to  describe  the 
metal  surface  and  its  interaction  with  the  adsorbate,  while  remaining  small 
enough  for  the  application  of  accurate  theoretical  treatments.  Beryllium 
has  been  widely  used  for  such  studies  because  it  is  the  smallest  aton  with 
a  closed-shell  ground  electronic  state  to  exhibit  metallic  character.^ 

The  toxicity  of  beryllium,  however,  has  limited  its  use  experimentally. 
Nevertheless,  two  recent  experimental  studies  have  investigated  tha 
adsorption  of  Hydrogen  on  beryllium  metal. The  former  study^^  was' 
concerned  with  H2  adsorption  on  polycrystalline  Be,  while  the  latter22 
involved  the  effects  of  H2,  H  and  D  on  the  (0001)  surface  of  beryllium 
metal.  The  experimental  data,  in  conjunction  with  theoretical 

calculations,  can  be  used  for  the  development  of  more  reliable  models  for 
adsorption. 

The  model  of  the  beryllium  surface  used  in  this  study  consists  of  a 
seven-layer  cylindrical  piece  from  a  wafer  having  Be-Be  internuclear 
distances  equal  to  those  in  the  bulk  metal.  Work  function  values,  the 
effect  of  api^  Mupling  on  the  cluster  and  net  charge  differences  between 
neutral  an4^||illl^eluaters  are  analyzed  to  determine  how  closely  the  model 
approximatelr  tmii  behavior.  Where  possible,  comparisons  between  theory  and 
experiment  are  made.  Finally,  adsorption  energies  are  reported  and  the 
effects  of  H  adsorption  on  the  Be  surface  are  discussed  in  the  context  of 
adsorbate  plus  substrate  as  a  supermolecule. 

II.  CALCULATIONS 
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The  model  system  is  taken  as  a  cylindrical  piece  from  a  Be  wafer  whose 
surface  corresponds  to  the  (0001)  hep  metal  face.  Table  I  defines  all  the 
cylinders  of  a  given  radius  which  may  be  formed  from  a  wafer  of  a  given 
thickness.  The  seven-layer  wafer  in  this  study  has  a  radius  w;.ich 

corresponds  to  a  one  unit  cell  step  along  the  a-direction.  Be-Be 
internuclear  distances  are  taken  as  those  in  the  bulk  hep  metal  (a«2.29A, 
c-3.58A).^^  The  resulting  cylinder  has  point-group  symmetry,  and  the 
adsorption  of  atomic  hydrogen,  one  on  the  top  surface  and  one  on  the  bottom 
surface  of  the  cluster  is  modeled  such  that  the  three-fold  symmetry  and 
horizontal  mirror  plane  are  preserved.  The  Be-H  internuclear  distances  are 
varied  in  a  series  of  self-consistent  field  (SCF)  calculations.  The 
hydrogenated  Be^^  cylinder  is  shown  in  Fig.  1.  Results  obtained  using  this 
system  are  compared  to  those  calculated  by  treating  an  identical  bare  Be^^ 
cylinder. 


Calculations  were  accomplished  on  a  Cray  Y-MP  supercomputer  using 

programs  based  on  the  "equal  contribution  theorem"  for  two-electron 
14 

integrals.  Ab  initio  restricted  closed-shell  and  restricted  open-shell 
Hartree-Fock  calculations  corresponding  to  an  average  energy  of 
configuration  as  well  as  to  the  energies  for  pure  spin  couplings  were 


carried  out  on  nusierous  low-lying  states  of  each  cluster.  (The  average 

'  ■- . 

energy  of  cQitifvration  is  defined  as  the  weighted  mean  of  the  energies  of 

■:4-  ,  15 

all  the  oitfltiplets  for  the  configuration.  )  Multiconfiguration 


self-consistent  field  (MCSCF)  calculations  involving  two  configurations 


were  also  carried  out  for  selected  electronic  states.  The  calculations 


required  about  one  hour  of  Cray  Y-MP  time  for  each  geometric  orientation. 
The  following  basis  seta  of  contracted  Gaussian-type  functions  were  used 
for  beryllium  and  hydrogen,  respectively:  (3s2p) / (2alp] and 
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{4alp)/(28lp).  Ab  initio  effective  potentials  (CP)  were  used  to 

represent  the  Is  core  electrons  in  Be.^®  (Valence  SCF  energies  for  atonic 

Be  and  H  are  -0.95083,  -0.49928  hartrees,  respectively,  for  the  basis  sets 

used  in  this  study.) 

III.  RESULTS 


Total  valence  energies  for  six  low-lying  electronic  states  are  given 
in  Table  II.  These  results  are  reported  for  both  the  equilibrium 
hydrogen-to-Be-plane  distance  and  at  the  near  dissociative  limit  of  lOA. 

All  of  the  energies  correspond  to  triplet  spin  couplings  with  the  exception  | 

■  i 

of  state  2,  which  is  the  lowest-lying  closed  shell  singlet  state  fOMA.  > 
Table  III  shows  the  effects  of  electron  correlation,  proper  dissociating 
wavefunctions  and  spin  coupling  for  the  lowest-lying  electron  configuration 
calculated  at  the  minimum  0.86A  and  at  the 

dissociative  limit  of  lOA.  Table  IV  contains  energy  values  at  dissociation 
for  state  6,  a  low-lying  MO  configuration  having  the  largest  number  of 
open-shell  electrons. 

Total  SCF  energies  for  the  ground  state  of  H  and  of  each  cluster,  as 
well  as  calcplated  adsorption  energies  are  given  in  Table  V.  Theoretical 
work  func^Hk^N^  shown  in  Table  VI.  The  first  values  in  each  row  were 
calculatsdiS^Hfl  diffsrence  between  the  total  valence  SCF  energy  of  the 
neutral  And  <  that  of  the  ion  resulting  from  the  resioval  of  one 

electron  from  the  highest  occupied  orbital  of  the  cluster.  The  values 
immediately  following  are  due  to  Koopmans'  theorem  and  correspond  to  the 
negative  of  the  energy  of  the  molecular  orbital  (MO)  frcm  which  the 
electron  was  ionized.  Experimental  values  are  also  shown. 

Atonic  net  electron  charge  values  are  given  in  Table  VII,  as 
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19 

calculated  using  a  Mulliken  population  analysis.  Values  of  R  in  this 
table  correspond  to  the  optimum  H  atom-to-surface  distance  for  the 
lowest-lying  electronic  state.  This  distance  was  calculated  as  the  minimum 
of  a  curve  fit  to  SCF  energies  for  three  or  more  surface-adsorbate 
separations. 

IV.  DISCUSSION 


Only  two  of  the  six  low-lying  electron  configurations  found,  states  1 
and  6  of  Table  II,  dissociate  properly  at  the  restricted  Hartree-FocJt  level 
of  approximation.  The  remainder  of  the  states  lie  significantly  higher  in 
energy  at  dissociation.  An  analysis  of  their  Mulliken  populations 
indicates  that  the  net  charges  on  the  H  atoms  at  lOA  separation  (R^)  are  as 
high  as  0.5  electron/atom,  clearly  indicative  of  an  improper  dissociation 
liadt . 

All  of  the  open-shell  states  appearing  in  Table  II  are  coupled  as 
triplets.  Tables  III  and  IV  show  the  effects  of  spin  coupling  on  the 
energy  of  the  Be45H2  system.  Both  the  two  configuration  MCSCF  and  the  pure 
spin  triplet  calculations  (Table  III)  yield  the  same  energy  at  dissociation 
(-45.57065  hartrees) .  However,  the  improperly  dissociating  singlet  spin 
state  is  significantly  higher  in  energy  (144  kcal/raol  higher) .  The  average 
energy  of  configuration  calculation  also  yields  a  much  higher  energy  (by  94 
kcal/mol) ,  but  this  calculation,  as  mentioned  previously,  reflects  a 
weighted  average  of  the  energies  for  ail  the  multiplets  of  a  particular 
configuration.^^  Therefore,  the  higher  energy  is  a  consequence  of  the 
contribution  from  the  singlet  component  of  the  multiplet.  These  results 
indicata  that  the  unpaired  electrons  in  the  overall  cluster  must  be  coupled 
as  triplets  with  respect  to  each  other  to  insure  proper  dissociation  in  a 
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single  configura^on  treatment .  According  to  Mulliken  population  analyses, 
these  electrons  correspond  solely  to  hydrogen  Is  electrons  at  R^.  Although 
the  energies  are  not  identical,  the  MCSCF  result  also  agrees  well  with  that 
of  the  pure  spin  triplet  case  at  equilibrium,  the  discrepency  between  the 
two  being  approximately  6  kcal/mol,  which  is  attributed  to  electron 
correlation  effects.  However,  even  at  this  short  distance,  coupling  the 
unpaired  electrons  singlet  as  opposed  to  triplet  results  in  a  20  kcal/mol 
increase  in  energy. 

The  effects  of  spin  coupling  on  the  Be45H2  system  at  dissociation  were 
further  studied  using  state  6  of  Table  II.  This  is  the  highest  spin 
electron  configuration  (quintet)  found  to  dissociate  properly.  The 
energies  for  various  spin  states  of  this  configuration  were  calculated. 
These  results  appear  in  Table  IV.  Three  states,  the  quintet,  triplet  and 
first  singlet,  have  energies  which  are  nearly  identical.  The  second 
singlet,  however,  is  160  kcal/mol  higher  in  energy.  Again,  the  average 
energy  of  configuration  calculation  reflects  the  contribution  of  the  second 
singlet  state.  As  with  the  previous  results,  the  three  low-lying  states  in 
Table  IV  have  the  two  unpaired  electrons  corresponding  to  the  hydrogen  Is 
electrons  coupled  triplet  with  respect  to  each  other,  while  the  less  stable 
singlet  state  has  these  electrons  coupled  singlet.  The  coupling  of  the 
unpaired  electrons  associated  primarily  with  Be  orbitals  does  not  affect 
the  energy  significantly.  For  example,  the  triplet  spin  state  and  the  more 
stable  singlet  spin  state  have  the  unpaired  Be  electrons  coupled  singlet 
and  triplet,  respectively,  yet  the  energy  difference  between  the  two  is 
only  2.4  Kcal/mol.  We  are  encouraged  by  this  behavior  since  it  is  expected 
if  the  Be^g  cluster  provides  an  accurate  model  of  the  bulk  surface,  which 
is  obviously  not  affected  by  spin  considerations. 

The  total  adsorption  energy  of  Be^5H2  is  19.72  kcal/mol -atom  (Table 
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V) .  This  energy  was  determined  two  ways.  First,  the  adsorption  energy  was 
calculated  as  the  total  energy  difference  between  the  lowest-lying 
electronic  state  found  at  equilibrium  and  at  dissociation.  Second,  the 
adsorption  energy  was  defined  as  the  difference  between  the  energy  of  the 
lowest-lying  electronic  state  found  at  equilibrium  and  the  sum  of  the 
energies  for  the  lowest-lying  electronic  state  found  for  bare  Be^^  plus  two 
ground  state  H  atoms.  Both  methods  yielded  values  of  19.72  kcal/mol-atom. 
Fig.  2  shows  a  plot  of  energy  vs  H  to  Be  surface  distance.  A  barrier  of 
24.68  kcal/mol  (12.34  kcal/mol-H  atom)  occurs  before  H  adsorbs  to  the  Be^^ 
surface. 

A  very  recent  experimental  study  by  Ray  et  al.  investigated  the 

effects  of  H  and  H2  adsorption  on  the  (0001)  surface  of  Be  metal  using  high 

resolution  electron  energy  loss  spectroscopy  (HREELS),  thermal  desorption 

12 

spectroscopy  (TDS)  and  workf unction  measurements.  This  group  found  no 
evidence  for  the  dissociative  or  molecular  chemisorption  of  H2  at  the 
temperatures  investigated,  in  agreement  with  previous  experimental  work.^^ 
They  also  found,  however,  that  H  does  adsorb  on  Be  and  reported  the  energy 
barrier  to  recombinative  desorption,  obtained  using  TDS,  as  ranging  from 
0.7  to  1.2  eV  (16.1  to  27.7  kcal) ,  depending  on  the  degree  of  H  coverage. 
This  barrier  energy  is  equal  to  the  sum  of  the  adsorption  site  binding 
energy  and  the  barrier  to  dissociative  adsorption.  This  value  is  in 
excellent  agreement  with  our  calculated  adsorption  energy  of  0.85  eV  or  20 
kcal/mol’atOB  (Table  V). 

The  convergence  of  both  methods  to  the  same  value  for  the  adsorption 
energy  indicates  that  the  lowest-lying  electronic  state  correlating  to  the 
dissociative  limit  has  been  determined  accurately.  Such  a  determination  is 
pivotal  because  properties  such  as  work  functions,  adsorption  energies  and 
charge  distributions,  which  indicate  the  chemical  nature  of  the  system,  are 
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sensitive  to  the  electronic  state  chosen.  In  our  previous  work,^^  we 

concluded  that  H  does  not  chemisorb  to  Be^^  since  no  shift  in  the 

workf unction  of  the  Be^^  surface  occurred  upon  H  adsorption.  Furthermore, 

the  net  charge  distribution  difference  between  Be^^^  vs.  Be^^  and  86^5112^ 

vs.  Be^2H2  were  nearly  identical,  indicating  that  the  ionization  process 

was  not  affected  by  the  presence  of  hydrogen.  In  the  present  calculations, 

however,  a  lower-lying  electronic  state  has  been  found  for  both  Be^^  and 

®®45^2'  Although  the  adsorption  energy  corresponding  to  this  state  is  only 

1.56  kcal/mol'atom  lower  than  was  determined  previously,  the  workfunction 

of  Be^^  is  now  4.4  eV  (Table  VI),  in  close  agreement  with  the  experimental 
20 

value  of  3.92  eV.  Adsorption  of  H  onto  Be^^  lowers  the  work  function  to 

3.14  eV,  a  shift  of  1.29  eV.  No  work  function  shift  is  observed 

experimentally  for  the  H2/Be  system. However,  a  direct  comparison 

between  the  experimental  and  theoretical  work  function  shifts  in  this  case 

cannot  be  made  since  the  experiment  involved  sputtered  H2  and 

polycrystalline  Be.  In  such  a  process,  H  adsorption  is  necessarily 

preceded  by  H2  dissociation,  which  requires  approximately  104  kcal/mol  at 
21 

room  temperature.  This  is  considerably  more  than  the  calculated  20 

kcal/mol  bond  formed  between  hydrogen  and  the  Be^^  surface.  Consequently, 
it  is  more  energetically  favorable  for  hydrogen  to  remain  as  H2,  unless,  of 
course,  the  surface  plays  a  heretofore  unclear  catalytic  role  in  lowering 
the  energy  needed  to  break  the  H-H  bond.  In  either  case,  the  focus  of  the 
present  calculations  is  to  study  the  interaction  of  H  with  the  Be  surface 
and  not  the  dissociation  of  H2.  For  the  H/Be  (0001)  system,  however,  the 
work  function  of  Be  reaches  a  minimum  (-0.56  eV)  for  a  coverage  of 
approximately  0.4  and  then  reaches  a  limiting  value  of  -0.44  eV  as  the  H 
coverage  is  increased.  This  is  also  in  good  agreement  with  our  calculated 
work  function  lowering  of  1.3  eV  (0.65  eV/atom) .  It  is  pointed  out, 
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however,  that  Ray  et  al.  believe  their  values  correspond  to  adsorption  of 
H  on  a  bridge  site,  whereas  the  present  study  investigates  H  adsorption  on 
one  type  of  trigonal  site  (eclipsed  site) . 

Net  charge  differences  appearing  in  Table  VII  indicate  that  for  both 
Be^^  and  Be^^H2,  the  largest  fraction  of  the  donated  electron  comes  from 
MO'S  involving  the  surface  atoms.  This  behavior  is  requisite  if  the  Be^^ 
cluster  is  to  provide  an  adequate  representation  of  the  bulk  metal  since 
the  wor)e  function  is  a  surface  property.  Adsorption  of  H  increases  the 
ability  of  the  surface  to  release  an  electron  to  some  degree  (0.48  vs  0.60 
for  the  bare  cluster  vs  the  hydrogenated  cluster,  repectively) .  However, 
there  is  no  net  charge  transfer  between  the  adsorbed  hydrogen  and  the 
surface,  indicating  an  essentially  covalent  interaction  between  the  two. 

In  this  study,  H  approaches  the  Be  surface  directly  above  the  center 
of  a  triangle  of  Be  atoms  located  on  the  surface,  the  center  of  which  is 
located  above  a  Be  atom  situated  in  a  layer  next  to  the  surface  (Fig.  1) . 
This  is  called  the  eclipsed  site.  There  are  three  other  possible  sites  for 
adsorption.  These  are  the  Be-Be  midpoint,  directly  overhead  and  open 
sites.  The  first  has  the  adsorbate  approaching  the  Be  surface  directly 
between  two  Be  atoms  located  on  the  surface.  The  second  Involves 
adsorption  directly  above  a  Be  atom  located  on  the  surface.  The  open  site 
has  the  adaocbtta  approaching  the  surface  directly  above  a  triangle  of  Be 
atoms  located  iili  the  surface,  directly  below  which  there  is  no  Be  atom. 

The  adsorption  of  H  onto  Be  clusters  modelling  the  bulk  surface  has 
also  been  studied  at  the  SCF  level  of  theory  by  other  groups.  In  the 
following  discussion  the  notation  Bej^(a,b,c, . . .)  will  be  used.  Here  n 
refers  to  the  number  of  atoms  comprising  the  Be  cluster  and  a,  b,  c,  etc., 
ate  the  number  of  atoms  in  the  surface  layer,  the  layer  directly  beneath 
the  surface,  the  second  layer  beneath  the  surface,  and  so  on,  respectively. 
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In  all  cases,  unless  otherwise  stated,  the  cluster  interatomic  separations 
were  frozen  at  the  bulk  metal  values. 

Q 

Cox  and  Bauschlicher  studied  the  adsorption  of  H  in  the  directly 
overhead  position  on  a  Be.^(7,0)  and  Be^Q(7,3)  cluster.  Adsorption  of  H  on 
a  bond  midpoint  site  was  modelled  using  a  Be^2(10,5)  cluster.  They 
obtained  adsorption  energies  of  65.5,  21.22  and  53.7  kcal/mol  for  the  three 
above-mentioned  sites,  respectively.  They  further  conclude  that  in  two 
previous  studies  of  H  adsorption  on  Be,  one  using  clusters  or  up  to  22 
atoms  and  one  using  a  Be^^  cluster  ,  the  chemisorption  energy  changed  by 
approximately  10%  in  going  from  a  22-atom  cluster  to  a  36-atom  cluster  and 
by  about  30%  in  going  from  a  Be^^  or  Be^^  cluster  to  a  Be^^  cluster. 

In  another  study,  Bauschlicher  and  Bagus^  reported  adsorption  energies 
of  64.0,  47.8,  and  35.8  kcal/mol  for  H  adsorption  on  the  directly  overhead 
site  of  Be.^(7,0),  the  eclipsed  site  of  Be.^(6,l)  and  the  bond  midpoint  site 
of  Be^Q(10,0),  respectively.  Calculations  by  Bagus  et  al.^  on  the 
Bej^(14,0)H  system  found  the  directly  overhead  site,  with  a  desorption 
energy  of  59.0  kcal/mol,  to  be  a  more  stable  site  for  H  adsorption  than  the 
bond  midpoint,  open  or  eclipsed  sites  (53.1,  56.1  and  56.1  kcal/mol, 
respectively).  H  adsorption  on  Be22(14,8),  on  the  other  hand,  is  predicted 
to  be  more  stable  for  the  open  site  (55.1  kcal/mol)  than  for  the  bond 
midpoint,  eclipsed  or  directly  overhead  sites  (53.4,  51.9,  and  31.4 
kcal/mol,  re^ectively) .  Finally,  the  most  stable  site  for  H  adsorption  on 
a  Be^g(14,8,14)  cluster  was  predicted  to  be  the  open  site,  also  (57.7 
kcal/mol) ,  with  the  bond  mi^oint,  eclipsed  and  directly  overhead  sites 
(43.7,  42.3  and  31.6  kcal/mol,  respectively)  giving  rise  to  lower 

adsorption  energies.  Clearly,  Be^g  and  Be22  9ive  rise  to  the  same  ordering 
in  the  stability  of  the  adsorption  energy  with  respect  to  bonding  site. 
However,  with  the  exceptions  of  the  directly  overhead  and  open  sites,  the 
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predicted  adsorption  energies  for  a  given  site  on  these  two  clusters  differ 
by  as  much  as  9.7  kcal/inol. 

The  results  reported  here,  along  with  those  cited  from  previous 
studies,  indicate  the  importance  of  cluster  choice  in  modelling  the  Be 
metal  surface  involved  in  H  adsorption.  It  must  be  ascertained  that  the 
cluster  model  adequately  represents  the  bulk  surface,  or  the  calculated 
results  will  correspond  to  H  adsorption  on  a  specific  electronic  state  of  a 
cluster  rather  than  on  the  (0001)  Be  surface.  This  explains  the 
discrepancies  in  adsorption  energies  among  the  various  studies  cited.  We 
believe  that  the  energies  arising  from  calculations  on  small  clusters 
correspond  to  H  adsorption  on  clusters  rather  than  on  the  Be  surface. 

There  are  two  critical  steps  in  establishing  a  valid  model  of  H 
adsorption  on  a  surface.  First,  the  cluster  representing  the  surface  must 
be  large  enough  to  model  the  bulk  solid  accurately.  Second,  the  ground 
electronic  state  of  the  system  must  be  determined.  A  recent  study^^  of  H 
adsorption  on  the  Be  (0001)  surface  investigated  three  possible  cluster 
models  of  the  bulk  surface:  Bej^g(6,7, 6) ,  ,S,1)  and 

Be45 (6, 7, 6, 7, 6, 7, 6) .  The  first  and  third  modelled  adsorption  on  an 

eclipsed  site,  while  the  second  involved  adsorption  on  the  directly 
overhead  site.  Since  the  adsorption  site  in  Be^^  is  identical  to  that  in 
Be^j,  results  for  both  systems  should  agree.  However,  the  calculated  shift 
in  the  ionization  potential  of  Be^^^  due  to  H  adsorption  was  0.0  eV,  while 
that  corresponding  to  H  adsorption  on  Be^^  was  1.4  eV.  The  work  functions 
of  Be^j  were  calculated  as  5.0  and  4.4  eV,  with  the  experimental 

value  being  3.9  eV.  These  discrepancies  were  explained  using  total  net 
charge  differences  between  neutral  and  ionized  clusters.  In  the  case  of 
Be^g,  ionization  of  the  cluster  indicated  participation  in  charge 
redistribution  by  both  the  surface  and  middle  layers.  Since  electron 


99 


emission  occurs  from  the  surface  of  the  bulk  metal  and  should  not  involve 

the  middle  layers,  it  was  concluded  that  Be^g  is  too  small  to  model  the  Be 

(0001)  surface.  As  previously  stated,  net  charge  differences  for  Be^^  (and 

indicate  that  the  largest  fraction  of  the  donated  electron  comes 

from  MO's  involving  the  surface  atoms  (see  Table  VII).  Since  the  work 

function  is  a  surface  property,  this  must  be  the  case  if  the  Be^^  cluster 

adequately  models  the  bulk.  Furthermore,  the  calculated  work  function  of 

Be^j  is  only  0.5  eV  greater  than  the  experimental  value  of  3.92  eV.^^ 

.  23 

A  study  by  Pacchioni  et  al.  found  the  directly  overhead  site  on  a 
Be.^(7,0)  cluster  to  be  more  stable  (48.3  kcal/mol)  than  the  bond  midpoint 
or  open  sites  (35.7  and  11.7  kcal/mol,  repsectively) .  They  further  point 
out  that  this  finding  holds  only  for  single-layer  models  of  the  bulk  metal. 
As  more  layers  are  are  added,  the  bond-midpoint  and  open  sites  become  more 
favored  than  the  directly  overhead  site.^'^  This  is  in  agreement  with  our 
Be^^  (  directly  overhead  site)  and  Be^^ (eclipsed  site)  models  of 
adsorption.  While  H  is  not  predicted  to  adsorb  on  Be^^,^^  it  does  adsorb 
on  Be^^  by  19.7  kcal/mol  (see  Table  V).  Both  clusters  are  multilayered 
cutoffs  of  the  bulk,  with  Be^^  containing  5  layers  and  Be^^  containing  7 
layers . 

As  pointed  out  earlier,  the  Be^^  cluster  in  this  work  contains  Be-Be 
interatomic  distances  equivalent  to  the  bulk  metal  values,  and  the  cluster 
is,  in  effect,  a  cylindrical  "plug"  from  the  bulk  metal  (Table  I) .  This 

approximation  is  valid  since  experimental  evidence  indicates  that  the  Be 

24  8 

surface  does  not  reconstruct.  Furthermore,  Cox  and  Bauschlicher  studied 

the  effects  of  relaxation  on  Be.^(7,0),  BejQ(7,3)  and  Bej^5(10,5)  clusters 

due  to  H  adsorption.  Adsorption  energies  due  to  relaxation  effects  were 

calculated  as  66.9,  21.4  and  55.6  kcal/mol,  respectively,  while  the 

energies  arising  from  adsorption  on  an  unrelaxed  substrate  were  reported  as 
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65.5,  21.2  and  53.7  kcal/mol,  respectively  (refer  to  earlier  discussion). 

They  predict  that  neglect  of  relaxation  effects  results  in  only  a  2 

kcal/mol  error  in  the  calculated  binding  energies. 

The  final  important  point  to  consider  when  modelling  H  adsorption  on  a 

metal  surface  is  the  method  of  determining  the  ground  electronic  state  of 

the  system.  This  is  imperative  since  the  chemical  nature  of  the  system  is, 

in  effect,  mainly  determined  by  the  ground  electronic  state.  We  have 

already  discussed  the  errors  resulting  from  improperly  predicting  the 

ground  electronic  state  of  Be^j  and  Be^jH2  in  out  previous  work.  Another 

example  of  the  ambiguities  or  problems  which  may  arise  from  the  choice  of 

23 

ground  electronic  state  can  be  seen  in  the  studies  of  Paccioni  et  al.  and 

p 

Cox  and  Bauschlicher.  The  latter  reported  relaxation  effects  in  the 

geometry  of  Be  clusters  due  to  H  adsorj^tion  (refer  to  prior  discussion) . 

However,  this  group  chose  the  ground  state  of  the  Be.^  system  to  be  closed 

shell  in  order  to  avoid  the  effects  of  dangling  orbitals  on  the  energy  of 

23 

the  system,  Pacchioni  et  al.  attribute  this  shrinking  to  the  choice  of 

electronic  state.  This  group  found  that  the  relaxation  effects  occurred  in 

different  directions  for  different  electronic  states  and  that  they  also 

depended  on  the  choice  of  theoretical  method  (all  electron  vs 

pseudopotential  at  both  the  SCF  and  Cl  levels  of  theory) .  The  effects  of 

relaxation  on  the  adsorption  energy,  they  suggested,  should  be  investigated 

further  with  other  cluster  models. 

25 

Panas  et  al.  studied  the  cluster  convergence  of  chemisorption 

energies.  This  group  generated  molecular  orbitals  using  SCF  calculations 

and  included  correlation  corrections  via  a  one  reference  state  contracted 
Cl  procedure.  They  studied  H  adsorption  in  the  fourfold  hollow  position  of 
Ni(lOO)  using  seven  Ni  clusters  ranging  in  sire  from  Nij(4,l)  to 
Ni^Q (16,9,16,9)  to  represent  the  Ni  surface.  They  found  a  marked 
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difference  in  the  calculated  chemisorption  energies  of  the  clusters 
depending  on  the  method  used  to  determine  the  lowest-lying  electronic  state 
at  short  and  long  distances  (the  chemisorption  energy  being  defined  as  the 
defference  between  the  energies  of  these  two  states) .  This  group  concludes 
that  the  determination  of  the  chemisorption  energy  is  greatly  improved  by 
requiring  that  the  cluster  systems  first  be  excited  into  a  proper  bonding 
state  and  then  defining  the  chemisorption  energy  as  the  difference  in  the 
energy  corresponding  to  this  state  at  short  and  long  distances  plus  the 
energy  needed  for  excitation.  In  the  present  study,  no  restrictions  were 
imposed  on  the  symmetry  and  spin  of  the  wavefunction  at  short  and  long 
distances.  However,  applying  the  rules  outlined  by  Panas  et  al.  led  to  the 
same  value  for  the  adsorption  energy.  Thus  the  energy  value  reported  here 
is  expected  to  be  accurate. 

One  last  point  must  be  mentioned.  All  of  the  results  reported  in  the 

current  study  were  calculated  at  the  SCF  level  of  approximation.  It  is 

therefore  necessary  to  estimate  the  extent  to  which  electron  correlation 

26 

corrections  will  affect  the  current  results.  Rubio  et  al.  studied  the 

effect  of  electron  correlation  on  the  process  whereby  H  adsorbs  to  small 

(three  to  seven  atoms)  Be  clusters.  They  determined,  using  data  obtained 

at  both  the  SCF  snd  Cl  levels,  that  SCF  theory  is  qualitatively  accurate  to 

describe  the  interaction  of  H  with  the  clusters.  The  changes  in  the 

H-to-Be-surface  distances  and  in  the  calculated  adsorption  energies  were 

27 

quite  small,  A  study  by  Siegbahn  et  al.  of  H2  dissociation  on  the 
Ni(lOO)  surface  using  Ni  clusters  with  as  many  as  34  atoms  to  represent  the 
surface  also  determined  that  the  SCF  treatment  was  adequate  for  a 
qualitative  description  of  the  chemisorption  process.  Given  these  findings 
and  the  fact  that  the  Be45H2  system  is  quite  large  (92  valence  electrons) 
the  incranental  effects  of  electron  correlation  on  the  results  reported 
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here  are  expected  to  be  small. 


V.  SUMMARY  AND  CONCLUSIONS 

An  accurate  model  of  the  adsorption  process  requires  that  the  surface 

be  modelled  by  a  cluster  which  provides  an  adequate  representation  of  the 

bulk  structure.  Once  this  is  done,  it  is  equally  important  to  determine 

accurately  the  lowest-lying  electronic  configuration  of  the 

cluster-adsorbate  system.  We  conclude  that  provides  an  accurate  model 

of  the  Be  (0001)  surface  for  several  reasons.  First,  this  cluster  is  not 

affected  by  spin  considerations,  as  is  expected  for  solid  surfaces.  Spin 

coupling  affects  the  total  energy  only  when  electrons  corresponding  to 

hydrogen  orbitals  are  involved.  Second,  the  ionization  potential  of  Be^^ 

at  4.4  eV  is  close  to  the  value  of  3.9  eV  for  bulk  Be  metal.  Third, 

Mulliken  population  analyses  indicate  that  the  electron  emitted  during 

ionization  is  donated  from  the  surface  layers,  which  is  also  the  case  for 

bulk  metal  surfaces.  Fourth,  previous  studies  have  indicated  the  need  for 
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clusters  several  layers  thick  to  accurately  model  H  adsorption.  '  ' 

Be^j  is  a  seven-layer  cluster  based  on  hep  bulk  lattice  constants. 

The  lowest-lying. electronic  state  of  the  Be^jH2  system  has  been  determined. 
An  indication  of  the  reliability  of  this  determination  is  the  fact  that  the 
workfunctiona  of  the  cluster  and  the  bulk  metal,  at  4.4  and  3.9  eV, 
respectively,  differ  by  only  0.5  eV.  In  our  previous  study^^,  the 
improperly  chosen  ground  electronic  state  caused  the  predicted  cluster 
wotkf unction  to  be  1.25  eV  lower  than  the  current  value  and  0.77eV  lower 
than  that  of  bulk  Be.  A  second  indication  that  the  ground  state  is  correct 
13  the  convergence  of  the  calculated  adsorption  energies  to  the  same  value 
using  two  different  methods.  This  is  expected  if  the  true  ground 
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electronic  state  has  been  found  for  both  systems.  In  the  present  case,  the 
systems  were  Be^^  +  2H  and  Be^5H2  at  infinite  Be-H  separation.  Finally, 
the  agreement  between  the  MCSCF  calculations  and  the  calculations  on  the 
triplet  spin  states  indicates  that  the  bulk  of  the  differential  correlation 
energy  is  accounted  for  by  requiring  the  electrons  corresponding  to  H  to  be 
coupled  as  triplets.  The  resulting  electronic  states  can  thus  be  described 
adequately  using  single  configuration  SCF-MO  calculations. 

Finally,  the  results  for  H  adsorption  on  Be^^  indicate  that  H  adsorbs 
at  an  eclipsed  site  of  Be  (0001)  by  20  kcal/mol  and  that  the  adsorption 
process  lowers  the  workfunction  of  the  cluster  from  4.4  to  3.1  eV,  or  1.3 
eV,  Since  no  charge  transfer  is  observed  between  hydrogen  and  the  Be 
surface,  the  adsorption  process  involves  a  primarily  covalent  interaction 
between  the  two. 
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Table  I.  Be  Clusters  by  Coordination  Cylinder  and  hep  Layer. 
Cylinder 

Height  R  R-=a  R. 

(z-coord.)  i  j 

o»c/2 


3c/2 

c 

c/2 

0 

-  c/2 

-  c 
-3c/2 


-••c/2 

No.  Atoms  1133  1  7  9  15  7  13  27  31  7  19  33  45 


0  3  3  6 

11  11  77  77 

000  333  333  666 

1111  1111  7777  7777 

000  333  333  666 

11  11  77  77 

0  3  3  6 
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Table  II.  Restricted  Hartree-Fock  Adsorption  Energies  of  States  of 


State 

MO  Configuration^ 

Energy^ 

No. 

a. ' 

a-' 

So" 

e' 

e" 

A  " 

R 

R 

1 

2 

2 

eq 

1 

1 

2 

1 

4 

4 

2 

-45.601315 

-45.570651 

2 

2 

2 

2 

4 

4 

2 

-45.604328 

-45.294292 

3 

1 

2 

2 

1 

4 

2 

-45.625789 

-45.264307 

4 

2 

2 

2 

2 

4 

2 

-45.633513 

-45.124689 

5 

2 

2 

2 

1 

1 

2 

-45.613706 

-45.101965 

6 

1 

1 

1 

1 

4 

2 

-45.591720 

-45.502133 

^All  open-shell  states  are  triplet  coupled.  Occupations  correspond  to  the 
number  of  electrons  occupying  the  highest  energy  molecular  orbital  of  a 
given  symmetry. 

and  R^  correspond  to  ^“^Sgurface  separations  of  0.86A  and  lOA, 
respectively.  Energies  are  in  hartrees. 
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Energies  (hartrees)  correspond  to  State  1  of  Table  II.  and 

correspond  to  ^"^Sgujface  ^^P^^^tions  of  0.86A  and  lOA, 
respectively,  ^^adsorb  difference  between  energies  at 

these  two  distances  (kcal/moi*H  atom) . 

^Two-configuration  MCSCF  (see  text) . 

'Average  energy  of  configuration  (see  text) . 
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Table  IV.  Restricted  Hartree-Fock  Energies  of  Be45H2  at  Dissociation.^ 


a 

Average  Energy 

-45.35865 

Pure  Spin  Singlet^ 

-45.50602 

Pure  Spin  Singlet^ 

-45.25044 

Pure  Spin  Triplet^ 

-45.50213 

Pure  Spin  Quintet 

-45.50603 

^Energies  (hartrees)  correspond  to  State  6  of  Table  II.  H  to  Be 


separation  is  lOA. 

^Electrons  in  symmetry  orbitals  a^^'  and  32"  ate  coupled  triplet  with 
respect  to  each  other  and  singlet  with  respect  to  those  in  symmetry 
orbitals  a2'  and  e'  (See  Table  II) . 

^Electrons  in  symmetry  orbitals  a^'  and  e'  are  coupled  triplet  with 
respect  to  each  other  and  singlet  with  respect  to  those  in  symmetry 
orbitals  a2  and  a2”.  (See  Table  II). 

'^he  electron  in  symmetry  orbital  e'  is  coupled  singlet  to  all  the 
other  open  shell  electrons,  which  are  coupled  triplet  to  each  other  (See 
Table  II) . 


no 
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Table  V.  Energies  o£  the  Be^jH2  System. 


Species^ 

Energy*^ 

h; 

-  0.49928 

B«,5 

-44.57210 

®«45«2 

-45.63351 

-45.57065 

®®45«2  -  B«45«2 

19.72 

®'«5"2  'V  -  I®*45  *  2"' 

19.72 

^Rgq*0.86A,  Rj^-IOA,  where  R  is  the 
Distance  from  H  to  the  Be  surface. 

^Flrst  four  values  are  in  hartrees,  while 
the  last  two  are  in  kcal/mol-atom.  All 
values  correspond  to  the  lowest-lying 
state  found  at  that  distance.  All  States 
are  triplet  coupled. 
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Table  VI.  Work  Functions  of  Be^^  and  8^^112  Clusters. 


Cluster 

0(ev)  ® 

®®45 

4.62,4.43 

®«45«2  , 

3.39,3.14 

Be  (expt.)° 

3.92 

BeH  (expt.)^ 

3.92 

A«(Be45-Be45H2]^ 

1.23,1.29 

0 

Sfhere  two  values  appear,  these  are  Koopmans' 
theorem  and  values,  respectively. 

’^Ref.  19. 

^Ref.  11. 

Values  correspond  to  a  work  function 
lowering. 
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Table  VII.  Atomic  Populations  of  and 


Cluster. Z-coprd^ 

Atom 

Labe*^ 

f 

No. 

Atoms'^ 

Cluster 

Net  Charge 

per  Atom 

{Cluster]^ 
Net  Charge 

per  Atom 

Total 

Net  Charge  Difference 
((Cluster]^  -  Cluster) 

Be, 5  0 

BeO 

1 

1.01 

1.00 

-0.01 

c/2 

BeA 

6 

0.42 

0.40 

-0.12 

c/2 

BeC 

6 

-0.20 

-0.15 

0.30 

0 

BeB 

6 

-0.33 

-0.31 

0.12 

c 

BeH 

12 

-0.19 

-0.17 

0.24 

c 

BeD 

2 

1.36 

1.34 

-0.04 

3c/2 

BeF 

6 

-0.14 

-0.14 

0.00 

3c/2 

BeG 

6 

0.02 

0.09 

0.48 

0 

BeO 

1 

0.99 

0.99 

0.00 

c/2 

BeA 

6 

0.43 

0.41 

-0.12 

c/2 

BeC 

6 

-0.16 

-0.13 

0.18 

0 

BeB 

6 

-0.36 

-0.32 

0.24 

c 

BeH 

12 

-0.15 

-0.14 

-0.12 

c 

BeO 

2 

1.31 

1.29 

-0.04 

3c/2 

BeF 

6 

-0.18 

-0.17 

0.06 

3c/2 

BeG 

6 

-0.04 

0.06 

0.60 

H 

2 

-0.01 

-0.01 

0.00 

^c=3.58A 
^See  Fig.  2. 

'No.  of  symmetry  equivalent  atoms 

^R=»0.86A  for  Be  plane  to  H  distance.  Be  to  H  distance  is  1.58A. 
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FIGURE  CAPTIONS 


The  potential  energy  curve  for  the  approach  of  H  to 
a  rigid  Be^^  surface  having  internuclear  separations 
corresponding  to  bulk  metal  lattice  constants.  Total 
energies  may  be  obtained  by  adding  -45.00  h  to  the 
given  values. 
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